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Численное моделирование поведения кинематически нестабильных 
склонов при динамических воздействиях 





П. П. Гайджуров® >4, Н.А. Савельева ©. Е. В. Труфанова © 


ФГБОУ ВО «Донской государственный технический университет» (г. Ростов-на-Дону, Российская Федерация) 


24 эрр-161@уапдех.га 


Введение. Предложена концепция оценки динамических параметров системы «основание — ослабленный слой — 
блок» с учетом физической нелинейности материала при кинематическом способе возбуждения колебаний. В 
соответствии с данным подходом учет физической нелинейности материала основания и блока осуществляется с 
помощью модели Друкера-Прагера. Ослабленный слой моделируется ЗО пружинными конечными элементами. На 
примере динамического расчета системы «основание — ослабленный слой — склон» осуществлена процедура 
верификации предлагаемой методики. 

Материалы и методы. Вычислительные эксперименты выполнены с помощью программного комплекса АМЗУЗ 
Месватса| в сочетании с нелинейным решателем, базирующемся на процедуре Ньютона-Рафсона. Для 
дискретизации расчетных областей применены объемные конечные элементы $045. Для моделирования 
смещения блока относительно неподвижного основания использованы комбинированные упруго-вязкие элементы 
СОМВИ\ 14. 

Результаты исследования. Разработана инженерная методика динамического анализа напряженно- 
деформированного состояния пространственной системы «основание — ослабленный слой — блок» при 
кинематическом способе возбуждения колебаний. На числовых примерах исследована точность и сходимость 
предлагаемой методики. 

Обсуждение и заключения. На основании выполненного математического моделирования показано, что 
разработанная методика позволяет оценить риски возникновения реальных оползневых процессов, обусловленных 
внешними нестационарными воздействиями. 

Ключевые слова: метод конечных элементов, модель Друкера-Прагера, метод Ньютона-Рафсона, динамическое 
моделирование, склон, оползневой процесс. 


Для цитирования: Гайджуров, ЦП. П. Численное моделирование поведения кинематически нестабильных склонов 
при динамических воздействиях / П. П. Гайджуров, Н. А. Савельева, Е. В. Труфанова // Адуапсеа Епетеегте 
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Введение. В 70-х годах прошлого века началось широкое использование метода конечных 
элементов (МКЭ) при анализе статической устойчивости склонов и откосов [1-3]. Как показала практика 
геотехнических расчетов МКЭ в отличие от упрощенных методик позволил учесть такие важные факторы как 
реальную геометрию и послойную структуру исследуемых объектов, а также наличие противооползневых 
сооружений и физическую нелинейность материала с заранее неизвестной областью пластичности [4]. В настоящее 
время многие исследователи стали применять МКЭ для анализа динамической устойчивости склонов реального 
заложения, а также прогнозирования рисков обвалов горных пород в подземных выработках и со склонов, 
расположенных вдоль автомобильных и железнодорожных магистралей. Также одним из актуальных направлений 
горной динамики является конечно-элементное моделирование последствий землетрясений с учетом 
кинематической нестабильности конгломераций [5]. В настоящее время существуют следующие методы 
динамических расчетов в области геотехники: 

1. ВМ (зепо@ тедисйоп тео) — метод снижения прочности [6, 7]. Разработан для расчета запаса 
прочности горного массива в физически нелинейной постановке. В ЗКМ фактические прочностные параметры 


грунта итерационно уменьшаются путем деления на некоторый коэффициент больше 1: 








С ап 
Е = ; ф’=агсапт р | 
КИ 5АМ 
где с, ф — фактические значения сцепления и угла внутреннего трения грунта соответственно; с’, ф’ — 


сцепление и угол внутреннего трения грунта соответственно после их уменьшения относительно фактических 


значений; К, — коэффициент снижения сдвиговой прочности. Величина А соответствующая предельному 


5АМ 2? 
состоянию системы, определяет нижнюю границу прочностных параметров материала. 

2. ГЕМ (Пий едаШбтат тефо4) — метод предельной кинетики, базирующийся на принципе 
Даламбера [8, 9]. Ориентирован на анализ динамической устойчивости сочлененных массивов горных пород. 

3. ТЬЕМ (фи 1Тауег @етепе тефо4 ог ЕЕМ) — метод тонкослойных конечных элементов [10], в котором 
упругопластические элементы тонкого слоя используются для моделирования поведения кинематически 
нестабильных структур. 

Анализ результатов, полученных с помощью методов 5КМ, ГЕМ и ТГЕМ, показал, что в настоящее время 
отсутствует единая концепция математического моделирования поведения структурно нестабильных 
геотехнических систем при нестационарном внешнем воздействии. Это обусловливает актуальность разработки 
методики динамического анализа систем типа «основание — ослабленный слой — блок» в конечно-элементной 
постановке с использованием нового подхода к моделированию плоскостей скольжения. 

Материалы и методы. Уравнение движения механической системы в конечно-элементной формулировке 
представим в виде [11]: 


[МО ЦИН КИ [К ЦИ} =ЗЕз+4Е@)}, (1) 


где [М], [С], [К] — матрицы масс, демпфирования и жесткости ансамбля конечных элементов соответственно; 
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{/"(Е, а, {И} — векторы-столбцы соответственно узловых ускорений, скоростей, перемещений; {АР}, 
Е (Г)} — векторы-столбцы заданных статических и динамических нагрузок соответственно в момент времени #. 


В дальнейшем полагаем, что матрицы [М] и [К] согласованные. 


Для численного интегрирования уравнения (1) используем метод Ньюмарка [12]. При этом шаг 
интегрирования по временной оси ДЕ назначаем так, чтобы с достаточной точностью учитывались вклады 


физически значимых собственных пар. В дальнейшем рассмотрим кинематические способы возбуждения 
колебаний, заданные с помощью либо модельной сейсмограммы \{И’(Р}, либо с помощью модельной 
акселелограммы {И”"(Е)}. При таком способе задания динамического воздействия второе слагаемое правой части 
уравнения (1) будет равно нулю: {А (#)1=0. 

Рассмотрим способ возбуждения механических колебаний посредством модельной сейсмограммы. 
Функцию И’ (ГР) представим в виде [13]: 

Й’ (Г) = А!е`*' зт (04), (2) 

где А — начальная амплитуда; у — коэффициент затухания; 9 — угловая частота внешнего воздействия. На 


рис. 1. приведен график функции И” (2) для значений: А = 0,01553 м; у =0,7143; @=5с (. 


И’ (0,м | 
0,006 





Рис. 1. График модельной сейсмограммы 


Результаты исследования. В качестве первого модельного примера рассмотрим задачу о вынужденных 
колебаниях склона, расположенного на основании (рис. 2). Граничные условия задачи показаны на рис. 3, где 
буквой 5 обозначена исследуемая точка. 


уз Ослабленный 





Рис. 2. Геометрия склона Рис. 3. Расчетная схема склона 


Соотношения геометрических параметров склона и основания (рис. 2) представлены в таблице 1. 
Таблица 1 
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Механические характеристики материала склона и основания: модуль деформации Е = 21 МПа; 
— : С. = -— = : 
коэффициент Пуассона у = 0,3; удельный вес у = 1702 кг/м”; сцепление с = 45 кПа; угол внутреннего трения 


ф = 15°. 
Для моделирования склона и основания используем объемные конечные элементы ЗОГ.ТО45 программного 


комплекса АМЗУ$ Месвашса|. Конечно-элементная модель для варианта с параметрами: /, = 2Н,, Н, = 10 м, 


отнесенная к глобальной декартовой системе координат, приведена на рис. 4. 

Конечно-элементная сетка построена таким образом, что на поверхности контакта смежные узлы 
основания и склона имеют одинаковые координаты, но разные номера. Это принято для того, чтобы разместить в 
этом месте ослабленный слой. Кинематическое воздействие в форме модельной сейсмограммы (2) задаем на 


каждом шаге интегрирования {, в виде узловых перемещений И’(!,) на торцовых поверхностях модели с 
параметрами: Х=Ои Х=А-+/; +[.. 
Моделирование ослабленного слоя (рис. 4) выполняем с помощью упруго-вязких комбинированных 


конечных элементов СОМВГМ14 [14]. Двухузловой элемент СОМВГМ14, состоящий из пружины жесткостью К и 


демпфера жидкостного трения с коэффициентом демпфирования с,, приведен на рис. 5. В рассматриваемом 


случае этот элемент работает только на растяжение-сжатие. 





Х 


Рис. 4. Конечно-элементная модель Рис. 5. Комбинированный конечный 


системы «основание — склон» элемент СОМВПУ14 


В каждом узле контактной поверхности (рис. 4) вдоль глобальных осей Х, У, ГД вводим элементы 
СОМВГУ14. Параметры комбинированных элементов: 
к, =30 кН/м; А, =, =9,44-10 ’ кН/м; с, = 0,5. 


В этом примере и далее вводим допущение о естественном недеформированном состоянии системы 
«основание — ослабленный слой — склон». Для вычислений используем нелинейный решатель комплекса АМЗУЗ 
Месваптса|. 

Результаты конечно-элементного моделирования в виде визуализации деформированного состояния 
системы «основание — склон» с учетом максимального горизонтального смещения и распределения амплитудных 
горизонтальных перемещений И’. (Г) показаны на рис. 6 и 7. Шаг интегрирования уравнения (1) ДЕ = 0,01 с. Как 


видно, введение ЗО) упруго-вязких элементов позволяет смоделировать эффект кинематической нестабильности 
механической системы «основание — ослабленный слой — склон» при кинематическом способе возбуждения 
колебаний. 

й’. (1), м 


-002535 
.004178 
.005822 
. 007465 
.009109 
.010752 
-012395 
-014039 
.0156582 
- 017326 





Рис. 6. Визуализация 


и Рис. 7. Распределение перемещений И’. (1) 


относительно основания 
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Амплитудное значение перемещения в точке 5 составило й = 1,/ см. Для варианта склона 


шах 


[, =Н, (рис.2) И. = 1,1 см. Графики колебаний основания и склона в исследуемой точке 5 (рис. 3) в 


тах 


направлении оси Х приведены на рис. 8. 
-- = 
И, И (0), м 


























0,02 
0,01 
0 
-0,01 
9 Е 
с 


Рис. 8. Графики колебаний в точке 5 основания И”. (Г) исклона 7, (1) 


при кинематическом возбуждении с помощью модельной сейсмограммы 


На основании приведенных графиков видно, что начиная с момента времени Е > 1,5 с наблюдается 
рассогласование колебаний основания и склона. 
Рассмотрим поведение системы «основание — ослабленный слой — склон» (рис. 3) при возбуждении 


колебаний с помощью модельной акселерограммы 7" (1). С этой целью дважды продифференцируем 
выражение (2). В результате получим: 
Й”" (Е) = Ае*' [ у’ зш(0) 1-2 7с03(01)01—зт(01)0“1- (3) 
—2 узш(01)+2с0$(01)6]. 
График функции (3) для параметров: А = 0,01553 м; у =0,7143; 90 =5с ' показан на рис. 9. 


-0,2 
Рис. 9. График модельной акселелограммы 
Кинематическое воздействие в форме модельной акселерограммы (3) по аналогии с сейсмограммой (2) 
задаем на каждом шаге интегрирования #1, в виде узловых ускорений Й”"(#.) на торцовых поверхностях модели с 
параметрами: Х=0 и Х = +1; +1, . На рис. 10 представлены графики колебаний в исследуемой точке 5 (рис. 3) 


при кинематическом воздействии в форме модельной акселерограммы. Сравнивая графики колебаний, 
приведенные на рис. 8 и 10, устанавливаем, что они практически совпадают. Это свидетельствует о корректности 
разработанной конечно-элементной модели, позволяющей описывать поведение системы «основание — 
ослабленный слой — склон» при различных способах нестационарного кинематического воздействия. 
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Рис. 10. Графики колебаний в точке 5 основания И”.. (Г) исклона 7, (1) 


при кинематическом возбуждении колебаний с помощью модельной акселелограммы 


В качестве второго модельного примера рассмотрим задачу о вынужденных колебаниях склона с 
кинематически нестабильным клиновидным включением (рис. 11). В связи с симметрией конфигурации, в 
расчетной схеме учитываем только 1/2 часть склона и включения. Граничные условия для принятой расчетной 
схемы показаны на рис. 12. Здесь буквой 5 обозначена исследуемая точка, принадлежащая одновременно 


основанию склона и клиновидному включению. 
у 







Ослабленный 
слой 
Направление 
эх скольжения 
Й’ =0 
р й=0 52 
У 
Рис. 11. Схема склона с клиновидным Рис. 12. Расчетная схема для задачи «склон — 
включением клиновидное включение» 


Конечно-элементная модель склона и клиновидного включения показана на рис. 13. Как и в предыдущем 
примере, в данном случае используем элементы 5ОТ.1245 и СОМВГ\М14 с теми же характеристиками материала. 


ь 


а) 6) в) 


Рис. 13. Конечно-элементная модель: а — склон; 6 — клиновидное включение; 





в — склон с клиновидным включением 


Функция, описывающая модельную акселелограмму, имеет вид: 


3 

— ки 

Й’ (1) =А-со$ бЕ, С 

ти сз 

где А — амплитуда ускорения; 9 — частота внешнего воздействия. На рис. 14 приведен график И’ (Г) при А = БК 
2,5 м/с’, 09 = 2Гц. Значения ускорения И (&.) на 1-том шаге интегрирования уравнения движения (1) 7 


прикладываем к узлам поверхности модели с координатой Х = 0 (рис. 12). 
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ИР. (В), м/с? 





0 0,05 0,1 0,15 0,2 


Рис. 14. График модельной акселелограммы И (1 ) 


Результат моделирования в виде распределения амплитудных значений перемещений И’. (Г) показан на 
рис. 15. Шаг интегрирования ЛЕ = 0,01 с. Графики колебаний основания склона и клиновидного включения в 


исследуемой точке 5 (рис. 12) в направлении оси Х приведены на рис. 16. 
Как видно из рис. 15 при заданном кинематическом воздействии происходит разрыв сплошности массива 
склона по ослабленному слою и клиновидное включение смещается относительно основания склона вдоль оси Х. 


Й’, (#),м 


—-02252- 
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—-017224 
—-014573 
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Рис. 15. Распределения И’. (1) в 1/2 части склона с клиновидным включением 
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Рис. 16. Графики колебаний в точке 5 основания склона И/^ (1) и клиновидного включения И’ (1) при кинематическом 
5 5 


возбуждении колебаний с помощью модельной акселелограммы 


«Дрейф» И”. (г) на рис. 16 обусловлен тем, что данная конечно-элементная модель не имеет связей, 


препятствующих смещениям вдоль оси Х. Как показано в [11], решить проблему «дрейфа» можно путем вычитания из 
значений перемещений И’. (Г) и И, (1) перемещения основания склона, которое представляет собой смещение «как 


жесткое целое». Отметим, что полученные амплитудные значения перемещений И’, (1) и И’, (1) позволяют оценить 
динамические параметры системы «склон — ослабленный слой — клиновидное включение». 


Заключение. Разработана и верифицирована конечно-элементная модель для исследования 
динамического поведения кинематически нестабильных склонов в трехмерной постановке с учетом физической 


Гайджуров П. П. и др. Численное моделирование поведения кинематически нестабильных склонов при динамических воздействиях 





нелинейности материала. 
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Введение. Рассмотрены сварные соединения в металлических крупногабаритных сооружениях (например, в 
конструкциях корпусов кораблей), подверженных малоцикловой усталости. Отмечено характерное для них 
появление мягких прослоек, которые при рабочих нагрузках значительно пластически деформируются. 
Деформация структуры металла с повреждениями, особенно в виде трещин, снижает прочность и надежность 
элементов конструкций и соединений. Предварительная деформация отрицательно влияет на пластичность, 
поэтому многое зависит от остаточной пластичности трескающегося материала. При этом с уменьшением 
остаточной пластичности снижается такой важный показатель надежности, как сопротивление материала 
распространению трещины — вязкость разрушения. 

Статья посвящена разработке модели, включающей аналитические зависимости для оценки трещиностойкости 
металлических конструкций и их сварных соединений с мягкими прослойками по пределу трещиностойкости 
при всех размерах трещин. 

Материалы и методы. Использованы теория и методы линейной механики разрушения материалов, 
структурно-механический подход. Результаты расчетов — анализировались и — сравнивались с 
экспериментальными данными и другими аналитическими решениями. Численный эксперимент выполнен для 
широко используемых в промышленности сталей феррито-перлитного класса Ст10, Ст50, 22К, СтЗсеп и др., а 
также легированных закаленных на среднюю и высокую прочность З0ХГСА, 37ХНЗА и др. 

Результаты исследования. Получены аналитические зависимости для расчета относительного предела 
трещиностойкости по трем основным известным механическим характеристикам материала мягкой прослойки 
сварного соединения. 

Обсуждение и заключение. Полученные результаты могут быть использованы для оценки трещиностойкости 
предварительно деформированных элементов конструкций и сварных соединений (в том числе с мягкими 
прослойками), работающих на поперечную нагрузку. Итоги экспериментальных данных и аналитических 
расчетов показаны в безразмерном виде, что позволяет получать инвариантные результаты по отношению к 


пределу трещиностойкости. 


Ключевые слова: сварное соединение, мягкая прослойка, механическая неоднородность, деформация, предел 


трещиностойкости, длина трещины, контактное упрочнение, относительное удлинение. 
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Введение. В сварных соединениях, подверженных циклическим нагрузкам, часто возникают и 
развиваются трещины, которые могут стать причиной опасного состояния конструкций и их разрушения. В [1] 
приводятся примеры формирования таких трещин в корпусах судов. Кинетика образования в сварных швах 
макротрещин и их развитие зависят от множества факторов, что требует проведения комплексных 
исследований трещиностойкости сварных соединений и разработки соответствующих математических моделей. 

На сопротивление металла распространению трещин существенно влияют появляющиеся после сварки 
так называемые мягкие прослойки, подверженные пластическому деформированию. Для разных участков 
сварного соединения характерна большая или меньшая степень неоднородности механических характеристик 
металла [2]. Она зависит от многих начальных условий технологического процесса. В их числе — скорость 
охлаждения зон сварного соединения и обработки свариваемой стали, а также содержание легирующих 
элементов, углерода и др. О. А. Бакши рассматривал некоторые проблемные вопросы вязкого и хрупкого 
разрушения механически неоднородных сварных соединений, влияния дефектов и микротрещин на процессы, 
происходящие в структуре металла при нагружении, и пр. Исследователи решили ряд частных задач по 
увеличению прочности и долговечности сварных соединений с прослойками [3—6]. 

При установившемся режиме сварки ширина зон термического влияния (ЗТВ) и механические свойства 
металла незначительно меняются по длине сварного соединения. Поэтому целесообразно рассматривать 
неоднородность характеристик в разных зонах в сварном соединении при поперечной нагрузке. Для 
закаливающихся или прошедших термическую обработку (термоупрочненных) металлов характерно появление 
закалочных структур или разупрочненных зон соответственно. При этом в зависимости от свойств стали в 
исходном состоянии и параметров технологического процесса могут образовываться мягкие зоны. Здесь предел 
текучести от и предел прочности о» металла меньше, чем в соседних прилегающих (так называемых твердых) 
областях от и оз. Таким образом, мягкая зона будет вовлекаться в пластическое деформирование между 
областями с характеристиками о. > от, 0, > 0» раньше соседних, а степень ее влияния на прочность 
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соединения будет зависеть от схемы приложенной нагрузки‘. Деформирование мягкой прослойки снижает 
сопротивляемость распространению трещины — инициируется интенсификация хрупкого разрушения, при 
котором критическая длина трещины уменьшится в разы [7]. 

Цель данного исследования — разработка аналитических зависимостей и математической модели для 
оценки трещиностойкости сварных соединений с мягкой прослойкой в результате ее предварительного 
пластического деформирования. 

Допустим, что некоторое ограничение линейных размеров сварного образца (кроме толщины) под 
поперечной нагрузкой заметно не влияет на рассматриваемые в структуре металла процессы и явления (кроме 
образования мягкой прослойки), вызванные термическим циклом сварки. 

Материалы и методы 

Расчетные предпосылки. Рассмотрим изменение физико-механических характеристик металла мягкой 
прослойки при растяжении стыкового соединения поперек шва достаточно большой протяженности (рис. 1). 
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Рис. 1. Схема зависимости характеристик о, ‚ би Ц от относительной толщины прослойки 


В упругой стадии нагружения прослойка и соседние участки рассматриваемой области деформируются 
однородно. При достижении уровня напряжений, соответствующих пределу текучести от, в ней образуется 
пластическая деформация, а соседние участки остаются в упругом состоянии. Для упрощения примем угол 
наклона грани прослойки В = 0 (рис. 1). При дальнейшем повышении нагрузки и деформации растет 
коэффициент поперечной деформации и (коэффициент Пуассона) мягкой прослойки. В итоге он будет больше, 
чем у соседнего участка металла. По мере развития пластической деформации в прослойке и -> 0,5, ав 
областях с упругим состоянием металла и = 0,3. Из-за неодинаковой поперечной деформации в сварном шве 
возникают касательные напряжения. Максимальные их значения концентрируются в районе плоскостей раздела 
мягкой прослойки и соседних участков. Касательные напряжения препятствуют поперечному сужению мягкой 
прослойки в направлении толщины листа. Чем уже прослойка (т. е. чем меньше отношение ге = № /5, рис. 1), 
тем меньшее поперечное сужение она получает к моменту возникновения истинных разрушающих 
напряжений о. Так возникает эффект контактного упрочнения, который обеспечивает повышение несущей 
способности. 

Увеличение разрушающей силы имеет пределы. Соседние с мягкой прослойкой более прочные участки 
при определенных условиях тоже будут пластически деформироваться. Чем большей прочностью обладают 
соседние зоны и чем меньше относительная толщина прослойки 2, тем заметнее эффект контактного 
упрочнения. Его усиливают более твердые прослойки, находящиеся вблизи относительно мягких. 

Относительное поперечное сужение ф в мягкой прослойке и абсолютное удлинение образца 6 зависят 


от г и свойств металлов в соединении. В широких прослойках, когда еще нет контактного упрочнения, ф 
остается постоянным при уменьшении г (рис. 1). А 6 постепенно падает по мере уменьшения г в общей 
ширине сварного образца. В области контактного упрочнения 4 резко уменьшается, так как растут поперечные 
касательные напряжения, препятствующие сужению прослоек в направлении толщины материала. 
Относительное удлинение образца 6 при этом также сначала уменьшается. Когда же реализуется существенная 
жесткость напряженного состояния в мягкой прослойке, д возрастает, поскольку в пластическую деформацию 
вовлекается металл твердой прослойки. Отметим, что при сварке термоупрочненных сталей дбь мягкой 


прослойки оказывается выше, чем 6-. Но ее пластическая деформация значительно больше, и остаточная 
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пластичность после нагружения до напряжений о, > о > от может стать критической. В результате 
значительно уменьшается длина критической трещины, которая может возникнуть под действием дальнейших 
рабочих циклических напряжений на стыке мягкой и твердой прослоек. При этом существенно возрастает 
опасность хрупкого разрушения и следует оценить его возможность (предел трещиностойкости”” [8, 9]) в 
области с дефектом в виде трещины. При этом учитываются условия плоской деформации, контактного 
упрочения и пластически деформированной мягкой прослойки. По результатам исследования деформирования 
слоев композитного материала методом конечноэлементного моделирования [3] можно заметить, что 
наибольшая деформация мягкой прослойки образуется около стыка ее с твердой прослойкой. Особенно четко 
это проявляется при средних значениях относительной толщины прослойки ге. 

Иногда при статической нагрузке удается найти оптимальное значение ге и достичь равнопрочности 
соединения. Однако имеющиеся в мягкой прослойке дефекты или несплошности после некоторой ее 
деформации могут развиться в трещины с критической длиной при условии, что оставшаяся часть 
докритической деформации относительно мала. Степень контактного упрочнения в случае плоской 


деформации определяется по формуле’ [2]: 
К. =. (1) 


Временное сопротивление сварного соединения с мягкой прослойкой: 

05 = 6, 'К», (2) 
где ов — временное сопротивление мягкой прослойки. 

Естественно, что для обеспечения равнопрочности при статической нагрузке требуется как минимум 
выполнение условия 0% = оъ. В этом случае предельное (минимальное) значение коэффициента контактного 
упрочнения определяется отношением временных сопротивлений: 

о 0. (3) 


Используя (1) и (2), найдем предельные значения г при которых достигается равнопрочность 


пред> 
соединения в условиях плоской деформации: 
1 


пред — ЗК, - п) (4) 


Однако необходимо иметь в виду, что повышение прочности сварного соединения с мягкой прослойкой 


ге 


ограничивается истинными напряжениями отрыва для мягкого металла. 

При сварке упрочненных сталей ширина разупрочненной зоны зависит от погонной энергии и 
регулируется выбором подходящего режима сварки. Для каждой толщины металла, способа и режима сварки 
есть определенная ширина разупрочненного участка, при которой обеспечиваются максимально возможное 
контактное упрочнение, а также равная прочность сварного соединения и основного металла. 

Исследование участка разупрочнения при аргонодуговой и электронно-лучевой сварке металла 
толщиной 4,5—8 мм показывает: 

— твердость на нем определяется только температурой нагрева; 

—щ прочность соединений зависит не от уровня твердости разупрочненного участка, а от его 
ширины Й (рис. 1). 

При этом следует учитывать, что участок разупрочнения плавно переходит в более прочные участки 
зоны термического влияния. 

Изучение кинетики и механизма разупрочнения в ЗТВ позволили установить, как на прочность 


сварных соединений влияют свойства участков, окружающих мягкую прослойку. Экспериментально обоснован 


СС 


новый критерий, определяющий временный предел прочности сварного соединения оър, 


— это относительная 


ширина мягкой прослойки ВБ: /Бу. Предложенная эмпирическая зависимость имеет вид: 


Бо Бо" 
0<° = ом + 0,Зом (>) + 0,30“ (>) (5) 
Ь Ь. 
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где Ор — предел прочности сварного соединения, МПа; о» — предел прочности мягкой прослойки, МПа; Бу — 
ширина прослойки при равной прочности сварного соединения и основного металла, м; Б, — текущая ширина 
прослойки, м. 

Из последнего уравнения следует, что при ВБ, = ББ. мягкая прослойка упрочняется максимум на 60 %, 
т.е. о, = 1,60ь, так как дальнейшее упрочнение не имеет смысла. 

Результаты экспериментов указывают на снижение выносливости сварных соединений с мягкой 
прослойкой, образованной после сварки термоупрочненных сталей У4 и Уб. Это говорит о существенном 
снижении несущей способности [10]. У этих сталей пределы выносливости сварных образцов без концентрации 
напряжений практически совпадают с пределами выносливости таких же образцов из неупрочненной стали. То 
есть, несмотря на термическое упрочнение сталей (способствующее повышению прочности), вибрационная 
прочность сварных соединений на базе 10° циклов оказывается практически одинаковой. Следовательно, 
можно предположить, что пластическая деформация в мягких прослойках вызывает раннее образование 
макротрещины, критический размер которой [+ резко ограничен происшедшими пластическими деформациями. 

Трещиностойкость пластически деформированного материала мягкой прослойки до значений о < бь 
может оцениваться двухпараметрическим критерием разрушения — пределом трещиностойкости, который 


.. `. 4 
рассчитывается по известной аналитической зависимости [2, 11]: 


бс\- 
к = |1- (52), (6) 
Ов 
где о; — критическое напряжение разрушения в брутто-сечении; Кс — критический коэффициент 


интенсивности напряжений (постоянная материала мягкой прослойки). 

После пластической деформации мягкой прослойки о. будет зависеть от характеристики дефекта. 
Примем, что эта характеристика — критическая длина [‹ зародившейся трещины. Тогда ос = К. / те. 
Однако, чтобы воспользоваться (6), необходимо знать степень деформации прослойки, так как от нее в первую 
очередь будет зависеть Кс. Величину [‹ будем считать неизвестной. Ее нужно определить в зависимости от 
деформации мягкой прослойки. 

Авторы работ” [11], исследуя структурный параметр разрушения и определяя его для дисковидной 


трещины, приходят к выражению для оценки предела трещиностойкости в одномерном виде: 


Ке _ У21И- т) (1) 
К, — агссоз()' 

гле | =а/(а-+ а) — безразмерный параметр, заданный диапазоном изменения (0 <1п< 1); а — радиус 
дисковидной трещины. 

Из последнего равенства следует, что для больших трещин характерно относительно малое значение 4, 
а. Разрушение происходит с малыми пластическими деформациями в нетто-сечении, а К; = Кс. Обратная 
ситуация складывается при относительно небольших трещинах, если для разрушения понадобится большая 
пластическая деформация в нетто-сечении, ит = 0, Кс «< К\с. 


Представим теперь, что мягкая прослойка подвергается циклической нагрузке. В этом случае 
используем концепцию эффективного порогового коэффициента интенсивности напряжений Ков для 
вычисления предельной амплитуды цикла нагружения о. при уровне средних напряжений 0, > 0щ + 0; > 


от. С этой целью рассчитаем длину начальной макротрещины Ё, соответствующей предельной амплитуде 


напряжений в пластически деформированном материале до ож: 





А о 
п \07, 07, (8) 


= е —= 
при К+но = Кио," = —1, 
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где о;, — предельная амплитуда напряжений; о,;‚ — циклический предел текучести; 7 — коэффициент 


асимметрии цикла; и — коэффициент Пуассона. 
Значение длины начальной макротрещины [ следует подставить в расчетную зависимость для 


вычисления предела выносливости [12]: 
т. 0,5 
О И 
бр = ор | ПЁ (”) + 0,81. (9) 


Курт 

Решением будут искомые значения предельных амплитуд. Характеристики, входящие в (9), 
вычисляются согласно [12] или’. Далее не составляет труда определить вклад предварительной пластической 
деформации мягкой прослойки в эквивалентную поврежденность с позиции оценки длины макротрещины. 

Результаты исследования 

Построение математической модели. Опишем предлагаемую последовательность разработки 
аналитических зависимостей. Оставшиеся пластические свойства элемента конструкции или ее части можно 
оценить по деформации, используя относительное удлинение для области материала. Будем считать, что на 
малом участке мягкой прослойки имеет место равномерная пластическая деформация. Она связана с 
критическим относительным удлинением до разрушения 65(10) известной зависимостью [13]: 

е = ш(1 +6..). (10) 

Положим, что материал мягкой прослойки пластически деформируется. Рассмотрим случай, когда 
пластическое деформирование происходит только в первых циклах нагружения. То есть область мягкой 
прослойки вытягивается пластически под внешней поперечной номинальной нагрузкой, а далее она работает 
под упругими циклическими напряжениями. Так развивается усталостная трещина. В мягкой прослойке ее 
может спровоцировать высокая пластическая деформация и (или) дефект. Из-за снижения предельной 
деформации материала мягкой прослойки предварительное пластическое деформирование приводит к 
уменьшению критического раскрытия в вершине трещины. По результатам работы [14] можно заметить, что 
очень близка к линейной зависимость критического раскрытия вершины трещины 6. (КРТ) при деформации 
стали вытягиванием. Этот факт наблюдается как для алюминиевого сплава, так и для аустенитной стали. Итоги 
упомянутых экспериментов позволяют предположить, что при предельной равномерной пластической 
деформации этих материалов д. = 0. 

От степени предварительного деформирования существенно зависит предел выносливости. Такую 
зависимость сложно получить теоретически [15], потому что до области средней предварительной 
пластической деформации материала предел выносливости достигает минимума, а после возрастает и может 
достигать показателей как у недеформированного материала. Таким образом, сам факт появления трещины, 
видимо, играет основную роль в усталости. Нераспространяемые трещины даже при большой предварительной 
пластической деформации только в редких случаях ведут к разрушению. Это очевидно из примеров 
повреждений судовых конструкций, приведенных в [|]. 


Учтем, что 6. =0 в пределе пластической деформации, и используем линейную аппроксимацию 
зависимости 61 (=): 








61с 
6(=) = 6 — Е. (11) 
тах 
В вершине трещины КРТ определяется известной зависимостью линейной механики разрушения: 
Кс 
61е ЕЛ | 
АЕ (12) 
Коэффициент Л [12] находится по формуле: 
— 2 + 
а А-В ег | (13) 
1,24по. \Атсе) 
Здесь Е — модуль упругости материала; Кисе — напряжение микроскола деформированного 


материала; Р — коэффициент, учитывающий повышение первой главной компоненты напряжения для случая 
сложнонапряженного состояния; 4 — коэффициент, показывающий отношение второй главной компоненты 
напряжений к первой; тт — коэффициент упрочнения. 

Для указанных характеристик и коэффициентов существуют вполне определенные расчетные 
зависимости по основным механическим характеристикам материала. Критический коэффициент 
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интенсивности напряжений в (12) рассчитывается по формуле: 


(14) 





где 4, — средний диаметр зерна исходного материала мягкой прослойки. Подставим последнее выражение и 
(13) в (12). Сделаем необходимые сокращения. Учтем, что в (12) и (11) можно пренебречь малой составляющей 
упругих деформаций. В итоге получим: 


ЕО} 
51. (&) = 5(1 — 22а, (1 Е тов) (15) 


Деформация = вычисляется по зависимости (10). Найдем зависимость К\с(=). Из (12) имеем Ку с (=) = 
\/Ё0.61с(=)/Л. Подставим известные выражения и получим: 


1 
Ютсе О т" 2Ов 16 
те) 6.18па,0? (1— (16) 
от ы 20? ( т 












Кас (Е) = 


ИЛИ 








1 
0,618па,2 (Ютее Вт" 
Кис (5) = | "ее .( 


от (1-24) 


После упрощения можно получить простую зависимость критического коэффициента интенсивности 








) (17) 


напряжений К\ с от текущей величины 6 вытягивания мягкой прослойки: 








Кс (8) = К, = Ки |1 (8) 
где 6 — удлинение (вытягивание) материала мягкой прослойки, доли. 
В формуле (18) К. (6) — предел трещиностойкости К.. Его относительное значение Кс (6)/Клс 


показывает, как снижение вязкости разрушения материала зависит от удлинения (вытягивания) элементарной 
области мягкой прослойки. Для простоты расчетов можно использовать усредненное удлинение мягкой 
прослойки под полезной нагрузкой. Отметим, что локальность пластического удлинения в сварном соединении 
мягкой прослойки на ДА относительно Й может существенно влиять на ее усталостную прочность. Для 
упрочненного основного материала это влияние, наоборот, будет незначительным (рис. 1). 

При 6 = 6: отношение К/К: должно давать ноль, поэтому приравняем (18) к нулю и запишем 
относительно 6 расчетную формулу: 

б- = ехр(100/в,) — 1. (19) 
Здесь о, подставляется в МПа. 

Анализ результатов и проверка адекватности полученных зависимостей. Можно признать 
удовлетворительными результаты проверки согласования расчетных характеристик по полученным 
аналитическим зависимостям со справочными данными конструкционных сталей феррито-перлитного класса в 
состоянии поставки. Речь идет о тонколистовом прокате из стали Ст10, Ст50, 22К, СтЗсп, 37ХНЗА, ЗОХГСА 


и др. (таблица 1). 


Таблица 1 
Расчетные и справочные характеристики сталей 


Го |8 [09 | 15 | 356 [ >И 
ЕО 


Гав Г 39 [09 [9 [203 [2% 
ЕСО ОЕ ОС ООО ВО ОСЛО 
ВтхнаА [ча | 9 [ом | 7 [06 | 8 [3 [ № о 
ВохгсА | 175о | 1аво [009 | 3 [0 [5 [6 [7 
ВохгсА Го Г пов Гоа [4 6 [5 


"Расчет по (19); “расчет по т = (0,75 - 190, (4 + 1.4ф-)/во2]}/19 [105 т (——) /(200 + 0,59%2)| [15]: 


***расчет по (14). 
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Выполненный ранее анализ [12] подтверждает, что коэффициент А, рассчитанный по (13) при КРТ в 
вершине 0,6184», согласуется с экспериментальными данными для сталей 22К и Ст10. В первом случае ^, равен 
0,22, во втором — 0,23. 

Рассмотрим предварительную деформацию тонких листов при плоском напряженном состоянии с 
условием, что зерно получит относительное удлинение е, в продольном направлении. В этом случае средний 


его диаметр в том же направлении будет больше в (1 +е,) раз, а в поперечном уменьшится в (1 = е„). То 
есть, если в продольном направлении средний диаметр возрастет в 2,618 раза, то в поперечном он уменьшится 
в 1,618 раза. Здесь поперечное направление совпадает с направлением движения трещины. Легко видеть 
линейность зависимости, по которой определяется увеличение среднего диаметра зерна. Из этого следует, что 
при равномерном одноосном пластическом деформировании критическое раскрытие в вершине трещины тоже 
должно снижаться линейно, в зависимости от относительного удлинения е,. При прокатке стальных листов 
такой выраженной линейной зависимости не обнаруживается. Очевидно, это связано с поперечными 
деформациями, для которых характерна приплюснутость зерен. 
Диаграммы изменения К. ‹ как функции от 6 для сталей (таблица 1) приведены на рис. 2. 









К,„МПА.-м“” 120 — 22к 
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20 | 
0 | 
0 5 10 15 20 25 30 35 40 
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Рис. 2. Изменение трещиностойкости в зависимости от предварительного удлинения сталей 


Критический коэффициент интенсивности вычислен по (14) с использованием данных таблицы 1. Это 
позволяет определить уровень К1с при пластической деформации вытяжкой и сопоставить его со значениями 
других сталей. Так, сравнение стали ЗОХГСА, обработанной на величину предела текучести о. = 1360 МПа, и 
стали 10 показывает, что при использовании низкопрочной стали 10 уровень Кс можно снизить только при 
больших деформациях удлинения (около 25 %). На практике в конструкциях так бывает достаточно редко. 
Однако при целенаправленной вытяжке или холодной прокатке возможно существенное снижение К1с, и это 
следует учитывать. 

Если стали ранжировать по предельной пластичности и сопоставить для них максимальные значения 
К1с, то окажется, что пластичность убывает от одной марки стали к другой гораздо быстрее, чем меняется 
критический коэффициент интенсивности напряжений. Так, К1с для сталей ЗОХГСА и 10 в исходном состоянии 
поставки отличается в -2 раза, а относительное удлинение 6 — в -5 раз. Появление пластичности в вершинах 
различных концентраторов напряжений неизбежно, поэтому для обеспечения необходимой надежности 
элементов из стали ЗОХГСА целесообразно контролировать пластические деформации в конструкции в -5 раз 
точнее. Кроме того, необходимо учитывать, что снижение вязкости в -2 раза приведет к снижению критической 
длины трещины в -4 раза. Причина — низкая вязкость и пластичность стали ЗОХГСА, для эксплуатации 
которой придется в -5раз увеличить коэффициент запаса прочности. Последнее не означает, что 
высокопрочные стали как конструкционный материал хуже сталей средней прочности. Проблема заключается в 
выборе оптимального сочетания вязкости, пластичности и прочности. Полученные в настоящем анализе 
расчетные зависимости дают возможность оценить данные характеристики расчетным путем. Это важно для 
рационального выбора материала на этапе конструирования, определения оптимальных режимов сварки и 
термообработки. 

Сопоставим результаты расчетов по (18) с известной зависимостью (7)“ (рис. 3). Первая функция 
зависит от оставшейся пластичности стали, а вторая — от относительной длины трещины в концепции 
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структурного элемента. Чтобы удобно представить результаты для их сравнения, введем безразмерный 












параметр: 
1 =1-а/(а+а). (20) 
К/Ки | 7 [ о Ь=50 ММ 
: + 65-100мм 
х х 50-200 мм 
1,0 ых 0’ 
= т 

0,8 | 
0,6 т 
0,4 Е 
0,2 Е 

0 | | | | | | 1 

0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 
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Рис. 3. Зависимости относительного предела трещиностойкости от безразмерного параметра дисковидной трещины 1” и от 
безразмерного параметра предварительного равномерного удлинения материала 6*. Линии — аналитические зависимости 
(7), (18). Маркеры — экспериментальные данные для компактного образца? [11, 16] с размерами В 


После подстановки (20) в (7) получаем инвариантное выражение для среднего значения. Если известно 
|”, то можно найти отношение а/(а + а) по (20). Результаты сравнения показывают, что функции практически 
совпадают. В относительных координатах для них характерна инвариантность с хорошим приближением. Обе 
функции хорошо аппроксимируют экспериментальные данные предела трещиностойкости, полученного на 
образцах с характерными размерами В [13]. 

Рассмотрим условия плоской деформации при известной критической длине трещины [х‹. Для феррито- 
перлитных сталей с о. < 400 МПа ее приближенно можно оценить как [, = 2К?./п/о2 и легко найти 


критическое напряжение внешней нагрузки ос = Ке/\/пЁс: 


(21) 





Уравнение отражает закон уменьшения критических напряжений ос ниже значения предела 
выносливости материала с длиной трещины [Ё: в зависимости от величины предварительного удлинения 
материала. В этой же постановке снижение критической длины трещины в предварительно деформированном 
материале мягкой прослойки можно оценить по зависимости: 





= (==) (1- ша +5)}, (22) 
п‘ д 100 ( ) 
где [ — половина длины для сквозной тещины; с — критическое напряжение для катастрофического 


распространения трещины длиной [. 

При о =о_, из (22) следует, что длина Ё катастрофического распространения трещины прямо 
Ов 
100 
Полученные зависимости в системной совокупности представляют предлагаемую математическую 





пропорциональна значению (1 —— Ш + 5)) и константе материала мягкой прослойки [с. 


модель оценки трещиностойкости сварных соединений с мягкими прослойками. 

Обсуждение и заключение. При сварке термоупрочненных сталей для корпусных и других листовых 
конструкций образуются мягкие прослойки, металл которых пластически деформируется в первых циклах 
нагружения. Методами численного анализа с применением САЕ-средств'” можно определить деформации или 
увеличение ширины мягкой прослойки сварного соединения под действием рабочих нагрузок. Надежность 
работы соединения с мягкой прослойкой зависит от степени пластической деформации, ее упрочнения, 
жесткости напряженного состояния с учетом приграничных зон твердых прослоек, а при появлении трещины 


7 Матохин Г. В., Молоков К. А. Указ. соч. 
0 САЕ (от англ. сотрщег-а14е4 епошеегие — компьютерная инженерия). Так называют программы и программные системы для решения 
инженерных задач. 
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— от ее критической длины Ё. При вытягивании мягкой прослойки [ снижается прямо пропорционально 
константе материала мягкой прослойки [; и выражению, зависящему от степени вытягивания 6, т. е. от 
относительного предварительного удлинения. 

Получена простая аналитическая зависимость критического относительного удлинения до 
разрушения (19). Она обеспечивает результаты расчетов, которые удовлетворительно согласуются с 
экспериментальными данными для большинства конструкционных сталей (по состоянию поставки). Формула 
позволяет оценить справочную характеристику 65 образца по единственному известному параметру — 
временному пределу прочности. Для данной зависимости характерна теоретическая закономерность, которая 
хорошо согласуется с параметрами конструкционных сталей. Однако в некоторых случаях (специальные 


режимы термообработки или отпуска) она может дать относительную погрешность до - 30-40 %. 

Предложена аналитическая зависимость изменения предела трещиностойкости мягкой прослойки от ее 
предварительного пластического деформирования. Сопоставлены: 

— расчетные значения предела трещиностойкости; 

— экспериментальные данные, полученные на компактных образцах для испытаний; 

— более ранние аналитические решения. 

Результаты сравнения убедительно доказывают применимость предложенного аналитического 
решения. Рассмотрено действие закона снижения предела трещиностойкости при относительном 
предварительном пластическом удлинении вплоть до разрушения (65) и изменении параметра (1”) 
предварительно выращенной трещины в материале вплоть до [.. Для данного случая установлена 
инвариантность в относительных координатах при условии плоской деформации и отсутствия влияния 
ограниченности размеров образца или конструкции. Учет ограниченности размеров образцов и разных схем 
нагружения существенно меняет упомянутую зависимость ' [13, 17]. 

Предложенный метод оценки предела трещиностойкости сварных элементов с мягкими прослойками 
может быть использован в инженерных машиностроительных расчетах по минимизации металлоемкости, а 
также для корабельных конструкций, подверженных малоцикловой усталости. 
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Введение. Снабжение людей качественной питьевой водой всегда имело первостепенное значение. Но её 
транспортировка по трубопроводным системам зачастую бывает связана с проблемами, вызываемыми, 
например, температурой воды и окружающей среды, а также возможностью возникновения гидравлических 
ударов на отдельных участках труб. В болышей степени эти проблемы касаются систем, в которых 
используются трубы из полиэтилена. Температура — это один из наиболее важных факторов, влияющих на 
гибкость полиэтиленовых труб, и этот фактор оказывает влияние не только на проектирование водопроводных 
сетей, но и на капиталовложения в их развитие. Целью настоящих исследований явилось изучение влияния 
температуры воды и окружающей среды на свойства материала труб и скорость распространения 
гидравлической ударной волны в полиэтиленовых трубах. 

Материалы и методы. В проведенных авторами опытах применялся метод натурных исследований, когда для 
испытания используются специализированное оборудование и специально изготовленные для этих целей 
образцы. В данном случае использовались образцы труб из полиэтилена высокой плотности, которые 
подвергались испытаниям на растяжение на разрывной машине, причём каждый опыт проводился троекратно. 
В процессе проведения опытов образцы подвергались воздействию определённых температурных режимов (как 
внешних, так и внутренних), при этом исследовалось также влияние гидродинамического давления жидкости в 
трубе в результате изменения скорости жидкости в её секциях. Для этого в образцы подавалась жидкость под 
определённым давлением, чтобы выяснить влияние на трубы эффекта, известного как гидроудар. 


Результаты исследования. В ходе исследования было установлено, что значение модуля упругости 
полиэтилена высокой плотности РЕ100 уменьшается с повышением температуры воды, причём величина 
снижения при температуре 60 °С достигает 60,21 % по сравнению с его значением при температуре воды +4 °С. 
В результате опытов по определению изменения модуля упругости полиэтилена с повышением температуры 


было выведено экспоненциальное уравнение для расчета значения полиэтиленового коэффициента как 


—0,011 


функции времени Ё =1,312е с коэффициентом корреляции А” = 0,988, а для расчета значения скорости 


распространения гидравлической ударной волны было выведено экспоненциальное уравнение как функции 


2" с коэффициентом корреляции А” = 0,987. 


времени С = 275,9е 
Обсуждение и заключения. В результате экспериментов было установлено, что вредное воздействие на стенки 
труб оказывает такое явление, как гидроудар, возникновение которого по возможности необходимо избегать 
ещё на этапе проектирования водопроводной сети. В ходе исследования было определено, что с увеличением 
температуры значения модуля упругости полиэтилена уменьшались при одновременном понижении значений 


скорости распространения гидравлической ударной волны. 
Ключевые слова: полиэтиленовые трубы, температурный диапазон, гидроудар. 
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Введение. Температура протекающей жидкости и температура окружающей среды оказывают 
существенное влияние на свойства материала труб, изготовленных из полиэтилена высокой плотности. При 
повышении температуры материал трубы становится более гибким, что приводит к снижению скорости 
распространения гидравлической ударной волны, а значит, и к снижению давления. Промышленно 
изготовленные полиэтиленовые трубы из РЕ100 (НОРЕ) являются термопластами, так как их пластичность 
варьируется в соответствующих пределах, зависящих от рабочих условий. Эти пластмассы имеют ряд 
преимуществ, в том числе сравнительно небольшой удельный вес, они дают возможность работать при низких 
температурах, способны выдерживать умеренное давление, обладают высокой коррозионной устойчивостью 
(особенно к кислотам и щелочам), хорошими электроизоляционными свойствами, легко окрашиваются и 
просты при установке. Они изготавливаются сразу большой длины, что сокращает при их прокладке 


количество соединительных элементов, которые считаются дорогими и хрупкими и вызывают проблемы, 
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связанные с утечками [1—6]. Полиэтиленовые трубы относятся к гидравлически гладким трубам, в которых 
потери энергии на трение в процессе протекания жидкости намного меньше, чем в металлических трубах, что 
приводит к экономии энергии перекачивания. Эти трубы характеризуются очень малым диапазоном 
шероховатости (К=0,001-—0,008 мм), что в 3—20 раз меньше шероховатости новых стальных труб. Температура 
является одним из наиболее важных факторов, влияющих на свойства гибкости полиэтиленовых труб, и этот 
фактор не только влияет на проектирование и капиталовложения в различные сети снабжения питьевой водой 
при стабильности потока, но и вызывает нестабильное его течение в виде гидроудара в полиэтиленовых 
трубах [7-10]. Международные стандарты определяют свойства полиэтиленовых труб в сетях питьевого 
водоснабжения при температуре 20-23 °С, в то время как европейские стандарты определяют их при 
температуре 10 °С. Модуль упругости полиэтиленовых труб относительно велик при низких температурах 
водного потока, значение его понижается, когда при гидравлическом ударе температура потока повышается. 
Поэтому гидравлический удар может оказывать значительное влияние на стенки полиэтиленовых труб в сетях 
питьевого водоснабжения [11—14]. Цель настоящих исследований — изучение влияния температуры жидкости 
и температуры воздуха, окружающего полиэтиленовые трубы, на плотность, другие свойства материала труб и 
на скорость распространения в них гидравлической ударной волны. 

Методы и материалы исследования. Гидравлический удар. Гидравлическим ударом называют 
увеличение или уменьшение гидродинамического давления жидкости в трубе в результате изменения скорости 
жидкости в еб секции. Этот удар может создавать большое давление, которое необходимо учитывать при 
расчете толщины стенок трубы. Сжимаемость жидкости и склонность стенок трубы к деформации придают 
гидравлическому удару упругость, поскольку он представлен в виде волн давления, распространяющихся по 
трубе в нестабильном потоке. Различают два типа гидроудара: 


1. Прямой гидравлический удар опасен и возникает, когда время остановки меньше периода 
гидравлической волны, {. <Ц. 
2. Косвенный гидравлический удар, который не опасен и срабатывает, когда время отключения больше 


периода гидравлической волны, 1 >Ц: 
ЕС. (1) 
где № — период гидравлической волны, с; [Г — длина трубы, м; С ‚ —_ скорость распространения 


гидравлической ударной волны, м/с; #. — время закрытия или открытия клапана, с. 


Прямой гидравлический удар возникает, когда на конце трубы внезапно открывается или закрывается 
клапан, соединенный с широкой поверхностью резервуара, или внезапно прекращается подача электроэнергии 
на насосную станцию. Изменение давления в результате прямого гидравлического удара определяется 


соотношением Жуковского-Леви: 
АР „ =р, С, И, (2) 


шах 


где ДР 


шах 


2 
—Щ величина изменения высокого или низкого давления, Н/м 5 


5 
р, — плотность жидкости при установившемся расходе, кг/м; С, — скорость распространения 


гидравлической ударной волны, м/с; /, — установившаяся скорость потока жидкости в трубе, м/с. 


Что касается непрямого гидравлического удара, то он возникает при медленном открытии или 
закрытии клапана на конце трубы, соединенной с широкой поверхностью резервуара, или при увеличении 
времени остановки работы насосной группы на насосной станции. Величина изменения давления в результате 


непрямого гидравлического удара определяется следующим соотношением: 
АР. „ =2р, -Ё.-И, / ЕТ, (3) 


Величина скорости распространения гидравлической ударной волны в полиэтиленовых трубах 


определяется по следующей формуле: 


(4) 
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где М — постоянная консолидации при О/е>40 

Для поддержки одной стороны трубы М =1,25-и, М =1-и/2, всей трубы — М =1- м”; для труб без 
опоры (с компенсаторами) — М=1. Здесь и — модуль Пуассона, составляет 0,45 для труб из НОРЕ [3—4]; Ё/, — 
коэффициент упругости воды, Н/м”; р ; — ПЛОТНОСТЬ ВОДЫ, кг/м; 2 — внутренний диаметр трубы, м. 

Скорость распространения гидравлической ударной волны в полиэтиленовых трубах находится в 
диапазоне 180-370 м/с [1, 9]. В таблице 1 представлены данные о физических свойствах воды при атмосферном 
давлении на еб поверхности. На рис. | показано уменьшение плотности воды с увеличением температуры, а на 
рис. 2 — увеличение значения модуля упругости воды с повышением её температуры. 


Таблица 1 


Физические свойства воды при атмосферном давлении на её поверхности 


томперурьес Гоа 
оз 109 
2.02 





Модуль упругости воды, 
ГН/м? 


1002 
1000 
998 
996 
994 
992 
990 
988 


Плотность воды, кг/м” 


986 
984 


982 
0 10 20 30 40 50 60 70 
Температура, °С 


Рис. 1. Изменение плотности воды с повышением температуры 
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322 Рис. 2. Зависимость значений модуля упругости воды от температуры 
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В таблице 2 приведен модуль упругости материала трубы из полиэтилена при 16 °С [2, 10]. 
Таблица 2 
Модуль упругости материала трубы из полиэтилена при температуре 16 °С 


Модуль упругости Модуль упругости Модуль упругости 
полиэтилена, № полиэтилена, МН/м* полиэтилена, №/Ё2 





150 000 1 030 22 000 000 


Лабораторные испытания. Из полиэтиленовой трубы высокой плотности НОРЕ был изготовлен 


21 образец. В ходе проведенных испытаний были измерены значения определенных характеристик, которые 


приведены в таблице 3. 


Таблица 3 
Технические характеристики труб из НОРЕ 


Вю | ваши | бы ок | БМР 
го бе 


Для проведения исследований были выполнены следующие действия: 

— выбраны и приняты значения температур, при которых производились испытания: 4, 10, 20, 30, 40, 
50, 60 °С: 

— для каждой из выбранных температур испытывалось по три образца; 

— в ходе испытаний определялись физико-механические характеристики, средние значения которых 
приведены в таблице 4; 

— образцы на растяжение нагревали в течение часа в цифровой электрической печи, затем растягивали 
с помощью специального устройства, полученные результаты представлены в таблице 4 и на рис. 3; 

— полиэтиленовую трубку диаметром 110 мм заполняли водой при определенных температурах до 
момента её разрыва и рассчитывали относительную деформацию по формуле = =АР/Д (таблица 4). 

Изменения упругих растягивающих напряжений в образцах труб из полиэтилена с внешним диаметром 
110 мм при номинальном давлении 10 Бар и их деформации, измеренные при выбранных ранее температурах, 
анализировались в течение года. Полученные данные представлены в таблице 5 и на рис. 4. 


Таблица 4 
Максимальные значения упругих растягивающих напряжений 
и относительной деформации образцов из полиэтилена в соответствии 


с утвержденными температурами 


ОТ ОЕ О ПОТ И ОЕ О ООО ОЕ О ОЕ ООС 


К. 29,576 | 26,074 | 23,656 | 22,340 | 19,716 | 18,426 | 17,229 
напряжений, МПа 





Средние значения относительной 0.0218 | 0.0226 | 0.0256 | 0.0275 0,0293 0,0306 0,0319 
деформации, % 


Таблица 5 
Максимальные значения изменения упругих растягивающих напряжений для образцов из полиэтилена в 
зависимости от изменения температуры в течение года 


Упругие напряжения, МПа 28,16 31,30 26,49 26,01 24/71 24,64 


Упругие напряжения, МПа 23,32 22,02 24,20 24,27 25,40 27.51 
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Рис. 3. Напряжения при растяжении полиэтилена при повышении температуры 
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Рис. 4. Значения напряжений при растяжении образцов из полиэтилена в течение года 


Результаты. Результаты обрабатывались по закону Гука следующим образом: 
#=Ар/Р=о/Ё, (5) 


где & — относительная деформация; АР — изменение диаметра трубы, мм; Р — диаметр трубы, мм; 
о’ — напряжения растяжения в стенке трубы, МПа; Ё — модуль Юнга, МПа; 

В таблице 6 приведены значения модуля упругости полиэтилена РЕ 100 и скорости распространения 
гидравлической ударной волны при при выбранных температурах. 

На рисунке 5 показано уменьшение значения модуля упругости полиэтилена с увеличением 
температуры, а на рис. 6 — уменьшение значения скорости распространения гидравлической ударной волны с 
увеличением температуры. 

Таблица 6 
Значения модуля упругости полиэтилена РЕ 100 и скорости распространения гидравлической ударной волны 


при утвержденных температурах 


ОС ТТ ОЕ ОТ ПОР О МОЕ О ОО ОЕ ОС 
а | в | 63 | 69 | 540 
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Скорость гидравлической 280 259 233 219 200 190 180 
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Рис. 5. Значения модуля упругости полиэтилена при повышении температуры 
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Рис. 6. Значения скорости распространения гидравлической ударной волны при повышении температуры 


Обсуждение и заключения 
На основании полученных данных, представленных на рис. 5, авторами выведено экспоненциальное 


уравнение для расчета значения коэффициента упругости полиэтилена как функции времени Ё =1,312е°”". 
Коэффициент корреляции для данного уравнения составляет А” = 0,988. 


На основании данных, представленных на рис. 6, получено экспоненциальное уравнение для расчета 


значения скорости распространения гидравлической ударной волны как функции времени С =275,9е°°". 


Коэффициент корреляции для данного уравнения составляет А” = 0,987. 


Проведя анализ результатов, полученных в ходе исследований, можно сделать следующие выводы: 

— высокая температура воды оказывает существенное влияние на свойства материала труб из 
полиэтилена высокой плотности РЕ 100 и на гидравлические свойства ударной волны; 

— значение модуля упругости полиэтилена высокой плотности РЕ 100 уменьшается с повышением 
температуры воды. Например, при повышении температуры от 4 до 60°С снижение значения модуля упругости 
достигает 60 %. Это приводит к значительному снижению скорости распространения гидравлической ударной 
волны. Величина снижения достигает порядка 35 %, что ведет к заметному снижению значения максимального 
изменения давления в результате гидравлического удара. 
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Введение. В статье проведен анализ применения композиционных материалов, как основного определяющего 
способа уменьшения массы планера летательного аппарата и беспилотного летательного аппарата. 
Перспективные методы неразрушающего контроля позволяют оценить техническое состояние данных 
материалов, а также с высокой точностью определить концентраторы напряжений на планере летательного 
аппарата и беспилотного летательного аппарата с целью принятия решения о дальнейшей эксплуатации 
данного объекта контроля. 

Целью исследования являлось повышение точности и оперативности оценки трещиностойкости 
композиционных материалов путем применения акустико-эмиссионного контроля. 

Материалы и методы. Приведена номенклатура применяемых композиционных материалов при 
строительстве различных летательных аппаратов, в том числе беспилотных, представлены наиболее возможные 
вероятные дефекты данных материалов вследствие влияния эксплуатационных факторов. Произведено 
сравнение применяемых методов неразрушающего контроля композиционного материала и выборка наиболее 
подходящего согласно специфическим преимуществам. Проведен эксперимент по определению пределов 
прочности углепластика при помощи аппаратно-программного комплекса акустико-эмиссионным методом. 
Представлены результаты исследования в виде рисунков, проецируемых аппаратно-программным комплексом. 
Результаты исследования. Описано применение акустико-эмиссионного метода контроля композиционного 
материала. 

Обсуждение и заключение. Полученные в результате проведенного эксперимента результаты могут быть 
использованы в процессе определения пределов прочности различных композиционных материалов акустико- 
эмиссионным методом неразрушающего контроля для проведения оценки технического состояния в 
машиностроении, судостроении и авиастроении. Статья рекомендована научным сотрудникам, занимающимся 
проектированием летательных аппаратов и беспилотных летательных аппаратов. 


Ключевые слова: композиционный материал, стеклопластик, боропластик, углепластик, органопластик, 
летательный аппарат, беспилотный летательный аппарат, неразрушающий контроль, акустико-эмиссионный 
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Введение. Развитие современной авиационной техники (АТ) сопровождается созданием новых 
конструкционных композитных материалов (КМ) с перспективными механическими и физическими 
свойствами. Уменьшение взлетного веса конструкции является неотъемлемой задачей для авиационной 
промышленности. Конструкция планера большей части современных воздушных судов (ВС) выполнена на 
53 % из КМ, а беспилотных летательных аппаратов (БПЛА) — на 90 %. Исходя из этого, есть необходимость в 
оценке технического состояния КМ планера радиационными, тепловыми, акустико-эмиссионными (АЭ) 
методами неразрушающего контроля (НК). Специфической особенностью метода акустической эмиссии 
является возможность оценки развития различных дефектов. 

Целью исследования являлось повышение точности и оперативности оценки трещиностойкости 
композиционных материалов за счет применения акустико-эмиссионного контроля. 

Композиционный материал — искусственно созданный человеком материал, получаемый соединением 
в одну структуру разнородных составляющих и характеризующийся лучшими свойствами, по сравнению со 
свойствами каждого из компонентов. 

КМ имеет следующие отличительные признаки: 

— запроектированные состав и форма; 
— невозможность встречи аналога в природе; 
— состав определяется набором п компонентов; 
— свойства КМ определяются свойствами компонентов; 
— «служебные свойства» — отсутствие таких свойств у компонентов по раздельности; 
— неоднородность в микромасштабе; 
— неоднородность в макромасштабе. 
Основные преимущества и недостатки КМ приведены в таблице 1: 
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Таблица 1 
Основные преимущества и недостатки КМ 


Преимущества Недостатки 


Низкая плотность (р = 1,35-—4,8 г/см?) Сложная технология получения 


Высокая прочность (6, = 1750 МПа) Высокая стоимость КМ 


Высокая жесткость (Е = 270 000 Мпа) Невысокая прочность связи волокон с матрицей 


Жаропрочность Низкое сопротивление эрозии и деструкция 
Термическая стабильность 


Основные материалы, применяемые в авиастроении, приведены на рис. 1 [1]. 





Композиционные 
материалы 
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стеклопластики 


Металлы, 
Углепластики армированные Боропластики 
волокнами 


Органопластики 





Рис. 1. Классификация КМ 


Материалы и методы 

Высокопрочные стеклопластики 

Стеклопластики — КМ, состоящие из полимерной матрицы и стекловолоконного наполнителя, в 
котором присутствуют газообразные включения. Объем использования стеклопластиков различного назначения 
увеличивается. Это обосновано невысокой стоимостью и доступностью сырья, малой энергоемкостью 
производства стеклянного волокна, универсальностью и возможностью регулирования в широких пределах 
физико-механических свойств, возможностью создания гибридных материалов и конструкций на их основе. 
Изделия из стеклопластиков поддаются всем видам механической обработки. К основным примерам 
стеклопластиков, применяемым при создании ВС и БПЛА, относят: КАСТ-В листовой, ВФТ-С, СТ-911-ТА, СК- 
9ФА, СТМ-Ф герметичный, ВПС-19М герметичный [2, 3]. 

Углепластики 

Углепластики — КМ, состоящие из связующего и упрочнителей в виде углеродных волокон, нитей, 
жгутов, лент или тканей. Связующими выступают различные эпоксидные или фенолоформальдегидные смолы. 
Основные марки углепластиков: КМУ-1, КМУ-1У, КМУ-1В, КМУ-1Л, КМУ-2Л, КМУ-ЗЛ, КМУ-3, КМУ-4Л, 
КМУ-4Э, КМУ-6-41. 

Боропластики 

Боропластики — КМ, состоящие из полимерного связующего и упрочнителя — борных волокон. В 
качестве матрицы используют эпоксидные и полиамидные связующие, а в роли упрочнителей — борные нити 
или комплексные боростеклонити. Применение боростеклонитей облегчает технологический процесс 
изготовления боропластиков. К ним относят: КМБ-1, КМБ-1М, КМБ-1К, КМБ-2К, КМБ-ЗК. 

Органопластики 

Органопластики — КМ, в которых в качестве армирующего наполнителя выступают волокнистые 
наполнители. Синтетические волокна обладают хорошими текстильными свойствами. Из них можно получить 
широкий спектр различных структур: нити, жгуты, ленты, комбинированные ткани. Синтетические волокна 
незначительно теряют прочность при текстильной переработке. Они малочувствительны к повреждениям. 
Примеры органопластиков: 7Т, 7ТО, 5Т, ЭТ, 6ТКС, 6ТКЬ, 7ТКС, 8ТКС. 
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Металлы, армированные волокнами 

Металлы, армированные волокнами, — КМ, в которых в качестве упрочнителя выступают волокна 
бора, углеродные волокна, нитевидные кристаллы тугоплавких соединений. В качестве матрицы применяют 
различные металлы и сплавы, обладающие характерной пластичностью. 

Примеры таких материалов: алюминий-борное волокно (ВКА-1А), алюминий-угольное волокно 
(ВКУ-1[), магний-борное волокно (ВКМ-[), магний- угольное волокно, никель-фольфрамовая проволока 
(ВКН-1). 

Объемы использованных КМ в конструкции планера перспективного ЛА различны и составляют: 
крыло — 80%, оперение — 81%, фюзеляж — 31%, пилон — 34%, шасси — 23 %. Органы управления 
самолета Ил-96-300, лопасти несущего винта вертолета Ми-28, планер самолета МС-21 частично выполнены из 
КМ. 

Объем КМ планера БПЛА достигает 90 %. Яркими примерами являются «Крунк», «Дозор-600», 
«Иноходец». 

Применение КМ в планере ВС и БПЛА может значительно облегчить вес конструкции. По формулам 
(1), (2) можно определить изменение массы планера при применении в нем КМ: 





АтТил — Тил — Тил — Пе Е Ат) — (УРатти, ео = фЕ1 — =р))ттл, + Дт (1) 
Дтил = а фЕ(1 — =) Тил; (2) 
где ти„ — масса 1-й части планера из традиционных материалов; ти, — масса 1-й части планера из 
композиционных материалов; т — масса традиционных материалов; ф; — коэффициент, учитывающий 
отношение масс традиционных и композиционных материалов, значения которого лежат в диапазоне 
0 =; < 1; &; — коэффициент массового достоинства КМ, равный отношению массы элементов 1-Й части 
к 
планера из КМ к массе 1-й части из традиционных материалов (&; = о. 
ПЛ 
Если Ати > 0, то тк, < тл„ — планер из КМ легче планера из традиционных материалов. Применяя 


формулы (1) и (2) можно добиться снижения массы конструкции планера ВС на 28 %, ЫШТА — на 60 %. 
Вследствие влияния эксплуатационных факторов (рис. 2) возможно возникновение дефектов в КМ, 
представленных в таблице 2 [1]. 


Эксплуатационные 
нагрузки 


Факторы 


Природно- 
климатические 


Биологические : 
птицы, грызуны, 


микроорганизмы 


Климатические: Аэродинамические 
температура, механические, ударные, 


влажность, вибрационные, температурные 


атмосферные осадки акустические 





Рис. 2. Классификация эксплуатационных факторов влияния на структуру КМ [1] 
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Таблица 2 
Классификация и характеристика дефектов КМ 


Наименование 
дефекта 


Изображение дефекта Причины возникновения дефекта 


Поверхностная трещина 


Высокая скорость 
охлаждения, превышение допустимых 
Внутренняя трещина 

механических нагрузок, 


транспортировка 


Сквозная трещина 





Повышение содержания летучих в 
препреге. изменение скорости 
Раковина РетрУ® Р 
нагрева, увеличение времени и 


величины приложения давления 





Попадание инородных материалов 
Поверхностное 
при изготовлении препрега при его 





инородное включение 
расколе и выкладке 





Данные дефекты можно диагностировать за счет применения методов неразрушающего контроля, 
позволяющих определить техническое состояние планера ВС или БПЛА без уменьшения их пригодности к 
применению. 

К таким методам относят: 

— рентгеновский метод; 

— тепловой метод; 

— ультразвуковой метод; 

— метод акустической эмиссии (АЭ). 

Рентгеновский метод 

Рентгеновский метод основан на приеме и анализе проникающего ионизирующего излучения после 
взаимодействия с контролируемым объектом. Интенсивность излучения изменяется в зависимости от 
плотности материала, его толщины и наличия дефектов'. Для регистрации прошедшего излучения применяются 
рентгеновские пленки выбранной вариационной чувствительности, флуоресцирующие экраны и телевизионные 
установки. Часть конструкции планера ЛА или БПЛА помещается в рентгеновскую установку «СаФайр» для 
осуществления контроля. Результат проведения рентгеновского излучения проецируется на пленке, в которой 
дефект будет изображен цветом бледнее, чем фон всей конструкции. Основным недостатком рентгеновского 
метода является вредное воздействие на объект и субъект контроля [4, 5]. 

Тепловой метод 

Тепловой метод основан на приеме данных 0б изменении тепловых температурных полей 
контролируемых объектов, вызванных дефектами. Данный метод осуществляется посредством применения 
тепловизоров и авиационных фенов. 

Первоначально происходит нагрев 1/3 части конструкции планера БПЛА при помощи авиационного 
фена. На межмолекулярном уровне происходит передача энергии от более нагретой части планера к более 
холодной. Наблюдаем процесс перетекания энергии используя тепловизор. Окраска дефекта отличается от 
описываемого выше процесса (рис. 3) [6-8]. 


'Феррозондовый метод неразрушающего контроля деталей вагонов. Руководящий документ РД 32.149. Москва: Федеральное агентство по 


техническому регулированию и метрологии. 2008. 163 с. 
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Рис. 3. Изображение дефекта планера БПЛА путем применения теплового метода 


Недостаток данного метода — неспособность осуществления проведения контроля в необогреваемых 
помещениях. 

Ультразвуковой метод 

Ультразвуковой метод основан на регистрации упругих волн, создавшихся в контролируемой детали. 
Дефектоскопы УД2В-П, УДТ-40, УСД-50 позволяют обнаруживать трещины порядка 0,5-—1 мм. Для 
осуществления контроля необходимо установление контактной среды между контролируемой частью планера 
ЛА или БПЛА и приемником-преобразователем дефектоскопа [9-11]. Основной недостаток — вредное 
воздействие на субъект контроля. 

Метод акустической эмиссии (АЭ) 

Метод АЭ основан на процессе излучения акустических волн в процессе перестройки структуры 
материала. Данный метод является пассивным, так как основан на регистрации низкочастотных импульсов 
дефектов деталей под нагружением. АЭ метод, в сравнении с традиционными методами НК, обнаруживает 
глубоко залегающие трещины в конструкции материала (таблица 3). 

Таблица 3 
Сравнение методов НК 


Метод НК Длина обнаруживаемой трещины, мм 


000000 
00.001 





Для диагностики был создан экспериментальный образец аппаратно-программного комплекса АЭ 
диагностики, который позволяет в реальном масштабе времени регистрировать развивающиеся дефекты КМ в 
низкочастотном диапазоне (рис. 4). Для повышения эффективности АЭ контроля специалистами Военного 
учебно-научного центра Военно-воздушных сил «Военно-воздушная академия имени профессора 
Н. Е. Жуковского и Ю. А. Гагарина» (ВУНЦ ВВС «ВВА»), в число которых входил и автор данной статьи 
А. В. Попов, был разработан теоретико-вероятностный подход к оценке информативных параметров эмиссии. 
Установлено, что на ранних стадиях деформирования поток сигналов АЭ КМ от микродефектов, случайным 
образом распределённых по объему ВС, имеет пуассоновский характер. С ростом нагрузки объединение 
микродефектов в трещину КМ нарушает пуассоновское распределение. Разработанный специалистами центра 
способ оценки процессов накопления повреждений КМ в ЛА и БПЛА основан на оценке изменения 
распределений числа актов АЭ на фиксированных интервалах времени в процессе деформирования КМ. 
Деформация КМ приводит к образованию макродефекта, характеристики потока импульсов становятся 
зависимыми. Объединение микродефектов ОК разрушает гипотезу пуассоновского распределения. Данное 
явление позволяет построить параметрические инварианты, справедливые для пуассоновского характера, и на 


этой основе оценить отклонение анализируемого процесса от пуассоновского (3): 
= Зт?[х|-т“[| =1 (3) 


На основании уравнения (3) получим несколько выражений для определения степени отклонения 
потока импульсов АЭ от пуассоновского (4)—(7): 
Ь = т[х3] *т?[х] — Зтй[х| — тёе[х| — т?[х| =1 (4) 
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В = т[х3] — 3т? [х]| — т[х| -1=0 (5) 
1 = т[х?] * т? [х](3 — п" [х]) =1 (6) 
те а + т?[х] =3 (7) 


Построение инвариантных зависимостей на основании формул (3)—(7) производит экспериментально 
разработанный АПК (рис. 4), позволяющий осуществить многоканальную регистрацию, обработку и анализ 
значений нагрузок и деформаций; осуществлять видеонаблюдение и запись процесса испытаний и показаний 


аналоговых приборов при помощи видеокамер. 





Рис. 4. Аппаратно-программный комплекс акустико-эмиссионной диагностики: | — высокочувствительный 
пьезоэлектрический датчик ОТ-300; 2 — предусилители виброакустического сигнала; 3 — аналого-цифровой 
преобразователь; 4 — ПЭВМ для обработки данных 


Высокочувствительные пьезоэлектрические датчики ОТ-300 устанавливаются на контролируемую 
поверхность КМ. Информативные параметры о состоянии КМ с пьезоэлектрических датчиков поступают на 
предусилители виброакустического сигнала °. После усиления сигнал поступает на аналого-цифровой 
преобразователь, в котором происходит его преобразование из механической работы в электрическую. ПЭВМ 
служит для кластеризации и обработки принятых данных. 

Результаты исследования. При помощи АПК и разрывной машины РМ-1 на базе ВУНЦ ВВС «ВВА» 
производили контроль углепластика КМУ-1, применяемого в производстве самолета ИЛ-96-300 и БПЛА 


«Орион». Основные характеристики КМУ-1 представлены в таблице 4: 
Таблица 4 


Характеристика углепластика КМУ-1 


Предел растяжения 0..ст= 1020 МПа 


Предел сжатия 0.„.т= 400 МПа 
Предел изгиба о„.‚= 1100 МПа 


Предел сдвига осдвиг= 30 МПа 


Модуль Юнга изгиба Е„.‚= 145х10_3МПа 


Модуль Юнга сдвига Е двиг= 3,5х10-ЗМПа 





Модуль Юнга растяжения Ераст= 180х10 МПа 


Для проведения испытаний поместили КМУ-1| в разрывную машину, предварительно установив на КМ 
датчики АЭ ОТ-300. Производили нагружение КМУ-1 при помощи АПК. Результаты нагружения представлены 


на рис. 5 (а, 6, в): 
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Рис. 5. Результаты нагружения КМК-1: а) осциллограмма сигналов АЭ КМУ-1; 6) режим работы АПК; 
в) проекция динамики нагружения РМ-1 КМ КМУ-1 


На рис. 5 а изображена осциллограмма сигналов АЭ КМУ-1. Максимальное значение амплитуды АЭ 
ОК соответствует разрушению матрицы и волокон ОК ( Ижах = 0,25 В). 

На рис. 5 б представлен режим работы АПК. Здесь выделены три зоны — красная, желтая, зеленая. 
Красная зона характеризуется наличием ярко выраженных дефектов КМУ-1. Численное значение инварианта в 
пределах от 0,66 до 1 при { = 7-9 мин. Это свидетельство критически активного дефекта — стадия разрушения. 
Желтая зона — зона активного дефекта и стадии образования трещины (1=0,33—0,66; { = 4—7 мин). Зеленая — 
доброкачественная зона, характеризующая стадию микротрещин (1=0—0,33; += 9—4 мин) [12, 13]. 

На рис. 5 в представлена проекция динамики нагружения РМ-1 КМ КМУ-1[. Вследствие увеличения 
нагрузки (Р = 25,75,150 кгс/мм^) растет значение инварианта и амплитуды АЭ КМУ-1. 

Обсуждение и заключения. Применение КМ является перспективным методом оценки технического 
состояния при проектировании различных ВС и БПЛА. Одно из основополагающих преимуществ КМ — 
облегчение массы конструкции по сравнению с традиционными материалами. Контроль и оценку технического 
состояния ВС и БИЛА необходимо производить перспективными методами НК, способными выявить 
различного рода дефекты на ранней стадии их развития. Этим требованиям отвечает метод акустической 
эмиссии. В перспективе можно установить экспериментальный образец аппаратно-программного комплекса на 
многофункциональный истребитель пятого поколения Су-57 для осуществления диагностирования 
конструкции летательного аппарата в полете. Данная акустико-эмиссионная система позволит 
идентифицировать дефекты, возникающие в элементах конструкции воздушного судна на малой стадии их 
развития. Экспериментально разработанный образец аппаратно-программного комплекса позволяет оперативно 
в реальном масштабе времени регистрировать развивающиеся дефекты авиационных материалов, применяемых 
в конструкции летательных аппаратах и беспилотных летательных аппаратах. 
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Введение. Статья посвящена анализу конструктивных решений гидравлических делителей потока, 
используемых для синхронизации гидравлических приводов рабочих органов технологических и мобильных 
машин. Выявлены потребности рынка в многопоточных беззолотниковых дроссельных делителях потока с 
управляемым коэффициентом деления, таких как многоосные шасси транспортных средств. 

Целями работы были анализ возможности и обоснование целесообразности разработки дроссельного 
четырехпоточного делителя потока беззолотникового типа с чувствительными элементами типа трубка 
Вентури. Решение должно обеспечить синхронность перемещения (вращения) более трех рабочих органов 
технологических и мобильных машин. 

Материалы и методы. Исследованы патенты существующих конструкций гидравлических делителей потока. 
Рассмотрены системы, в которых требуется деление потока гидравлической жидкости более чем на два 
исполнительных органа. Для конструкции дроссельного трехканального делителя потока беззолотникового 
типа предложен вариант модернизации, который позволяет делить поток на четыре ветви. 

Результаты исследования. Обоснована актуальность разработки многопоточного дроссельного делителя 
потока беззолотникового типа для применения в промышленных и мобильных машинах. Рассмотрены два типа 
четырехпоточных делителей, указаны их слабые стороны. Отмечено, что разработка многопоточного 
дроссельного делителя потока на основе конструкций, созданных в 1989 и 1991 годах, позволит сократить 
число гидронасосов, избавиться от последовательного соединения двухпоточных делителей. Таким способом 
можно снизить потери давления в гидросистеме и реализовать адаптивное управление гидродвигателями 
многодвигательных мобильных машин. Показана возможность на базе трехпоточного дроссельного делителя 
получить делитель-сумматор на четыре потока путем добавления выпускной камеры, связанной с камерой 
мембраны через канал, входящий в сопло трубки Вентури. Описан принцип работы такого оборудования. 
Обсуждение и заключения. Рассмотрены принципы построения дроссельных делителей потока 
беззолотникового типа. Предложен вариант модернизации для увеличения количества каналов деления с трех 
до четырех. Однако для подтверждения работоспособности такой конструкции необходим численный анализ 
различных режимов работы делителя — расчет приведенной объемной жесткости его рабочих полостей. 
Полученная информация может быть использована для модернизации гидравлических узлов технологических и 
мобильных машин, повышения их надежности, технологичности и коэффициента полезного действия. 
Определены вопросы, которые необходимо решить в дальнейших исследованиях. 


Ключевые слова: дроссельный делитель потока, трубка Вентури, гидравлическая трансмиссия, многоосное 
шасси, транспортная платформа, синхронизация. 
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Введение. Условие абсолютной синхронизации — соблюдение пропорциональности линейных и 
угловых перемещений гидродвигателей и всех производных от них по времени. Часто бывает необходимо 
синхронизировать: 

— однородные линейные или вращательные движения двух и более гидроприводов, 

— один двигатель с поворотом второго. 

В станкостроении и машиностроении при синхронизации приводов важно учитывать ошибки 
синхронизации по положению. В статье рассмотрены способы синхронизации гидродвигателей в 
многодвигательных системах с делителем потока, а также примеры применения таких систем в 
промышленности. Помимо этого обсуждается применение многопоточного дроссельного делителя потока 
беззолотникового типа в качестве синхронизирующего устройства многоосной колесной сцепки 
гидравлической транспортной платформы. 

Фундаментальные труды о гидромеханических системах известны с древности. Самый ранний из 
сохранившихся — трактат Архимеда «О плавании тел». Сегодня применение сложных многоуровневых и 
разветвленных гидросистем достаточно исследовано, однако остается актуальным вопрос управления и 
регулирования в гидроприводе" . [1—4] для повышения его КПД и уменьшения габаритов системы. 


' Урекин В. С., Истомина Ю. В. Разработка гидравлических систем автоматического управления приводами технологического оборудования // Инновации 
технических решений в машиностроении и транспорте : сб. ст. П Всерос. науч.-тех. конф. Пенза : Межотраслевой науч.-инф. центр, 2016. С. 259—262. 

 Лазута И. В., Лазута Е. Ф. Нелинейная корреляция релейной автоматической системы управления гидравлическим приводом // Образование. Транспорт. 
Инновации. Строительство : сб. науч. тр. нац. науч.-практ. конф. Омск : СибАДИ, 2018. С. 64—69. 
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В разряд сложных гидросистем входят такие, где один гидронасос питает два или более независимых 
гидродвигателя. При этом независимо от нагрузки на выходные звенья должна обеспечиваться синхронность их 
движений. Требуемая точность синхронизации определяет выбор ее способа. Простейший из них — 
обеспечение жесткой механической связи гидродвигателей. Однако для систем со сложным расположением 
исполнительных органов в пространстве он не применим. Наиболее распространенные методы: 

— с помощью связанных валом гидромоторов или гидронасосов (объемная синхронизация); 

— с последовательным соединением полостей гидродвигателей; 

— с гидравлическими следящими системами с распределителями и гидротахометрами, 

— с гидромеханическими и электрогидравлическими следящими системами. 

Наиболее целесообразные и простые средства синхронизации в мобильных машинах и 
технологическом оборудовании — это дроссельные делители и сумматоры потоков” ” [5-8]. Гидравлический 
делитель потока разделяет один поток рабочей жидкости на два или более. При этом учитывается заданное 
соотношение и требуется обеспечить синхронное (согласованное) движение выходных звеньев независимо от 
их нагрузки. Сумматор потока действует аналогично делителю, но в обратном направлении. Он легко 
интегрируется в гидравлическую систему, дешев и удобен в эксплуатации. Сумматоры потока обеспечивают 
синхронизацию в широком диапазоне скоростей движения выходных звеньев гидродвигателей. 

В классической конструкции дроссельного делителя в качестве запорно-регулирующего элемента 
используются золотниковые пары [5]. Это обусловливает характерные недостатки таких делителей. Например, 
высокая себестоимость прецизионных золотниковых пар объясняется сложностью изготовления, большой 
долей таких высокоточных операций, как хонингование, шлифование, электроэрозионная обработка. 
Селективная сборка и доводка золотниковых пар делает их практически неремонтопригодными, что тоже 
повышает их эксплуатационную стоимость. Вышеизложенное также указывает на восприимчивость системы к 
качеству гидравлической жидкости, к которой предъявляют весьма высокие требования [6]. 

Известны попытки конструировать делители без золотниковых пар” ® ”°. Исследовались возможности 
улучшения характеристик золотниковых делителей. Следует отметить, что перечисленные конструкции 
делителей потока обеспечивают синхронизацию только двух исполнительных механизмов, и это значительно 
сужает сферу их применения. Цель настоящего исследования — поиск перспективной конструкции 
многопоточного дроссельного делителя, отвечающего требованиям эксплуатации мобильных машин. 

Материалы и методы 

Многодвигательные гидромеханические системы. Известны гидросистемы, в состав которых входит 
более двух рабочих органов (как реверсируемых, так и не реверсируемых). В зависимости от условий 
эксплуатации они требуют постоянной или кратковременной синхронности движения с постоянным или 
изменяемым соотношением скоростей. К таким системам, в частности, можно отнести: 

— ходовые гидроприводы мобильных машин, 

— вытягивающие прессы для крупного листового проката, 

— штамповочные станки для габаритных деталей, 

— группы гидродомкратов для подъема фундаментов зданий, 

— гидроприводы складывания рулевого оперения в аэрокосмической технике, 
— системы позиционирования лепестков рефлекторного зеркала радиотелескопа. 

Применение в таких системах известных дроссельных делителей потока беззолотникового типа, 
несомненно, улучшает эксплуатационные характеристики. В то же время число синхронизируемых рабочих 
органов ограничено двумя на один делитель потока. Данное ограничение вынуждает конструкторов 
использовать последовательное каскадирование делителей потока. Вследствие этого на выходе из каскадов 
деления растут потери давления (до 30%) и суммируется ошибка синхронизации каждого делителя. 
Немаловажно и увеличение объема, занимаемого гидроаппаратурой. В целом, такой подход технически не 
оправдан и затратен в финансовом отношении. 


°Яцухин Ю. А. Исследование гидравлической системы синхронизации движения рабочих органов сельскохозяйственных машин : автореф. 
дис. ... канд. тех. наук. Ростов н/Д, 1972. 25 с. 

* Рыбак А. Т. Изыскание рационального типа запорно-регулирующего элемента дроссельного делителя потоков гидроприводов синхронных 
механизмов сельскохозяйственных машин : автореф. дис ... канд. тех. наук. Ростов н/Д, 1989. 23 с. 

° Делитель потока : патент 1041773 СССР : Е 15 В 11/22 / Р. П. Кириков, Э. Б. Шерман, В. Л. Кузик [и др.|. № 3336116/25-06 ; заявл. 
03.09.1981 ; опубл. 15.09.1983, Бюл. № 34. 3 с. 

° Делитель расхода : патент 1566897 Франция / № 150410 ; заявл. 03.05.1968 ; опубл. 1970. 

7 Дроссельный делитель потока : а. с. 1151725 СССР : Е 15 В 11/22 / А. Т. Рыбак, Ю. А. Яцухин, В. Н. Негодов [и др.]. № 3583155/25—06 ; 
заявл. 26.04.1983 ; опубл. 23.04.1985, Бюл. № 15. 3 с. 

*Ео\и агяаег : патент 10578221 США / №15/730,196 ; заявл. 11.10.2017 ; опубл. 03.05.2020. 
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В этой связи стоит отметить конструкцию дроссельного делителя потока мембранного типа”. Его 
можно использовать в реверсивных гидросистемах, в которых нужно динамично менять соотношение 
скоростей исполнительных органов. Кроме того, данная конструкция позволяет создать на ее основе делитель 
на 3, 4 и более независимых потоков. 

Рассмотрим применение дроссельного делителя потока беззолотникового типа на примере 
гидрообъемной трансмиссии колесного привода многоосных транспортных шасси. 

Для конструирования трансмиссионных систем современных полноприводных автомобилей, как 
правило, применяются блокированная и дифференциальная схема. Блокированная трансмиссия обеспечивает 
жесткую связь ведущих колес. Дифференциальная характеризуется разветвлением одного входящего потока 
мощности между ведущими мостами и колесами через специальный дифференциальный узел. Существуют и 
комбинированные трансмиссионные схемы, в которых одна часть колес объединена дифференциальной связью, 
а другая — блокированной (например, блокировка нескольких дифференциалов). Для трансмиссионных систем 
невозможна принудительная регулировка: 

— подводимой мощности к одному или нескольким колесам, 

— одного моста независимо от других. 

Наличие в трансмиссионной системе шасси гидропривода позволяет индивидуально подводить 
мощности к колесам и обеспечивает бесступенчатое регулирование в широком диапазоне значений. 
Конструкция гидропривода компактна, устойчива к внешним атмосферным воздействиям. К тому же у 
гидропривода, как правило, болыпой диапазон регулирования, что важно для автомобиля высокой 
проходимости. Еще одно преимущество гидрообъемного привода колес — возможность быстрого 
реверсирования. Это важно при прохождении препятствий путем раскачивания. Гидрообъемный привод колес 
обеспечивает длительное движение на минимальной скорости с высоким тяговым усилием, что также 
повышает проходимость автомобиля. 

Гидравлическая схема ходовой части для одной оси такой машины в общем виде представлена на 
рис. 1. 
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Рис. 1. Гидравлическая схема одной колесной оси мобильной транспортной машины 


Рассмотрим систему гидрообъемной трансмиссии. Двигатель внутреннего сгорания 6 передает 
крутящий момент на реверсивный гидронасос 5, который соединен трубопроводами и гидрораспределителем 4 
с делителем потока 1 и гидромоторами 2. Валы гидромоторов связаны с ведущими колесами. Зачастую 
аксиально-поршневые или радиально-поршневые гидродвигатели конструктивно соединены с колесом. Такая 
компоновка образует единый узел «гидромотор — колесо». Во время работы двигателя гидронасос 
обеспечивает гидромоторы ведущих колес гидродинамическим напором жидкости, который в дальнейшем 
преобразуется в механическую работу. В зависимости от конструктивного исполнения гидроагрегатов рабочее 
давление в системе находится в диапазоне 25—50 МПа. Эксплуатация узла «гидромотор — колесо» повышает 
проходимость шасси в результате непрерывного потока мощности и плавного изменения крутящего момента. 

Компании ЗЭстеиейе и Катаг (Германия) выпускают модульные гидравлические транспортные 
платформы: средние — На Сот@ К22 и тяжелые самоходные — ПиегСот@ К25 (рис. 2). 

Модули средних серий шириной 2750 и 3000 мм рассчитаны на нагрузку 23 и 25т на одну ось при 
скорости до 10 км/ч. Модули тяжелых серий оснащают 2—6-колесными осями. Ширина платформы — 3000 и 
3100 мм, нагрузка на одну ось — 36—45 т. Самоходные модули приводятся в движение силовыми установками 
Ро\уегВоо$ег. Варианты их объединения: 


” Дроссельный делитель потока : патент 1670191 СССР : Е 15 В 11/22 / Ю. А. Яцухин, Л. П. Колосов, А. Т. Рыбак [и др.]. № 3841241/29 ; 


заявл. 02.01.1985 ; опубл. 15.08.1991, Бюл. № 30. 3 с. 
° Многопоточный делитель-сумматор потоков : патент 1742530 СССР : Е 15 В 11/22 / А.Т. Рыбак, Ю. А. Яцухин, В. И. Антоненко. 
№ 4694140/29 ; заявл. 23.05.1989 ; опубл. 23.06.1992, Бюл. № 23. 2 с. 
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— одна платформа, на которую укладывают мостовое строение, 
— две тележки с опорно-поворотными столами. 
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Рис. 2. Два соединенных 6-осных самоходных модуля П\егСотё1 К25 компании Катае (Германия) 


В отечественном машиностроении можно выделить перспективную конструкцию многоосного, 
полностью гидравлического шасси, разработанную Центральным научно-исследовательским автомобильным и 
ы И 
автомоторным институтом «НАМИ». Это трехосный автомобиль «Гидроход-49061» (рис. 3). 





Рис. 3. Автомобиль НАМИ «Гидроход-49061» 


В конструкции этой машины реализована идея «гибкого колесного привода», которая предусматривает 
индивидуальный гидронасос для каждой из трех колесных осей. Транспортер построен на базе шасси «ЗИЛ- 
4906». Оснащение шасси гидрообъемной трансмиссией позволяет реализовать все возможные типы колесных 
трансмиссий: блокированный, дифференциальный, с индивидуальным регулированием — и их комбинации. 
Промышленность заинтересована в развитии производства объемно-гидравлических колесных шасси высокой 
проходимости, что стимулирует новые исследования в этой области. 

Делители потока широко применяются в схемах гидравлических трансмиссий для обеспечения 
блокирования дифференциала, что позволяет увеличить проходимость препятствий. Делитель потока 
гарантирует, что оба синхронизируемых колеса всегда будут создавать тяговое усилие, даже если одно из них 
находится на мягком или скользком грунте. Но для синхронизации и управления потоком в 4- и 6-осной 
конструкции необходимо последовательно соединить несколько делителей, причем с изменяемым 
соотношением деления по управляющему сигналу. 

Результаты исследования. Рассматриваемые четырехпоточные делители можно отнести к одному из 
двух типов. В первом случае речь идет о конструкциях с золотником. Это каскад из трех делителей потока в 


' Прочко Е. И. Методы построения систем силовых гидрообъемных приводов колес полноприводных автомобилей : автореф. дис. ... канд. 
тех. наук. М., 2007. 20 с. 
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12 
одном корпусе. В качестве примера можно привести делитель потока Росаш ЕО-М4 (Франция) ^. На рис. 4 
изображены его внешний вид и гидравлическая схема. 


А2 В2 
(2 





Рис. 4. Четырехпоточный делитель потока Роб1аш ЕО-4М: внешний вид (а), принципиальная гидравлическая схема (6) 


В конструкциях второго типа не предусмотрена возможность менять соотношение деления по 
внешнему управляющему сигналу. Это ограничивает сферу применения делителей в управляемых 
многодвигательных гидросистемах. Разработка многопоточного дроссельного делителя потока на основе 
конструкций, созданных в 1989 и 1991 годах“ г. позволит сократить число гидронасосов, избавиться от 
последовательного соединения двухпоточных делителей (соответственно, снизить потери давления в 
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гидросистеме), а также реализовать адаптивное управление гидродвигателями многодвигательных мобильных 
машин [77]. 

Так, например, можно получить делитель-сумматор на четыре потока (рис. 5) путем добавления 
выпускной камеры 16, связанной с камерой мембраны 10 через канал, входящий в сопло трубки Вентури 3. При 
этом левые камеры 15, 16 и правые камеры 17, 18 можно связать подвижным плунжером 13 между камерами 
1би 17. 
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Рис. 5. Четырехпоточный делитель с чувствительными элементами типа трубка Вентури 


Опишем принцип работы делителя-сумматора потоков с чувствительными элементами в виде трубок 
Вентури. Во время деления потока рабочая жидкость направляется на вход делителя-сумматора 1. Далее потоки 
разделяются между чувствительными элементами 2, 3, 4 и 5. Они выполнены в виде трубок Вентури, и их 
горловины соединены с прилегающими рабочими камерами 6, 7, 8 и9. При этом камеры 6, 7 и 8, 9 — 
исполнительные элементы конструкции. Они отделены друг от друга мембранами 10 и 11 и жестко закреплены 
в теле подвижных плунжеров 12 и 14. Подвижный плунжер 13 обеспечивает связь между левой и правой 
частями делителя. После трубок Вентури потоки жидкости проходят через осевые и радиальные проточки 
плунжеров 12 и [3, а далее попадают на переменные сопротивления 15, 16 и 17, 18 и к выходным каналам 19, 
20 и21, 22. 

Допустим, гидродвигатели на всех ветвях испытывают одинаковую нагрузку. Соответственно, в ветвях 
расходы одинаковы, а, значит, одинаковы и расходы через соответствующие им трубки Вентури. В горловинах 
трубок Вентури будут одинаковые скоростные напоры, поэтому давления в камерах управления будут также 
равны. С обеих сторон будет одинаковым давление на стенки мембранных элементов исполнительных 
механизмов, и это приведет их в уравновешенное состояние. 

В случае возрастания нагрузки (например, на гидродвигателе присоединенном к выходному каналу 19) 
сократится расход в соответствующей ветви и соединенной с ней трубкой Вентури 2. Как следствие, 
уменьшится скоростной напор в ее горловине, а, значит, давление в ней увеличится. В свою очередь, это 
приведет к повышению давления в рабочей камере 6. В результате нарушится равновесие мембраны 10, и она 
будет сдвигаться вправо. Плунжер 12 также станет перемещаться вправо, сужая живое сечение переменных 
сопротивлений 16 и 17, которые через подвижный плунжер 13 связаны с сопротивлением 16. В результате в 
соответствующих ветвях и трубках Вентури 3 и 4 расход уменьшится. Плунжер 14 начнет перемещаться вправо 
и уменьшит живое сечение переменного сопротивления 18. Этот процесс будет продолжаться до 
уравновешивания расходов во всех ветвях делителя. 

В режиме суммирования потоки для объединения подводятся к выходным каналам делителя- 
сумматора, и его рабочий процесс идет в обратном порядке. 
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Стоит отметить, что еще на стадии разработки необходимо знать надежность и функциональные 
свойства предложенной конструкции делителя-сумматора. Лучшие способы получить необходимые данные — 
моделирование и численные эксперименты ” [9—13]. 

Обсуждение и заключения. Широкое применение многодвигательных разветвленных гидравлических 
систем в промышленной и мобильной технике часто предполагает синхронизацию их исполнительных 
механизмов. В зависимости от условий эксплуатации и требований к функциональности используются решения 
на основе сложных конструкций из регулируемых объемных гидромашин или дроссельных делителей 
потока [5—8]. В случаях, когда необходимо задействовать механизм согласования на короткое время, наиболее 
экономичным, на наш взгляд, способом синхронизации является использование дроссельного делителя потока. 
Рассмотренные гидромашины (рис.2,3) в полной мере соответствуют понятию многодвигательной 
гидросистемы. Применение дроссельных делителей потока в конструкции ведущих колесных осей таких машин 
обеспечит реализацию дифференциальной схемы без механической трансмиссии. Это повысит проходимость и 
позволит эффективно управлять гидроприводами колесной пары. 

На рынке не представлены многопоточные дроссельные управляемые делители потока. Это указывает 
на востребованность их разработки. 

Дальнейшая работа по проектированию четырехпоточного делителя-сумматора потока описанной в 
статье конструкции будет направлена на создание математической модели. Ее базой станут основы теории 
объемной жесткости ® '” "8 [14—16]. Это позволит уточнить физические размеры элементов делителя-сумматора 
и рассчитать допустимые режимы работы. 
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Использование инструментов интеллектуального анализа данных и машинного обучения становится все более 
распространенным явлением. Их полезность особенно заметна в случае больших наборов данных, когда 
информация, которую необходимо найти, или новые взаимосвязи извлекаются из информационного шума. 
Развитие этих инструментов означает, что исследуются наборы данных с гораздо меньшим количеством 
записей, обычно связанных с конкретными явлениями. Такая специфика чаще всего приводит к невозможности 
увеличения количества случаев, а это может облегчить поиск зависимостей в изучаемых явлениях. В статье 
рассмотрены особенности применения выбранных инструментов к небольшим наборам данных. Предприняты 
попытки представить методы подготовки данных, методы расчета производительности инструментов с учетом 
специфики баз данных с небольшим количеством записей. Предложены избранные автором методики, которые 
помогли выйти из тупика в расчетах, т. е. получить результаты, намного хуже ожидаемых. Необходимость 
применения методов повышения точности прогнозов и точности классификации была вызвана небольшим 
количеством анализируемых данных. Эта статья не является обзором популярных методов машинного 
обучения и интеллектуального анализа данных, тем не менее собранный и представленный материал поможет 
читателю сократить путь к получению удовлетворительных результатов при применении описанных 
вычислительных методов. 


Ключевые слова: машинное обучение, интеллектуальный анализ данных, искусственные нейронные сети, 
ассоциативный анализ, автоматическая классификация. 
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Введение. В эпоху всеобщего доступа в Интернет все больше и больше устройств взаимодействуют 
друг с другом или с централизованными базами данных. Рекламодатели превосходят друг друга в 
эффективности персонализированной рекламы. Все это заставляет бурно развиваться группу инструментов, 
известную как искусственный интеллект. Объем данных, который необходимо обработать, чтобы получить 
нужную информацию, огромен, поэтому количество публикаций по алгоритмам, позволяющим быстро 
извлекать информацию из информационного шума, очень велико. Наиболее часто при этом приходится 
сталкиваться с информационной перегрузкой. Ученые из разных областей знаний знакомы с проблемами, 
связанными с анализом данных. Нередко сбор данных об изучаемых явлениях требует дорогостоящих 
устройств, установок и испытаний. Само исследование также может быть длительным. Это означает, что в 
научно-исследовательских базах данных о причинах и следствиях анализируемых явлений часто может 
содержаться всего несколько десятков или несколько сотен записей. Преимущества инструментов машинного 
обучения и интеллектуального анализа данных, в т. ч. возможности поиска значительных зависимостей между 
многомерными входными и выходными данными, дают возможность исследователям использовать эти 
инструменты для определения ранее не обнаруженных взаимосвязей изучаемых процессов и явлений. 
Недостаточное количество записей в созданной базе данных, описывающих какое-либо явление, может снизить 
ценность полученных результатов анализа. В статье представлены разработки автора, в которых инструменты 
машинного обучения и интеллектуального анализа данных использовались для исследования материалов и 
анализа процессов, когда количество входных данных было большим по сравнению с количеством 
выполненных тестов (т. е. записей в базе данных). Собранные примеры приложений были расширены за счет 
включения методов подготовки данных и методов оценки точности прогнозов и классификации, чтобы 
облегчить работу и помочь быстрее достичь ожидаемых результатов людям, которые намерены использовать 
инструменты машинного обучения для анализа собственных исследований. 

Анализ явлений, описываемых многими переменными. Перед любым исследователем непременно 
встают вопросы: какие входные значения принимать для анализа как влияющие на изучаемое явление, а какие 
параметры измерять на выходе. Очень полезен бывает статистический подход в исследовании, но он имеет 
существенный недостаток, который состоит в том, что можно анализировать только одну пару функций. При 
этом стоит, конечно, определиться, что такое статистика. Согласно [1], «Статистика — это наука о методах 
проведения статистического обследования и методах анализа его результатов». Предметом же статистического 
обследования являются отдельный набор — объектов, который — называется — статистическим 
сообществом (населением), или несколько статистических сообществ. Статистику можно разделить на три 
основные части: описательная статистика, распределение случайных величин и статистический 
вывод (рис. 1) [2]. 
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Рис. 1. Статистические управления [2] 
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В случае, когда на результат процесса влияет множество переменных, с помощью статистических 
методов очень сложно найти такие комбинации значений входных переменных, которые существенно влияют 
на изменчивость выходных данных. И самое востребованное — как эффективно управлять процессом или 
явлением, чтобы на выходе получить желаемый результат. С помощью инструментов интеллектуального 
анализа данных гораздо проще найти взаимосвязи между многомерными входными и выходными данными. 
Интеллектуальный анализ данных очень точно определяется самим названием книги [3], которое можно 
перевести как «Обнаружение знаний из данных». Существует определение интеллектуального анализа данных, 
сформулированное в 2001 году как анализ (часто огромных) наборов данных наблюдений с целью обнаружения 
неожиданных взаимосвязей и обобщения данных оригинальным способом, чтобы они были понятны и полезны 
своему владельцу [4—5]. Для этих нужд разрабатываются методы и алгоритмы, благодаря которым поиск 
вышеупомянутых соединений происходит быстрее и эффективнее. Методы интеллектуального анализа данных 
можно разделить на: 

— обнаружение ассоциаций (правила ассоциации); 

— классификация и прогнозирование; 

— группировка; 

— анализ последовательности и времени; 

— обнаружение характеристик; 

— интеллектуальный анализ текстовых и полуструктурных данных; 

— изучение контента, размещенного в Интернете; 

— изучение графиков и социальных сетей; 

— интеллектуальный анализ мультимедийных и пространственных данных; 

— обнаружение особенности [6]. 

На этой основе были разработаны методы, обычно называемые искусственным интеллектом, благодаря 
которым выполняются наиболее часто выбираемые задачи интеллектуального анализа данных. Несмотря на 
развитие информационных технологий и возрастающую вычислительную мощность компьютеров, до сих пор 
практически невозможно проверить все возможные комбинации многомерного ввода и вывода сложной 
системы [7]. Использование методов искусственного интеллекта тем более оправдано, чем сложнее проблема и 
не известны механизмы, управляющие ею, как показано на рис. 2. 
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Рис. 2. Предлагаемые условия использования методов искусственного интеллекта [7] 
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основе [8]. Вместо строгого поиска возможных комбинаций используется метаэвристика. Если необходимо 
использовать вышеупомянутые инструменты, следует решить, пользоваться ли специализированным 
программным обеспечением или создать его самостоятельно с помощью общедоступных модулей (так 
называемых «движков»), реализующих алгоритмы искусственного интеллекта. Независимо от принятого 
решения в основе будут лежать данные, которые будут анализироваться. 

Подготовка данных. Можно выделить следующие этапы подготовки данных к анализу: 

— очистка данных; 

— интеграция данных; 

— выбор данных; 

— консолидация и преобразование данных [6]. 

Такая подготовка должна выполняться независимо от размера базы данных. Для небольших наборов 
данных их правильная подготовка даже более важна, чем для больших. Примером может служить сравнение 
двух наборов данных: один с 10 000 записей, а другой со 100 записями, где 5 % записей относятся к 
повторяющемуся явлению (повторяемость еще не обнаружена). Когда две записи содержат ошибочные данные, 
то в первом случае можно найти повторяемость в 4,8 % случаев вместо 5,0 %. Во втором случае повторяемость 
обнаруживается только в 3,0 %. Разница существенная. 

Очистка и интеграция данных. При очистке данных из базы удаляются в основном записи, 
содержащие неполные данные. В больших базах данных удаление, например, двух записей существенно не 
повлияет на результаты, полученные на последующих этапах. При небольшом количестве наборов данных 
потеря даже одной записи может существенно повлиять на полученные результаты анализа. По этой причине 
отсутствующие значения не могут быть заменены, например, средним для всей генеральной совокупности 
(один из методов увеличения данных) или ее части (аналогичной описанию в записи, которую следует удалить), 
как это делается для больших баз данных Причина та же, что и описана выше — замена одного 
отсутствующего признака в описании явления может существенно изменить результаты, если анализируется 
небольшой набор данных. Однако записи, удаленные в процессе, не должны удаляться безвозвратно. На 
последующих этапах может оказаться, что в окончательно принятой модели данная особенность не будет 
учтена, а изначально удаленная запись будет содержать полные данные — это будет полезно для анализа. 

Второй важный этап очистки данных — это статистический анализ каждой характеристики (столбцов в 
базе данных) отдельно и ее корреляция с выходными данными. Рекомендуется представлять статистику 
основных характеристик анализируемого процесса (количество записей, среднее арифметическое, медиана, 
минимальное и максимальное значение, стандартное отклонение, квартили значений характеристик) также и 
для функции или функций, описывающих выходные данные. Диаграммы «рамка-усы» очень 


удобочитаемы (рис. 3). 
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Рис. 3. Пример диаграммы типа «кадр-усы» [9] 
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На таком графике легко прочитать, например, что для 50 % образцов с 21,9-25,0 % содержания глины с 
пылью прочность была выше 9 МПа, но при этом минимальная прочность для этого типа образцов была ниже 
2 МПа. Для таких образцов прочность была ниже 6,5 МПа. Анализ базовой статистики может облегчить 
решение об исключении из анализа записей (т. е. образцов или исследуемых явлений), для которых измеренные 
значения несовместимы со всеми другими случаями. Значительное несоответствие может быть результатом 
ошибочного измерения или того факта, что на измерение повлиял другой фактор, который вообще не 
принимался во внимание (он не учитывался и не измерялся ни в одном из случаев). По этим причинам все 
записи, отклоненные из базы данных, должны быть описаны, а также должны быть указаны причины 
отклонения [10]. 

Другой случай. Если обнаружено, к примеру, что решение об отклонении записи может быть принято 
только после того, как будут выполнены все или часть вычислений. Об этом идет речь в статье [11]: на основе 
наборов из 95 ускорений (а) стандартизированного молота (ударяющего по испытуемому стальному элементу), 
измеряемых каждые 0,01 мс с использованием искусственных нейронных сетей, была предпринята попытка 


отнести испытываемый стальной элемент к одному из девяти классов (рис. 4, 5). 
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Рис. 4. 95 ускорений молота для 67 испытаний стали [11] 


Анализируя ускорения на рис. 4, можно сказать, что один из тестов во временном диапазоне от 0,01 до 
0,31 мс отстает от других, но по-прежнему ведет себя как другие образцы. Только предварительная 
классификация по четырем группам марок стали показала, что через 0,31 мс в испытании № 29 были получены 
результаты, которые показывают резко отклоняющийся характер результатов также через 0,31 мс (рис. 5). В 
испытании №29 был исследован стальной образец, для которого во всех других испытаниях ускорение 
изменило знак с отрицательного на положительный между 0,424 и 0,450 мс. Для испытания № 29 знак 
ускорения изменился в пределах 0,460-—0,495 мс, то есть за время, подходящее для другой группы марок стали. 
Только этот вывод позволил достаточно хорошо обосновать отклонение из анализов испытания № 29. В 
результате была повышена первоначально полученная точность классификации до девяти марок стали по 
результатам 67 испытаний, равная 80 %, до 95 % (после отказа от теста № 29 и повторного использования 


искусственных нейронных сетей). 
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Рис. 5. Фрагмент диаграммы с рис. 4 с предварительной классификацией на четыре группы классов сталей и 
выпадающим тестом № 29 также через 0,31 мс [11] 


Интеграция данных — это объединение данных об одном и том же явлении из разных источников в 
одну базу данных. Пример интеграции содержится в работе [12], в которой прогнозируется задержка 
строительства участков скоростных и автомагистралей в Польше. Независимыми переменными, на основании 
которых делались прогнозы, являются данные о построенных объектах (предусмотренных законом о доступе к 
информации главного управления национальных дорог и автомагистралей), данные о предприятиях, 
реализующих эти объекты (собранные в регистрационный суд), Интернет, агентство бизнес-аналитики, 
макроэкономические данные (с источником в публикациях центрального статистического управления). 
Значение зависимой переменной (необходимое для «обучения» искусственной нейронной сети) — количество 
дней, на которое откладывается завершение каждой из проанализированных дорожных инвестиций, искали в 
публикациях в прессе ив Интернете. Собранная информация была использована для интеграции в базу данных 
о реализации 128 строительных проектов. В Польше в 2009-2013 годах было построено 156 участков 
скоростных и автомагистралей, но получить полную информацию о них было практически невозможно. После 
анализа хода этого строительства были отклонены и те случаи, когда возникали неожиданные нарушения 
(например, в виде протестов экологов, которые не учитывались в анализах как независимая переменная). Это 
уменьшило количество дел на 28, но обеспечило полноту, целостность базы данных — основы расчета. 

Выбор данных. В больших наборах данных их размер является существенной проблемой — большое 
количество записей приводит к неэффективной и длительной работе программного обеспечения. В базах 
данных с небольшими размерами записей программному обеспечению поиска отношений ввода-вывода может 
быть недостаточно, чтобы найти эти отношения. Бывает, что изучаемое явление можно описать многими 
параметрами, но в базе данных мало случаев (записей) с описанными параметрами явления. Таким образом, 
выбор данных означает необходимость выбора лишь нескольких независимых переменных, на основе которых 
будут выполняться классификация или прогнозирование выходного значения с использованием искусственного 
интеллекта (также известного как машинное обучение). При выборе независимых переменных может оказаться 
полезным следующее: 

— изучение взаимной корреляции линейных независимых переменных, а также корреляции со 
значениями на выходе; 

— анализ основных компонентов; 

— эмпирический поиск оптимального набора независимых переменных. 

Корреляционное исследование. Исследование линейной корреляции Пирсона между парами 
независимых переменных и между каждой из них и зависимой переменной может быть представлено в форме 
таблицы с числами, а также графически, в виде так называемых «тепловых карт» (пары независимых 
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переменных) [13]. Переменные наиболее сильно коррелируют положительно, а интенсивный синий цвет на 
рис.6 демонстрирует наименьшее значение коэффициента Пирсона. Сильная положительная или 
отрицательная корреляция, считанная с тепловой карты, не обязывает удалять переменную, сильно 
коррелированную с другой, это всего лишь предположение, потому что она сильно положительно коррелирует 
с 712 (коэффициент корреляции между ними составляет 0,88), и в то же время 715 не коррелирует с выходом 
(обозначено как \у, коэффициент корреляции равен 0,03). 
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Рис. 6. Примерная «тепловая карта» независимых переменных от п] до 05 и выходных данных 


Хотя вычисляется линейная корреляция, и фактическая связь между независимыми переменными (или 
независимой переменной и переменной, зависящей от выпуска) может быть не линейной, вычисление этих 
линейных корреляций часто подсказывает, какие переменные не включать (если есть необходимость их 
уменьшения). Такая проверка была сделана среди прочего в работах [9, 12]. В [12] количество зависимых 
переменных было сокращено, а в [9] для анализа была принята новая переменная как сумма значений двух 
сильно положительно коррелированных независимых переменных (это также было технически оправданно). 
Переменная, которая имеет сильную отрицательную корреляцию с другой независимой переменной, также 
может быть удалена из базы данных. 

Анализ главных компонентов. Анализ главных компонент (РСА) выполняется для независимых 
переменных — выходное значение не учитывается” В результате получаем рейтинг, показывающий, какая из 
независимых переменных больше всего влияет на изменчивость наборов независимых переменных. Можно 
представить, что каждая из независимых переменных осуществляет измерение многомерного пространства. 
Независимые переменные связаны между собой, они образуют множества (записи в базе данных, описывающие 
явление). Результатом РСА является ответ на вопрос, какая из независимых переменных наиболее ответственна 
за то, что расстояния (в многомерном пространстве) между точками (наборы независимых переменных, 
описанные в записях — координаты точек) являются наибольшими. Переменные, оказывающие наименьшее 
влияние на разброс данных, это те, которые могут быть удалены из анализа в попытке уменьшить количество 
независимых переменных. Примеры эффективного применения анализа главных компонентов для повышения 
производительности инструментов машинного обучения можно легко найти, например, в работах [14—16]. 
Однако стоит иметь в виду, что РСА не принимает во внимание значение зависимой переменной. 
Следовательно, нет уверенности в том, что именно независимая переменная, которая также оказывает 
наибольшее влияние на прогнозируемое значение (зависимая переменная), вызывает наибольшую 
изменчивость в наборах данных. 

Эмпирические исследования. И корреляционное исследование, и анализ главных компонентов не 
дают абсолютной уверенности в том, был ли выбор независимых переменных оптимальным. Оптимальный — 
это значит наиболее точный прогноз или максимально возможная доля точных классификаций для базы данных 
и выбранного инструмента машинного обучения. Инструменты искусственного интеллекта чаще всего 
применяются, когда их пользователь подозревает, что существует связь между вводом и выводом (между 
наборами независимых переменных и эффектом их совместного появления — зависимой переменной). Когда 
эти зависимости невозможно описать строго (функцией многих переменных), когда изучаемые процессы и 
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явления сложны (рис. 2), то использование инструментов машинного обучения может оказаться единственным 
способом, чтобы узнать об этом. Поэтому трудно ожидать, что какой-либо вспомогательный инструмент точно 
укажет, какие из зависимых переменных следует использовать для «обучения» искусственного интеллекта. 
Следовательно, одним из методов поиска оптимального набора входных данных (зависимых переменных) 
является эмпирическая проверка результатов инструмента искусственного интеллекта на различных наборах 
зависимых переменных. Есть два основных способа действия: вперед и назад. Вперед — это означает выбор 
двух зависимых переменных, на основе которых результаты прогноза или классификации являются 
наилучшими. Выбрать первый порой очень просто, сложно представить прогноз задержек без указания 
планируемой продолжительности [12, 17]. Чтобы выбрать вторую переменную, проверяем работу инструмента 
на каждой созданной паре зависимых переменных (их иногда называют предикторами) [18]. Когда выбрана 
лучшая пара предикторов, один из оставшихся предикторов добавляется последовательно. Это делается до тех 
пор, пока добавление какой-либо еще не использованной независимой переменной не улучшит результаты. В 
обратной процедуре первым шагом является использование всех независимых переменных, а затем 
последовательное удаление только одной, проверка того, какой предиктор был удален, точность прогноза и 
классификации увеличилась больше всего. Процедура продолжается до тех пор, пока удаление любого из 
предикторов не приведет к улучшению результатов. 

Консолидация и преобразование данных. Консолидация и преобразование данных состоит в том, 
чтобы они могли использоваться выбранным инструментом интеллектуального анализа данных [6]. Наиболее 
распространенной формой преобразования данных является их стандартизация, то есть такое преобразование 
значений независимых переменных и зависимой переменной, при котором они принимают значения из одного 
и того же диапазона. Стандартизация данных является результатом необходимости предоставить каждой из 
независимых переменных «равные условия для игры», которые будут включены в модель машинного обучения. 
В [12] дана формула: 

Чат << к а= (1) 


тахк (ао) 
где К — количество записей в базе данных; 
а; — 1-Й элемент переменной а до стандартизации; 
а; — 1-й элемент переменной а после стандартизации. 
Второй широко используемый тип стандартизации данных — это так называемая стандартизация «к 


нулевому среднему значению и стандартному отклонению единицы», определяемая следующей формулой: 
а =7— (2) 
где ау — среднее арифметическое значение переменной а до стандартизации; 

ау; — 1-й элемент переменной а до стандартизации; 

ал; — 1-й элемент переменной а после стандартизации; 

Обода — стандартное отклонение переменной а до стандартизации |7]. 

Другие типы стандартизации, также нелинейные, можно найти, например, в [19]. Однако следует 
помнить, что тип стандартизации данных может изменить результаты, полученные с помощью машинного 
обучения [20]. Таким образом, тип стандартизации данных может быть одним из параметров, которые 
инструмент настраивает для получения наилучших результатов. 

Бинаризация может быть вторым процессом преобразования данных. Она означает преобразование 
числовых значений переменной только в два значения (например, 0 и 1) по следующей формуле [22]: 

не 6 
а: >р>1 
где а; — 1-й элемент переменной а до бинаризации, 

ал; — 1-й элемент переменной а после бинаризации, 

р — параметр, выбранный пользователем. 

Бинаризация данных особенно полезна при поиске правил с использованием анализа корзины (также 
известного как анализ ассоциации)”. Этот тип анализа был создан для исследования содержимого корзин 
клиентов с целью увеличения продаж. Компьютерные программы с модулем анализа корзины работают 
наиболее эффективно, если переменные являются двоичными (данный товар присутствовал в корзинах 
покупателя или нет). Можно сформулировать множество научных задач, сама суть которых является двоичной, 
но в большинстве случаев описание явления включает числа, которые преобразуются в двоичную форму для 


применения анализа корзины [21—22]. Такое преобразование также может быть выполнено для переменной, 
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которая может принадлежать нескольким непересекающимся подмножествам. Затем два дихотомических 
подмножества (являющиеся суммами исходных подмножеств) создаются из первичных подмножеств. Тогда, 
если ао; принадлежит одному из них, то ал; равно 0, если другой равен 1. 

Измерения ошибок. До сих пор автором использовались общие термины, такие как «точность 
прогнозов», «правильность классификации», которые можно назвать качеством или эффективностью 
инструментов искусственного интеллекта. Однако, если проанализировать способы повышения качества их 
работы, то необходимо определить ошибки в результатах действия инструментов машинного обучения и итогах 
интеллектуального анализа данных. 

Ошибки прогноза. Инструменты машинного обучения в основном служат двум целям: для 
прогнозирования значений (регрессия) и для автоматической классификации. Использование одних и тех же 
общепринятых мер ошибок облегчает понимание работы, но также позволяет легче оценить ценность прогноза. 
Предполагается, что анализируется абсолютное значение ошибки. Следовательно, абсолютная ошибка 
(АЕ) [7, 18] может быть определена как 

АЕ = 6 -Ь| (4) 
где Б — прогнозируемое значение; 

Ь — фактически наблюдаемое значение. 

Относительная ошибка, выраженная как абсолютная ошибка в процентах (АРЕ), определяется как 


9 * 100% (5) 


АРЕ = |— 
Ь 
Чтобы иметь возможность оценить качество прогнозов, сделанных Сс помощью инструмента 





искусственного интеллекта (например, искусственной нейронной сети), некоторые данные не используются в 
процессе «обучения». После построения модели вводится этот набор данных, называемый проверочной 
выборкой, и машина делает прогнозы. Таким образом, прогнозируемые значения составляют дюжину, 
несколько дюжин или больше. Затем для оценки качества прогнозов можно рассчитать среднюю абсолютную 


о | 


МАРЕ = —® 


процентную ошибку (МАРЕ): 
ы-Ь; 
Ь: 
п 











100%) 
(6) 
где и — размер валидационной выборки. 

Наиболее распространенной мерой погрешности (проверочный тест) является среднеквадратичная 
ошибка (МЗЕ), определяемая как 


п г’_р 2 
МбЕ = ЕВ (7) 


При решении задач регрессии большинство инструментов машинного обучения благодаря 
использованию эвристических алгоритмов ищет отображение ввода и вывода, которое минимизирует МЪЗЕ. При 
указании качества полученных прогнозов наиболее распространенными являются МЗЕ или МАРЕ (или оба). 
Следует отметить, что в случае МАРЕ не имеет значения, рассчитывается эта ошибка для стандартизованных 
или реальных значений — МАРЕ то же самое. В случае с МЗЕ дело обстоит иначе. Эта ошибка чаще всего 
имеет разные значения для стандартизированных прогнозов и для прогнозов, преобразованных в истинные 
значения (без стандартизации). Следовательно, необходимо указать, для каких значений МЪУЕ были рассчитаны. 
Сравнение точности прогнозов (различных процессов, явлений с разными инструментами) на основе МЗЕ 
оправдано, если МЗЕ рассчитывается для стандартизованных значений. С другой стороны, с точки зрения 
практического применения важнее МАРЕ или максимальное значение АЕ. Полезность получаемых прогнозов 
также является важным вопросом [23]. Прогноз косвенных затрат на строительство со средней относительной 
погрешностью 6% можно считать очень точным и полезным прогнозом, но те же 6 % МАРЕ для прогнозов 
фондовой биржи делают их бесполезными [24]. Следовательно, размер ошибки, полученной в прогнозах, также 
следует оценивать с точки зрения полезности для лиц, принимающих решения, использующих прогнозы. 

Меры точности классификации. При использовании инструментов машинного обучения для 
автоматической классификации отнесение случая к неправильному классу можно оценить двумя способами. 
Во-первых, просто как ошибка. Однако одной информации о том, что инструмент правильно классифицирует 
90 % случаев, может оказаться недостаточно. Если есть несколько классов, которым назначены отдельные 
случаи (описанные в базе данных) (например, восемь), может случиться так, что для пяти классов 
классификация будет на 100 % правильной, а 10 % ошибок относятся к другим трем классам. Поэтому качество 
результатов классификации оценивается так называемой матрицей ошибок, пример которой представлен в 
таблице 1. 
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Таблица 1 
Результаты классификации тендерных процедур по валидационной выборке [25] 


Класс: Класс: 
в Класс: сговор Всего: для всех 
свободный от подозрение в ., 
очень вероятный классов 
сговора сговоре 


Численность в 


валидационной выборке 


Количество правильных 
классификаций 
Количество неправильных 
классификаций 
Доля правильных 
классификаций, % 


Доля неправильных 





классификаций, % 


Значительные различия в точности классификации отдельных подмножеств могут способствовать 
дальнейшему поиску еще более точной модели классификации. Матрица ошибок также может содержать 
информацию, к какому неправильному классу была неправильно отнесена данная запись из выборки проверки. 
Это влияет на вывод, основанный на прогнозах. Проанализируем пример из таблицы | со следующими 
предположениями: 

—[ два неправильно классифицированных делопроизводства из класса «свободный от сговора», 
отнесенных автоматическим классификатором к классу «подозрение в сговоре»; 

— четыре неправильно классифицированных случая из класса «подозрение в сговоре» были отнесены 
автоматическим классификатором к классу «весьма вероятный сговор»; 

— один неправильно классифицированный случай из класса «очень вероятный» был отнесен 
автоматическим классификатором к классу «подозрения в сговоре». 

В результате анализа можно констатировать, что до тех пор, пока классификатор не отнесет данную 
процедуру к классу «свободный от сговора», можно быть уверенным, что это производство не связано со 
стовором. Все дела, отнесенные к этому классу, были правильно классифицированы автоматическим 
классификатором (несмотря на точность классификации менее 100 %). Этот эффект использовался, например, 
в [11]. Следовательно, стоит проанализировать, к каким классам автоматически был отнесен данный случай. 

Ошибкой классификации в медицинских приложениях является разделение ошибок только на два 
класса, где так же важно не вводить лекарства здоровому человеку, как и не отказывать в лечении 
действительно больному человеку (принимая его за здорового) (рис. 7). 


Класс присваивается классификатором 
Положительный Отрицательный 


Количество Количество 
Положительный положительных верных, | отрицательных неверных, 
определенных, как ТР определенных, как ЕМ 
Отсортированный класс Количество Количество 


Отрицатель НЫЙ положительных отрицательных верных, 


неверных, определенных, определенных, как ТМ 
как ЕМ 





Рис. 7. Матрица ошибок классификации на два класса; серый фон указывает на правильную классификацию" [26] 
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Для п классифицированных случаев выполняется следующее равенство: 
п = ТР +ЕР +ТМ +ЕМ (8) 


Для интерпретации матрицы ошибок в форме, представленной на рис. 7, используются понятия 


5 
точности, прецизионности, чувствительности, специфичности, определяемые следующими уравнениями` [26]: 
ТР+ТМ 














точность = (9) 

прецизионность = —^ (10) 
реп _— ТР+ЕР 
ТР 

чувствительность = (11) 

ТР+ЕМ 

специфичность = — (12) 

ы — ЕР+ТМ 


Однако следует помнить, что указанные выше индикаторы применимы лишь в том случае, если 
инструмент классифицирует только два класса. 

Определение значимости обнаруженных ассоциативных правил. Благодаря использованию анализа 
корзины — одного из инструментов интеллектуального анализа — в данных обнаруживаются правила, которые 
можно записать как 

р > (13) 
где Ь — предшественник; й — преемник правила. 

Такое правило читается так: если был предшественник, то был и преемник. И предшественник, и 
последующий могут состоять из нескольких переменных, но наиболее часто ищутся правила, в которых 
предшественник описывается многими переменными, а последующий — одним (например, если давление 
упало утром, а температура в полдень превысила 30°С, потом днем была гроза). Такое правило (как в примере 
выше) не всегда работает. Следовательно, мерами качества обнаруженных правил являются не меры ошибок, а 
три параметра (пропорции) [6, 9], благодаря которым можно легко определить, что правило будет проверено в 
случае нарушения предшественника [3, 4, 6]: 

— поддержка — помечено как 5ир (от англ. зиррог); 

— уверенность — обозначается как сои} (от англ. сопЯ4епсе); 

— увеличение — отмеченный как ПИ (от англ. ПВ). 

Поддержка определяется следующим образом: 


вирь > В) =" 


14 
- (14) 
где п(Ь > А) — количество случаев, в которых возникновение предшественника сопровождалось появлением 
наследника; М — количество всех случаев в базе данных. 


С другой стороны, достоверность правила определяется следующим образом: 


сот! (Ь > В) = и 


(15) 
где п(Б)) — количество случаев, в которых было отмечено появление предшественника. 
Не менее важно приращение правила. Если его значение меньше 1, это означает, что найденное 


правило не объясняет появления наследника. Приращение определяется следующим образом: 


ПАЕСЬ > В) = "2 и 


(16) 
где Р(Й) — вероятность преемника (независимо от того, был ли предшественник или нет). 

Лучшему пониманию оценки качества обнаруженных ассоциативных правил будет способствовать 
пример, в котором в 12 процессах было сделано 10 наблюдений за возникновением предшественника и 
последующего (рис. 8). Для правила «если предшественник, то наследник» были рассчитаны поддержка, 


достоверность и приращение (для каждого из процессов), которые представлены в таблице 2. 


> РОЗаЕ ЗабузЕусхпе Оргозгато\уаше ОБИсхешо\е, там же. 
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№ № последовательного наблюдения 


пропосса ео 
ие ЕЕ 
поиожсныя МЕН ИИИ — 
явление предшественника ТТ 


2 
появление наследника 
явлению предшественника в БЕ БЕ вЫ 
З арааанив: Борак» |пранний ннаныия ЕЙ ЕЕЕЕЕЙ Е 
появление слдииа Г [Г | 
д появление предшественника 
появление наследника 
-.- .-`-.-`- `` = -[--- 
появление предшественника 
5 
появление наследника 
появление предшественника 
появление наследника 
появление предшественника 
7 
появление наследника 
появление предшественника 
8 
появление сльшие | [ГОГ 
появление предшественника 
9 
появление наследника 
появление предшественника 
10 
появление наследника 
— =>... чинно сн ею 9 ры са 35355-35: 
появление предшественника 
11 
появление наследника р? 
появление предшественника 
12 





Рис. 8. Наблюдения за появлением предшественника и преемника в 12 процессах 


Сравнивая процессы 5 и 6, следует отметить, что высокая степень достоверности правил не всегда 
важна. В процессеб преемник почти всегда присутствует, и обнаруженное правило не объясняет 
возникновение преемника (подъем<1). Это не относится к процессу 8. Приращение указывает на важность 
правила, в то время как его достоверность низкая. Тем не менее, каждое наблюдение за преемником 
сопровождается наблюдением за предшественником. В случае процесса 8, следовательно, стоит уточнить 
предшественник (например, добавив еще одну переменную). В этом случае, вероятно, другой параметр (пока не 
включенный в предшественник) влияет на наличие преемника (в процессе 8). 


Таблица 2 
Оценка правил с использованием зир, сопГи ПП для процессов на рис. 8 
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Важные аспекты применения выбранных инструментов искусственного интеллекта и 
интеллектуального анализа данных к небольшим базам данных 

Количество записей в базе данных и сложность инструмента. В статистике чаще всего 
предполагается, что небольшой размер выборки составляет менее 30 случаев, однако можно обнаружить, что 
предельное число, за пределами которого мы не можем говорить о малой выборке, составляет 100 [27—28]. 
Инструментам искусственного интеллекта нужны полные наборы данных (входы и выходы), чтобы иметь 
возможность найти наиболее точный способ преобразования одного в другой. Чем сложнее проблема, тем 
больше требуется наборов данных (записей в базе данных). В глубоком обучении, когда инструмент 
преподается не с наборами чисел, а с файлами (графикой, аудио, текстом), необходимы тысячи наборов 
данных. В [29] более 4000 стандартизированных изображений использовалось для прогнозирования прочности 
на сжатие. В других исследованиях образцов было немного (например, в [11] всего 66). Есть указание на то, что 
для искусственных нейронных сетей количество связей между нейронами должно быть в 10 раз меньше, чем 
количество записей в базе данных [30]. Небольшие базы данных требуют большего количества 
вычислительных испытаний и более точной настройки используемых инструментов. Тем не менее, правило, что 
чем сложнее проблема и ее модель, тем больше наборов данных требуется для обучения инструмента, остается 
в силе. Использование сложных моделей (например, искусственных нейронных сетей с более чем одним 
скрытым слоем и множеством нейронов в скрытых слоях) в небольших базах данных чаще всего приводит к 
ошибкам (прогнозирования или классификации), намного большим, чем в моделях с меньшей сложностью 
самого инструмента. Отсюда популярность методов уменьшения количества независимых переменных, 
описанных в разделе 3. Когда количество наборов данных слишком мало, уменьшение количества независимых 
переменных чаще всего повышает точность прогнозов и классификации. 

Тип вывода и качество прогноза и классификации. В [30] можно найти предположение, что 
искусственные нейронные сети могут более точно определить, будет ли, например, предсказанное значение 
больше, чем значение, данное пользователем сети. Ссылаясь на требование полезности прогноза [23], если 
полученные прогнозы недостаточно точны (т. е. ошибки прогноза слишком велики), можно решить, требуется 
ли точное значение. Например, при прогнозировании прочности материала возможно вместо прогноза лишь 
сообщить, что прочность не будет ниже предполагаемой прочности, указанной пользователем. То же самое и с 
проблемой классификации. Автоматическая классификация стали на 9 марок на основе всего 66 записей в базе 
данных не позволила получить точность классификации выше 80% [11]. Затем один процесс классификации 
был заменен восьмым, в результате чего результаты испытаний стали были разделены на две дихотомические 


подгруппы, как показано на рис. 9. 


Этап! Этап 2 Этап3 Этап 4 Этап5  Этапб  Этап7 Этап 8 
р 51 991 991 91 91 91 91 91 
: 52 5&2 9а2 942 9а2 992 93а? 942 
то аз 523 аз 9а3 9а3 аз 9а3 93 
9: 524 924 9а4 924 994 9а4 9а4 9а4 
Е 5а5 9а5 9а5 9&5 95 95 9а&5 9&5 
8 56 9а6 9&6 9а6 96 9а6 96 96 
Ее За7 927 997 927 927 97 97 9а7 
2 5а8 528  9а8 928 928 98 98 98 
ы 5а9 529 929 9а9 9а9 9а9 9а9 99 


Рис. 9. Восьмиэтапный процесс классификации на два дихотомических подмножества 
(на каждом этапе классификации) [11] 


Автоматическая классификация только двух подмножеств была проведена с помощью более простой 
модели (со сложностью, соответствующей количеству испытаний стали, т. е. 66). Использованный восемь раз 
процесс позволил повысить точность классификации с 80 до 95%. Из первоначально выбранных 
12 предикторов осталось только 6 [12]. Это позволило сделать прогноз с достаточно малой ошибкой, чтобы 


Хуберт Аныш. Применение инструментов машинного обучения и интеллектуальный анализ данных в отношении баз данных 





модель была полезной. Компьютерные программы позволяют прогнозировать, например, две зависимые 
переменные одновременно, но по причинам, описанным выше, лучшие результаты (меньшие ошибки) 
получаются при прогнозировании двух зависимых переменных отдельно, с помощью двух моделей, хотя 
результаты получаются на одном и том же наборе данных. 

Модификация вывода. Рассмотрим случай, когда на начальном этапе расчетов не получаются 
достаточно точные прогнозы или классификации. В дополнение к действиям, связанным с данными и 
сложностью в инструменте машинного обучения, использованном выше, можно рассмотреть возможность 
поиска другой зависимой переменной, той, на основе которой можно будет вычислить строго зависимую 
переменную, которую необходимо найти. При небольшом количестве записей в базе данных это может 
упростить модель, что, в свою очередь, может повысить точность прогнозов. Этот эффект был использован 
в [31], где прогнозирование с помощью искусственной нейронной сети является лишь частью процесса 
проектирования состава смеси. Еще одна процедура, которая может уменьшить ошибки прогнозов, — это 
прогноз относительного значения (вместо абсолютного значения). При оптимизации работы модели в [12] 
прогнозируемая задержка, выраженная в днях, была заменена задержкой, выраженной как пропорция 
количества дней задержки к запланированному количеству дней строительства. В указанном случае это не 
уменьшило ошибок в прогнозах. Другая возможность изменить тип вывода (т. е. прогнозируемую зависимую 
переменную) — заменить одно число несколькими значениями функций принадлежности, рассчитанных на 
основе теории нечетких множеств [32]. В [17] вместо количества дней задержки построения на выходе из 
искусственной нейронной сети использовались три значения функции принадлежности множествам: низкая 
задержка, средняя задержка и большая задержка. После уточнения прогнозируемых значений оказалось, что 
ошибки прогноза были меньше, чем при прогнозировании количества дней задержки [12]. То же самое было 
предпринято в [33] путем прогнозирования значений функции принадлежности на первом этапе вычислений, и 
только на втором этапе на основе этих прогнозов случаи были разделены на три подмножества. Однако в этом 
случае прямое использование искусственной нейронной сети в качестве классификатора привело к повышению 
точности классификации. 

Гибридные инструменты. При небольшом количестве кейсов в базе данных используемый 
инструмент машинного обучения не может быть очень сложным, поскольку слишком мало кейсов для 
успешного обучения модели. Гибридные модели могут быть средством от слишком большой ошибки прогноза 
или слишком низкой точности классификации. Вместо одного сложного инструмента используются два более 
простых. В вышеупомянутом примере [17] применение теории нечетких множеств фактически добавляет к 
модели два элемента, которые можно правильно «настроить». Схема модели представлена на рис. 10. 
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Рис. 10. Три модуля нейронечеткой модели [17] 


Первый модуль — это преобразование точных чисел в значения трех функций принадлежности. 
Второй — настроить сеть для прогнозирования этих трех значений с наименьшей ошибкой, а третий — 
повысить точность полученных прогнозов до точных цифр. Операции в каждом из этих трех модулей могут 
выполняться по-разному, поэтому для моделирования зависимой переменной можно использовать три 
инструмента, а не только саму искусственную нейронную сеть. Есть много примеров, когда гибридные модели, 
то есть те, в которых более одного инструмента используются совместно, дают более точные прогнозы, чем 
модели с одним инструментом [34—35]. Поэтому при анализе результатов исследований стоит рассмотреть 
возможность использования инструментов машинного обучения вместе с другими математическими 
инструментами. 
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Извлечение проверочного подмножества. Для обучения искусственных нейронных сетей из 
существующей базы данных выделяются три подмножества, содержащие как независимые, так и зависимые 
переменные: обучение, тестирование и проверка. Обучающее подмножество используется для обучения 
инструмента. Этот процесс продолжается до тех пор, пока МЗЕ не перестанет уменьшаться в тестовой выборке, 
затем процесс обучения сети останавливается. Дальнейшее обучение сети могло бы привести к лучшему 
согласованию инструмента с обучающими данными, но возможности обобщения были бы потеряны (ошибки 
МЗЕ для тестовой и проверочной выборки были бы намного выше) [30]. Такой эффект, называемый 
переобучением искусственной нейронной сети, схематично показан на рис. 12. 
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Рис. 12. Схематическое изображение переобучения искусственной нейронной сети 
и правильно обнаруженного тренда [30, 31] 


Подмножество валидации используется для оценки качества прогнозов или классификации. Эти 
наборы независимых переменных и зависимая переменная не видят сеть в процессе обучения. Поэтому 
прогнозы (или классификации) делаются для проверочного подмножества, и путем сравнения с известными 
зависимыми переменными могут быть вычислены ошибки, описанные в предыдущих разделах. Очень важно 
указать, к какому подмножеству относится вычисленная ошибка. Один из методов оценки того, не 
переоборудован ли инструмент, — это сравнение ошибок (чаще всего МЗЕ) для вышеупомянутых трех 
подмножеств, они должны быть на аналогичном уровне. Явно более низкая МЗЕ для обучающего 
подмножества может указывать на переобучение, а явно более высокая МЗЕ — на несовершенство 
инструмента. Явно более низкая МЕ для тестовой и проверочной выборки предполагает необходимость 
выбора других параметров искусственной нейронной сети или других зависимых переменных (или даже 
другого инструмента). 

Для больших баз данных подходит случайное разделение данных на три подмножества. В литературе 
встречаются предложения о том, что пропорции размера этих подмножеств должны быть в диапазоне 60:20:20— 
70:15:15 (преподаватель: тест: проверка) [36-37]. Для небольших баз данных случайный выбор подходит 
только в тех случаях, когда различные диапазоны зависимой переменной (или класса для классификации) 
одинаково численно представлены. Чаще всего это условие не выполняется. Тогда было бы хорошо обеспечить 
такое сбалансированное представление во всех трех подмножествах. Такая контролируемая разбивка данных 
использовалась в [11,25]. ВПП] из 66 испытаний стали одна марка (5 испытаний) была наименее 
многочисленной, а марка с 12 испытаниями — самой многочисленной. При случайном выборе тестов для 
подмножеств могло случиться так, что в обучающую подгруппу не было включено ни одного теста какой-либо 
марки стали, что, несомненно, помешало бы его автоматическому распознаванию. В [25] из 
249 проанализированных дел только девять, по мнению авторов, следует отнести к категории «весьма 
вероятный сговор». После раздела набора данных на подмножества четыре из этих процедур были назначены 
обучающими, одна — тестовой и четыре — проверочными. Благодаря этой процедуре три из четырех процедур 
из набора для валидации могут быть правильно классифицированы. 

В польской и англоязычной литературе можно найти примеры того, что авторы используют только 
концепции обучающего набора данных, тестового набора данных. Вероятно, это происходит по двум причинам. 
Некоторые компьютерные программы полностью (без вмешательства пользователя) контролируют, не 
переоборудован ли инструмент. Затем извлекается только проверочный набор (либо путем указания его 
количества или процентной доли в общих данных, либо путем выбора записей, которые будут использоваться 
для проверки). В этом случае подмножество проверки часто называют подмножеством тестов. Вторая причина 
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заключается в том, что некоторые инструменты машинного обучения (например деревья классификации, 
деревья С&ВТ°) защищены от переобучения другим способом, нежели контроль М$Е для тестовой выборки. 
Качество этих инструментов можно проверить с помощью тестового подмножества (точно так же, как 
подмножество валидации — не участвуя в процессе обучения инструмента). Обычно уже при чтении работ по 
машинному обучению становится ясно, были выделены два или три подмножества данных, поэтому различная 
номенклатура набора, используемого для оценки качества работы, не является проблемой. 

Проверка качества модели. Базы данных с небольшим количеством записей, используемые для 
построения модели на основе машинного обучения, делают модель менее устойчивой (т. е. дают существенно 
разные ошибки) при замене записей между подмножествами (обучение, тестирование и проверка). Чтобы 
узнать, работает ли построенная модель хорошо только с определенным разделением данных на подмножества, 
ее следует запустить на рандомизированных подмножествах. Используя предложения по пропорции разделения 
данных, их можно поделить на 5—7 подмножеств, и процессов построения модели должно быть выполнено как 
можно больше, так что на каждой итерации будет происходить обучение. Такая проверка качества модели 
называется перекрестной проверкой (на английском языке наборы данных при перекрестной проверке 
называются «свертками»). Ошибки МЗЕ или МАРЕ затем усредняются. Если, однако, в каком-либо наборе 
данных ошибки значительно отклоняются от среднего значения, следует поискать причину этого. 

Независимо от размера базы данных оценку ошибок также следует рассматривать через призму 
полезности результатов. Прогнозы с относительно большими ошибками, явления, которые нельзя точно 
описать, могут быть очень полезными и считаться важными. С другой стороны, прогнозы с теми же МЪЗЕ или 
МАРЕ другого явления могут не дать какой-либо новой информации об изучаемом явлении. Поэтому, помимо 
числовых значений точности прогноза или точности классификации, важно ссылаться на само изучаемое 
явление. 

Также важно проверить, не была ли построена другая модель для ранее изученного явления. Если да, то 
ссылка на эти предыдущие исследования (относительно уровня полученных там ошибок) также подтвердит 
качество нового патентованного решения. Если такие модели ранее не создавались, можно проверить качество 
нового решения, сравнив результаты с гораздо более простой моделью (например, на основе Мпсгозой Ехсе]| с 
надстройкой эо]уег). 

Существует вопрос, на который пока нет четкого ответа: означает ли небольшое количество наборов 
данных, что полученные прогнозы, классификации и правила ассоциации неактуальны (именно из-за их 
небольшого количества)? Например, сложное и строгое правило со следующими индексами 5ир = 0,01, 
сот] = 1, ПЕ = 100 для базы данных с 10000 записями означает, что каждый раз, когда конкретный 
предшественник встречается 100 раз, указанный предшественник всегда (5$1с!) является последующим. То же 
правило для базы данных из 100 записей означает, что был только один уникальный случай, когда указанный 
предшественник имел место. Затем его сопровождал конкретный преемник. В одном случае нельзя говорить о 
правиле, но можно о случае (случайности). Однако, даже для небольшой базы данных (например, со 
100 записями) правило 100-процентной достоверности, поддерживаемое тремя случаями, может уже указывать 
на повторяемость (при условии, что приращение для этого правила больше единицы). 

Существуют небольшие базы данных обычно потому, что они малы сами по себе, потому что 
подготовка более крупных дорогостоящая, требует очень много времени, они к тому же могут быть следствием 
природы явления (например, ограниченное количество объектов, построенных одной и той же компанией). В 
реальности эти базы данных не могут быть расширены (в несколько раз или в несколько десятков раз). 
Изучение таких баз данных при соблюдении условий, описанных в этой главе, может привести к появлению 
новых, пока не обнаруженных зависимостей. Построенные, правильно функционирующие модели могут быть 
неприменимы напрямую к другим подобным явлениям, но они могут эффективно указывать на методы поиска 
общих взаимосвязей в случаях, когда количество наборов данных значительно больше. 

Уникальные правила ассоциации или неожиданные автоматические классификации также могут 
указывать на области, на которых следует сосредоточить дальнейшие исследования описываемых явлений. 

Заключение. Обсуждаемые в статье проблемы и вопросы, связанные с поиском взаимосвязей между 
многомерными входными и выходными данными, представлены в небольшом количестве случаев. Данная 
работа поэтому не может являться исчерпывающим обзором по теме. Объем статьи и ограниченный опыт 
автора в применении методов машинного обучения и интеллектуального анализа данных не позволяют описать 
большинство используемых методов. Проблемы при расчетах (в основном на небольших базах данных), 
содержащиеся в упомянутых в тексте работах, и их решения систематизированы таким образом, чтобы 
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соблюдалась последовательность действий: от подготовки базы данных до расчетов, до обсуждения 
результатов. Невозможно четко указать, сколько записей в базе данных может свидетельствовать о «малой» 
или «большой» базе данных. Однако можно сказать, что для эффективного использования инструментов 
машинного обучения или интеллектуального анализа данных требуется как минимум несколько десятков 
записей. Применяя соответствующие процедуры (разработка входных данных, построение моделей), можно 
применять эти инструменты для успешного моделирования и изучения явлений, описанных всего в нескольких 
десятках случаев. На основе проведенных расчетов возможно сделать надежный вывод. 
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Гптгодисйоп. 'Та$ рарег за41ез фе сара иу о{ Чегет: бурез оГагайсла| пеига| пебмогКк$ (АММ) 0 ргед1с" Фе тоди[и$ 
ог @“азиспу оЁ рауетеи: 1ауегз ог Нехе азрвай рауетеп{ ипдег орегайио сопа1оп$. Тре Шпо уе деЙесцотащег 
(ЕУУО) аз з4еце4 10 зпишае фе дупапис га Ис |1оа45 ап теазиге фе Нехига! 6о\/1$ оп Фе гоа4 загсе 10 оМалш Ше 
даабазе ог АММ тоде{5. 

Маепта5 апа Мешо4. АтаЁс1а| пебмогК$ тбурез (Фе !ее4Ююг\ага ЪасКргорасайоп, 1ауег-геслигепь сазса4е Баск- 
ргорагайоп, апа Ейпап БасКргорагайоп) аге деуеоре ю 4ейпе Ше орита|! АММ тоде| изше МаЧаб зой\аге. То 
арргеслае Ше еРйслепсу оЁеуегу то4е|, \е изе4 е сопзиа‹еа АММ то4де!5 Юг ргед1сйпе Фе е]азис тоа[а$ уа[ез Юг 
25 пе\ рауетепе зесйопз фай \еге по{ изе4 ш Фе ргосе$$ оЁ гаш, уа|Чайоп, ог {езНпо фю епзиге #5 зацабИиу. ТВе 
еЁНстепсу шеазигез зисЬ аз шеап абзо[ще егог (МАЕ), ше соеЁЯсете оЁ шире деегишанопз В^, Воо{ Меап Запаге 
Егтог (КМЗЕ), Меап АБзопие Регсет Еггог (МАРЕ) уа[аез \еге оМаше4 Гог аП то4де[$ гези 6. 

Кезий5. Вазе4 оп Фе регюогтапсе рагатщегз, ( \уаз сопсТа4е4 Ша атопе Фезе а[2огИил$, фе Еее д-ГРог\уаг4 то4е] Ваз а 
Бецег регРюогтапсе сотраге4 {о Фе оШег гее АММ 1урез. ТВе гези $ ог Фе Бе5{ юг то4е|$ ууеге сотраге4 {о еасв оег 
апа фо Ше асша| дайа оМате4 о даегилте Пе Без те оч. 

О5сизяоп апа Сопсши5оп5. 'ТВе 4Негепсез Бебуееп Фе гези5 ог Фе юлпг Бе5ё по4е15 Гог Фе Юг фурез оЁ а! ог 115 
изе4 ууеге уегу зта|, аз беу зВо\е Ше с1озепезз Бебмееп ет ап4 фе асла| уашез. ТВе гезеагсВ гези $ сопйми Фе 
розу ор АММ-Базе4 по4е15 10 еуашае Фе @азис тоЧа[а$ оЁ рауетеп: 1ауегз зрее4Пу апа гейа у юг из И ш Фе 


угасбига| аззеззтепЕ о? (МОТ) Нехе рауетеп: даа а{ Фе арргориаже ите. 
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Введение. Данная статья посвящена исследованию способности искусственных нейронных сетей различных типов 
прогнозировать модуль упругости слоев нежестких дорожных одежд в условиях эксплуатации. Для моделирования 
динамических нагрузок дорожного движения и измерения чаш прогибов на поверхности покрытия для получения 
базы данных моделей ИНС была использована установка ударного нагружения Е\/О. 

Материалы и методы. Были разработаны типы искусственных нейронных сетей (сетей обучающихся по 
принципу прямого обратного распространения, послойного рекуррентного распространения, каскадного обратного 
распространения и обратного распространения Элмана) для определения оптимальной модели ИНС с 
использованием программного обеспечения МаНаб. Чтобы оценить эффективность каждой модели были 
использованы разработанные модели ИНС для прогнозирования значений модуля упругости для 25 новых участков 
дорожных одежд, которые не использовались в процессе обучения, проверки или тестирования, чтобы убедиться в 
их пригодности. Для всех результатов моделей были получены такие показатели эффективности, как средняя 
абсолютная ошибка (МАЕ), коэффициент детерминации В”, среднеквадратическая ошибка (ВМ$Е), значения 
абсолютная процентная ошибка (МАРЕ). 

Результаты исследования. Исходя из параметров эффективности, было сделано заключение, что среди этих 
алгоритмов модель распространения прямой связи обладает лучшей производительностью по сравнению с тремя 
другими типами ИНС. Результаты четырех лучших моделей сравнивались друг с другом и с реальными 
полученными данными для определения наиболее подходящего метода. 

Обсуждение и заключения. Различия между результатами четырех лучших моделей для четырех типов 
используемых алгоритмов были очень малы, так как они показали близость между ними и фактическими 
значениями. Результаты исследования подтверждают возможность моделей на основе ИНС быстро и надежно 
оценивать модуль упругости слоев дорожного покрытия для его использования в структурной оценке (МОТ) 
нежестких дорожных одежд. 


Ключевые слова: асфальтовые покрытия, искусственные нейронные сети (ИНС), дефлектометр падающего веса 
(ЕУУО), сеть обратного распространения, неразрушающий тест (МОТ). 
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Шеодисйоп. АП рауетеп: гоа4$ ул д&епогже оуег Ите, гесаг]ез$ оЁ пом \уе| 4ез1епе4 ог сопзиийе4 [1]. 
Реепогайоп оЁР фе рауетеп{ 15 аНецеа Бу мас уоште, сИтае соп41оп, БаПаше даашу, 1ауегз ф1сКпез$ ап ®е 
еРЙсепсу оЁ еагег гевабИканоп ап4 теайтеп р1ап. ОзцаПу, фе рауетепЕ соп@ оп з{ауз 2004 ш Ще Ёгз{ 50-75 % оЁ 
зегутсе Пе ап Фе ргосеззез оГ аепогайоп ргоогез$ $юо\Лу. ТБе деэтадайоп ргосеззез таке гар14 ргосгез$ опсе фе 
рауетепе $а{а$ 1$ десгеазеа [2, 3]. Ргорег таищепапсе ог гераг асну1ез тау $10\/ до\’а ог гезе{ дезотадайоп ргосеззез 1Ё 
9911764 а а зацае Ите. 

Тье змасвага| соп оп оЁ рауетеп сап Бе еуамае4 изше поп-аезгасиуе зигГасе 4еНесйоп {езИп2; пирч[$е 
]оа4 Чеу1сез засВ аз ЕаШие \УееЩ Оейесютаеег (ЕО) апа Неауу \Мез Ш ОПейесютеег (НУ). И 15 Фе то 
Кеадиеп Ну изе4 теазиние шз@атепе юг 15 об]есйуе. Вазе4 оп Фе теазиге4 рауетепЕ гезропзез ш аеНесйоп {е$5, 
шаепа] ]ауег пое$ сап Бе езитае4 изо Баск-са|1сшайоп [4, 5]. 

Тре ЕУ/О ]оа4 1е55 зипи=же аЁИс свагает$Ис$ засВ аз фуре, уоште, ап ите оЁ уешс]е 1оа4те сотгес|у. 
Трезе 4еу1сез арр!у ап пирч[5е 1оаа (Р) @гочеВ а таз$ ш Нее ЁаП оп а сисщаг р|=®е у а суйп4диса! габЪег Биг 
шоицтеа ип4дег фе аШие уе зузет [6, 7]. ТЬе демсе гесог4$ Фе уегаса| рауетепе деогтаноп изше ЧШегеп 
зепзогз |1осаёе4 аё уапоц$ Ч${апсез Нот Ше сепще оЁ Фе 1оадте сие [8,9], аз звоми ш Неиге 1. ТВе тшахилит 


Информатика, вычислительная техника и управление 


365 


БИр://уезиК-доп$ва.га 


366 


А4уапсеа Епятеетпия Кезеагсй 2021. Т. 21, №4. С. 364—375. 155М 2687—1653 





Чзр1асетеп 15 Кпо\’п аз фе реаК 4еЙесйоп, \мсВ оссигз ип4дег Фе оао ро1шЁ. Тга4юопа| те о4$ изе Ше 61е6ез 
уа[мез оГЕ\УР 4еНесНопз {ю Баск-са|сшае Ппеаг е]азис тоди[а$ Гог еасВ |ауег оРрауетеии [10]. 


ЕаШиз мере 
‘ 


р 5еп$0г$ | 





| ! | Фрей = 
1-02, 0) /: 
чи р ще # рауетепе 
о Зы —_ = —_ (еЙесноп Бом 


`` епзог деЙесноп 


Е!о. 1: Рауетет АеНесйоп ассог4те №0 Е\/О 1езйпе 


Випис Фе ра${ Ге\ уеагз, 11 Ваз Бееп об5егуе4 па{ Ше рауетеп{ адтиизгайоп$ изе пе\ апа Ч1Негеп тео4$ оЁ 
соПесйпе ап ргосеззше дайа юг гоа4 тапиепапсе [11]. ие ю Фе гар1А деуеюортепе о? шЮюгтайоп {есбпоозу ап 
агяНс1а! п\еШоепсе, Феге Бауе Бееп шире орробишиез$ юг паретещайоп ш 4еуеорше рауетеп тапасетеи 
зу$е115. 

Моге еНеспуе \уауз ю ад4гезз Фе рго ет о? зрес1 Фуше (сотр[ех, поп-Ппеаг, пииуапае) рагатеег$ звошША Бе 
соп$14Аегед. АММ 1$ опе оЁ Фе агайсла| ищеШоепсе фесбиаие$ ШФаЁ ргоу14е зоайоп$ {ю С1аз;1Исайоп апа геогез1оп 
ргоетл$. [ 1$ Кло\т аз опе оЁ Фе Без фесацез Гог даа пшшо фазК$. Ц Ваз а Нате\угогК Гог Чегет тасфлше |еагитао 
а|>огИ ип №0 регогт 1юхефег ул даа шриб5. АММ Теагиз Бо\и ю ргедс Фе ошриё Нот а $е{ оЁ ай щез. ТВе 
а|>оги тт 1еагпз ф0 Гюгеса${ фигте Фе {гаи ргосе$$, \у№аеВ пли саде даа \иЙ а Лагое дотпали, фо ауо14 а рго ет оЁ 
Раис Фе ехрецме4 ааа ое Фа гапое, УмсВ аЙес{$ Фе уаПАИу оЁ Фе гезиз [12]. ш аа4\юп, ве Недаепсу оЁ 
затрИи>® звошШ4 Бе зиЁйслет ю [еаги соттесЙу. [ Ваз Бееп обзегуе Кот а 10% оЁ тезеагсй а АММ ргоу14е4 соо 
ассигасу ш рауетеп{Е регюптапсе рге@сйоп. Те гоа| оГ АММ 1$ юЮ Йп4 зо[айоп$ фо ргоет$ ш а зппПаг таппег аб а 
Батап Бгаш 4оез [13]. 

ш 2004, Фе аифог$ забиие а Югишайоп Кг Фе геуегзе са|си]айоп ог ®е рауетет 1ауег то@щез изше 
агаНсла! пеига| пебхогк$ (АММ). ТВе гезеагсВ Баз зВо\уй ай Фе ргорозед АММ шефо4 пее4$ сопз1Аега у 1е5$ те 
сотрие4 ап офег тео4$ зисВ аз 1ауегеа еазИс Феогу, едшуаеиЕ [ауег 1сКпез$ (ЕГТ), апа Нице-еетепе тефод$, 
гезресиуе]у. ТВе АММ 1$ изе4 ш зпищайоп ав а 61еВ гае Бесалзе ( сап |еаги сотр ех попПпеаг ге]айоп$ мрз [14]. 

НаШ Сеуапе{ ап фе сопиТбийпе аифогз (2008) деуеоре4 Ъаск-са!си]айоп то4е]з$ Базе оп агайс1а| пеига] 
пебхогК$ (АММ№ ) Юг рге@сипе Ше еазИс тодиа[а$ ог Ше РогЧап4 сетепе сопсгае (РСС) Тауег ап Фе соеЁИсепЕ о? 
забога4е геасйоп Юг Фе рауетепе Ючипдайоп. АММ-Базе по4е]$ Вауе Бееп гашеД 10 езитае Фе 1ауег тодие; ми 
деЙесиуе Базш жа (ЕУУО) апа Фе рауетепЕ збгастаге @1сКпез$. Те гезеагсВ шЧсаез Фа фе АММ то4е|5 сап ргед1с{ 
Бе 11214 тодще оРГраутшс 1ауег$ \И №11 ргес1$1оп [15]. 

Моге зе \еге сопдице оп Фе деуеюортепе оЁ тоге ассигайе ап еНесиуе то4е!5 уу а[еогИиа$ ог 
орипитайоп ап4 Вуб4 зу%ета$ [16-19]. 

О. Вейгап ап Фе сопилриНие аафогз$ (2014) соПемед даа Нот Не 1е3{5 ю тесайси]=ме 1ауег то4е]5 Ъу 
агаНс1а] пеига| пебмотК по4е!5. ТБе гези5 ргоуе4 фе еЁсепсу ог Фе АММ то4е]5 ш сасайпе Фе рауетет5 
рагатаегз [20]. 

Маоуици М ап4 Нао У’апе Би а тоде! Юг са1си]айие Фе е]азисйу тоди$ юг Нехе рауетепЕ 1ауегз из 
Бо Фе сепейс а1оогт ап Фе агийсла| пеига1 пеб\уогк зубет Базе оп фе Шпе уеюШ деЙесютаег даа. ТБе 
гези 5 оР Фе АММ-СА тоде! зВо\е4 геазопае ассигасу ут Фе даа геол$ете4 ш Фе ГТРР 1е3{ ажабазе, \Веге еге 
у\еге по Ме Ч1Шегепсез Бебмееп Фе ргедщеа уашез оЁ Фе е|азйс шодиаз Гог Фе азрвай зи! {асе 1ауег ап Ше геол${еге4 
ш Фе ГТРР даа [21]. 

[1 1$ шуезисайоп, зеуега| апа[узез \уеге регГогте4 10 аейпе Фе Без розз1е агспиесвиге а1опх м Пеагиао 
пез ап Фе фуре оГ Ме АММ тоде! ю шсгеазе Ше Югесазипе сара Пйез ог АМ№. ТВе изе4 дэжаБазе шс[адез ул4е 
гапоез ог деЙесйоп уа[аез оБМашеЯ Нот ппрасй 1оа4 {е55 сопаисе4 оп ех15Ипе гее-Шауег рауетепе зузет$ оп Фе 
гоа4$ пебмогК Бу Фе Зае Сошрапу Ви5$1ап Н1еб\’ауз$ Нош 2014 ю 2018. Ц \аз ичй|7еа аз ап ехрегитена| Баз1$ г 
таштие агиЙсла| пеига] пебмотК тоде[5. 
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Маена[5 апд Мео4$ 

ТВе ргорегИе$ о? Ве изе4 зесйоп$. 

ш ТаЫе 1, Фе ассера Ме Шшии5$ оЁ фе рауетеп{ |ауег рагатеет$ изе4 ш Ба|Й4ше то4е| 5 аге тепнопе4 Юг 
сайси]айпе Фе е]азисйу тшодиаз Юг Ще рауетеп: 1ауегз (загЁасе, Базе, ап4 забогаде). 


ТаЫе 1 
[1115 ОР сеотейтез апа ргорегиез оР таёета15 Гог рауетеп зесйоп$ 


Маепа| буре Гауег 1сКпез$ (пл) Ро155оп’5 Кано Гауег е]азис тоЧи[а$ (МРа) 
Азрвай Сопсгае Т(АС) 190 : 220 У = 0.35 Е, = 900 : 4500 

Вазе-[ауег Т(В)= 350 : 460 У = 0.35 Е,= 80 : 500 

зибота4е-|ауег Т($)= У = 0.35 Е. = 40: 150 


Васк-са1саноп то4е]$ Базе оп АММ 5 арргоасв. 

[ш {1$ могК, \е изе4 Фе БасК-ргорагайоп а[=отИт Рапсйоп №0 зоуе Фе рго ет о{ Ше попПпеаг РлисНоп 
таррше, уПеге и Ваз шов еЁсепсу Бебуееп АММ а[еогииа$ [22—25]. Еафегтоге, АММ пебмогК$ оР 1$ фуре аге 
дейпе аз Фе пеига! пебмуогК$ оГ шиШауег Ееед-Ююг\аг4. ТВе шаЧ1опа| агсИиесвиге оЁ агАйсла| пеига! пебмогК$ 15$ 
ргезегуе4 ш 1$ аеогИитл. ТБе эасфаге оЁ 1$ а[еогИта со1$1${$ оЁ шри ап ошри6 гергезещеа Бу пеигоп$, ап 
Бебмееп Фет, Феге аге соппесйопз$ изе4 №0 НапзРег е \е1о{5 олуеп ю еасб се ассог4ше ф0 5 еНест. ТВе БаскК- 
ргоразайие аеогИ фт 1$ сБагацетте Бу из ау ю сВапее е пеигоп \е1ю5 {0 тедисе фе 41Шегепсез Бебмееп Фе 
с0а1|$ ап Фе опий уа[ез оГ Фе а!еогт изше Фе еггог гедасйоп {1есбиламе [26]. ТВе Нпа| зе{ оЁ поде Мазез апа 
соппесйоп \е12 5 1$ Кпоууй УВеп Ше еггог гае 15 тедисе4 1ю регил1551е ши [27]. 

Тре пебмогК 15 ташед Бу Чегет а[оогИа$ УИ Фе гаш дабазе. 

— Еее4д-Еог\уага Моде] 

Тре Еее4-ГРюг\ууага пебмогКк шуо]уез а 1еа$ гее 1ауегз (фе шриё 1ауег, @е м14еп 1ауег, апа Ше ошрий Тауег), апа 
( ау шсгеазе фо Вауе тоге ап опе №144еп 1ауег ассог4те фю Фе пебмогК пее4$. Аз И 15 Чеаг Кот Фе паше оЁ Фе 
пебхотК, фе шЮгтайоп Ваз опе Апесйоп Нот Ше 1приё ю Фе Ша4еп 1ауег ап Феп о Фе ори! Тауег, аз ш Е1отге 2. 

— Гауег Веситет Моде] 

ТБе угастаге$ оЁ Ше гесиггег( пепга| пебмогК ап Ше Еее4-РЮг\агА пе\огК аге зипПаг, Бай Фе гесиггеп( пепга1 
пебмотК 15 итаче ш фай Феге 1$ а зреслИс ЕееФаскК 1оор ю еасВ 1ауег ш Фе пебмо!К ехсер( Юг Ше 1а31 Тауег, аз зВо\уй ш 
Е1еиге 3. ТЬ1$ ГееаБаскК 1оор реги е пебхогК {о Вауе ап ипитие4 упатс геасйоп {ю шсотше ите зепез даа. 

— Сазсаде Еог\уата Мебмогк Моде] 

Тре Сазсаде Гог\уагА пебмог!К 1$ зийПаг 0 а Геед-Юг\’агА пебмотК, уп Фе ошу ЧШегепсе Бешх фа и шсааез а 
ПоК Кот Фе шриё фо еасВ |ауег ап Нот еасВ |ауег фо Фе ЮПоулизх 1ауегз. Аз звоми ш Е1елге 4, Фе Сазсаде Тог\уаг4 
по4е! ргодисез ПиК$ Нот Фе г {№0 Ще зесопа Гауег, Кош Гауег 2 ю Гауег 3, апд Нот Фе Игз ю Фе Шиа Гауег. 
ТБезе пеё\могК$ а150 ргом1Ае шриё соппесйопз$ Юг аП Тауег$, упеге а4Чопа| ПпК$ сап аилсКу пиргоуе Ше Теагип? 
ргосез$ оЁ Ше пебмогК то]. 

— Ейштап Меига! Мебхмогк 

Тве Ейтпап пеига| пебмогК згасваге оуегАез Фе Геед-ГЮюг\ууага пебмогКк Бу Ваушс а 1ауег саПе4 Ше сощехЕ 1ауег 
ш ааа оп {0 Фе шриь ошри апа ыА4еп 1ауегз. Тбе лисйоп ое сощех( Тауег 1$ ю юге Фе ошриЕ оЁ Фе ша4еп 1ауег 
ш еасВ Нит5Беа сусе ап4 тецзе и аз при © Ше шадеп 1ауег ш Фе пехё Кегайоп © епзиге фай равегиз аге сепегае4 оуег 
ите, аз звома ш Е1еиге 5. Ейтап пебмогК$ а[5о тедисе Ше еггог гаёе ш Фе ори ю Ше реги 1Ые шииз Бу изо Фе 
Баск-ргорасайоп Геайте, аз 15 Ше сазе ш фе Ююг\уга!г4 ргоразайоп пебмогК. 


Х. Меига! Меблогк Тгайтила (ппёгаитоо!) — | >. $ Мета! Мещмокк Теаиупд (ппёгаииоо 


ы Меига! МеблогЕ 
Меига! МеблогЕ 


ОСифрыт Гауег 


г иииитииоииииинтниинииииЕтиннниий, 
м Ошри 























Е!с. 2: 5гасйие оЁТеед-Юг\аг пебмотК Е!о. 3: уасйише оЕ]ауег гесиггеп( пебмогК 
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-&. Меига! Мевмогк Ткайтилд (пийгайтос!) — О Хх. &. Меига! Меблогк Тгайтила (ппёгайтоо!) — | х 





Меига! МейлюогК Меига! Мевлгогк 


Е!с. 4: Збгасвге оЁ сазсаде пебмогК Е!о. 5: угасйше оЕЕйпап БасК-ргорагайоп пебхогК 


Тргочев фе соплтапа оЁ “пиюоГ” ш Фе МаЧНаБ ргоэтат, ме паретеще4 Ше сгемеа АММ то4е|з. ТБе ааа фа 
уаз изе4 ш Ба|Й4те Ше п1о4е]$ Беопх {о а этопр оЁ азрВай зесйоп$ ор Фе М-4 шеб\ау оЁ Фе Визз1ап гоа4 пебмогК. ТБе 
5е{ оГ ташише жа изе4 ш 15$ за4у шс[аде4 е гези5 оЁ сай1сиайп» фе еазис тодий о? е згасага| |ауег$ оф поп- 
[1214 гоа рауетепт, саге ой ш а зресла1те4 зоЙ\аге расКаге зиррПед \ий ап Е\УО Риштах Воск 1юадте ип оп 
555 рауетепё зесйоп$, тапаге4 Бу Фе 5 Сотрапу Кизз1ап Н1еб\ауз. \Ме изе4 Юютг фурез оЁ агайсла|! пепга| 
пебуогКз, урасВ \еге (ее4-Гюг\мага, 1ауег гесиггеп®, сазсаде, ап Ейтап) Баск-ргораганоп ю се{ Фе Без{ гези $. Зеуега1 
сеП$ \уеге з@есе4 ш е В144еп Тауег Юг еасВ фуре оЁ тоде| ю зи4у 5 еНесй оп Фе Наше ргосез$. Тве ргостат 
Чту14е4 Ше даа а гапдотп УЕ 70 % Юг Фе гашиш ргосез$, 15 % Юг уемйсайоп, апа 15 % Юг Те {е5йпс ргосез$. Е1еиге 
6 5По\уз Ше агсбйесвиге оЁ Ве агиЙсла| пеига1 пебмотК тоде]. 

\Те сагле4 оцё @гее %ерз ю оМаш Фе орита| питбег о{ пеигоп$ апа №144еп 1ауегз ш АММ то4е[5: ме ташеа 
Фе то4е|! у ЧШегепе Б144еп 1ауегз (Нг5 %ер), езитае Фе гези5 оГ Ше регогтапсе то4е| (зесопа $ер), апа 
сотраге4 Ше ргефсе4 уаез о? Пе {ее даа ю Ще 1агое{ уаез$ (та %ер). 


три даа 






мае ьные = 
АХ 2 ОСНО 
ххххххххх) 


Е!с. 6: Меига| пебуотК агспиесиге 0 даегилше Фе е]азис поди оЁ Фе згисга]1 [ауегз оР Фе Нех!е рауетей5 


О1, 02, 03 .... 010 аге Ще геза 5 оЁ Фе теазигед деЙесйоп уашез ип4ег 1е зепзог$ — сеорвопез; Р 1$ Ше 
зиг{асе {етрегаиге оЁ Фе рауетеп; Р 1$ Ше ргеззиге оп Ше рауетепЕ; Т1 апа Т2 аге Фе ЧсКпез$ оР Ве [ауегз о азрвай 
сопсгее апа Фе Ч1сКпе$$ ое Базе |1ауег оР Фе рауетепЕ; ЕТ, Е2, ап4 ЕЗ — Фе е[азйс то4иП оЁ азрвай сопсгае 1ауегс, 
Фе Базе, апа Фе заБбога4е, тезресиуеу. 

Тре Мо4е!5 Еуалайоп. 








= Тре сотгеспез$ оЁ Ще уашез оЁ Фе ргед1сйоп геза5 Юг еуегу то] 15 са]сшаеа изше Фе теап абзопжще еггог 

Е (МАЕ), Ше шире деегиитанопз соеЁйсеги В”, Фе шеап гоо{ зацаге еггог (КМ$Е), ап4 Фе теап абзопие регсетасе 

я еггог (МАРЕ), \шсьЬ аге д&егите4 Нот Фе ЮПо\лие Югич1аз (1, 2 ап 3): 

э 1 

ю МАЕ= - т [А,- Е (1) 

я 

— т 2 

о ВМ«Е= рее (2) 
п 

р МАВЕ= - Ут 45 *100 3 

Е а - уе А, (3) 

-5 
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М/Пете А, 1$ Ше асша] уаше ш репо 1; Е, 15 Фе ехреце4 уае ш репоа 1; ап п 1$ Фе питфег оЁ Фе 1ю{| репоч. 
МГ ИВ гезресе ю Фе %айзИса! ш1сез, МАЕ, ВМЪЕ апа МАРЕ, зтаПег уашез изиаПу ш@сае Б1оТег ассигасу гези 5. ш 
{15 апа[у15, Ше МАЕ, В^, ВМЗЕ ап4 МАРЕ уаюез Юг еуегу тоде! аге омате4 @гоиэВ сотрагше 1е ргед1се4 гези 
аса1п${ бе асфла| уа1аез. 

Вези! 5. Еоиг буре оЕ ММА \еге деуеореа \ий Ююиг Чегет пеигоп пипбегз ш е В144еп 1ауег$ фо зее \мсЬ 
опе 15 тоге зийае ю изе ш Фе Ююгесазипс ргосез$. ТВе {а питфбег оГ АММ тоде|[$ сепегае4 \аз з1жееп. АП то4е15 
уеге Наше ип4ег уатой$ сопЧ1юоп$, шсаше 10, 15, 17, апа 20 пепгоп$ ш Фе Б144еп 1ауегз. ТЬе уааез оЁ Фе 
ЮюПо\ше МАЕ, В^, ВМ$Е апа МАРЕ та1сез ог а| пло4е!$ \еге са1сиае4 аз зпо\п ш ТаЫе 2 ю аззезз фе то4е! 
регогтапсе. 


ТаЫе 2 
Еуамайоп оР фе то4е|$ регогтапсе 


ргефсноп Е(АЗ) оЕРеед- ргефсйоп Е(АЗ) ог ргеФсйоп Е(АЗ) ог ргеФсйоп Е(АЗ) о Таует- 
Гогу’аг4 то] сазсаде то4е]5 Е]тап то4е|$ гесиггеп то4е|5 


пт 1590 |100 [095 [109 


МА 
аа аа а а а 


ргтефсНоп Е(Базе) оЁ ргефсйоп Е(Базе) оЁ ргтефсйоп Е(Базе) оЁ ргтефсйоп Е(Базе) оЁ 


Гее4-Гог\уаг4 то4е]$ сазса4е то4е[5 Е тап по4е[5 [ауег-геситеп то4е|5 


| 1.00 | | 1.00 | 1.00 | 1.00 | 1.00 | | 1.00 | | 1.00 | 0.99 | 1.00 | | 1.00 | Ро ЕИ-ТОЙ | 1.00 | | 1.00 | | 1.00 | 


МА 


ргефсноп Е(5аБ-ога4е) | ргефсйоп Е(5аБ-огаде) оР | ргедсйоп Е(заб-ота4е) ог | ргефсйоп Е(5аБ-отаде) оЁ 
ог Гее4-ГРюг\уаг4 по4е|$ сазсаде по4е]$ Ейлап п0о4е15 ауег-гесиггеп" п104е]5 


пар | 0.99 | 0.95 | 0.87 | 0.99 | 0.93. | 1.00 | 0.97. 
МА 


Тре Без{ Юг п1о4е[5 {Ваё ехргез$ фе Юг {урез оЁ агайсла| пеига| пебмотК$ \еге з@еце4 Базе оп Ше уашез оЁ 
Фе апа1уйса| рагтатаег$ тепйопей ш е ргеутоц$ {аЫе Юг 1е гези $ ое то4е[5 © сотраге ет апа Кпо\ Ше ежеп 
ог ег ппрасе. \УТегеаз, Фе Без гези{$ \еге дие ю Фе то4е|$ уп вхо Ша4еп 1ауегз апа 14-15-3 эгасваге юг аП 





еуеореа 1урез. АП то4е[ 5 оок е тахипит пипабег оЁ гереййоп$ ю сошр!ае Ше ташише ргосез$, \меВ \ете 1000 
гереййоп$. Ви Феу АШете4 ш Фе репо4$ 1аКеп ог Нашишо, аз Чаше ог Фе то4е|$ епае4 аЁег (54, 157, 67, апа 173) 
зесопа$ Гог фе Юг по4е|5, гезресиуеу, аз зво\и ш Е1еигез (6-9). Веше а\уаге оЁ баф Фе нашио ргосез$ юр Е Ше 
тахитит пипбег оГтерейНоп$ ог Фе Ите гае 1$ ехсеедеч. 
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$ Меига! Меблогк Тгаитипа (ппёгаи оо!) — 9 х 


Мега! МелгогК 


НЫ еп 1зуег Оцтри Г ауег 





АюогИ тт 
Оаеа Сиугяог: Вапдогп (Оп егапа) 
Теа о: ЕемепЬегд-Магдцаг Е (апп) 


Реоггпапсе: Меап 5диагед Ёггог (пе) 
Сасшаноп::  МАТЬАВ 


Ргодге$5 

росы о оенвок 10% 
Тиле: 
Ре{оггпапсе: 1.21е+06 м 0.00 
бтафепЕ 1.90 +06 [ООО ООООООООООООИИИИИИ | 1.00е-07 
Ми: 0.00100 1.00е+10 


Уайдабоп Спескз о ю 


Е1с. 6: Тгашше уидо\у оЕ Ееед-ГРог\уагА то4е| 


«$ Меига! Мевмогк Ттаититд (ппегатоо[) — О х 





Меига! МенлгогК 





АщогИР т 
Оаа Отуехюют: АапФогл (ЧпиЧегато) 
Тгаипит о: ГеуепЬегд-Магоцаго Е  (+гап!ип) 


РеГогплапсе: Меап 5дуцагеЯ Егог (ип5е) 
СасцаНоп$: МАТАВ 


Ргодгез5 
Еросн: о 100 | 100 
= О: 0% 


Ре{оггпатсе: 5.7 [ево | 0.00 
тает злТе+0з [ПО | 1.006-07 
Ми: 0.00100 1.00е+10 
МабЧаНоп Сеск: о О | 1000 


Е1с. 8: Тгашше уп доу ОЕ сазсаде то4е| 


4 Меуга! Мебмогк Тгатпито (пига оо) — х 


Меига! МевлогК 





Аююог т 
ОСаЁа Отуяюг: Капдогп (9пл9егап9) 
гаи: Гемепьегд-Магоиаг* (+гашпИип) 


РегРоггпапсе Меап %диагед Ёггог (ип5е) 
СасцщаНнопя: МЕХ 


Ргодге55 

Ерось: о 1000 енот] 100 
Пите Г 0 |] 
Регоггпапсе: 28.1 [иена | 0.00 
бгадтеп б.17е+03 [ПО б ии — | 1.00е-07 
Ми: 0.00100 1.00е+10 


Уайданоп Спескя 10 


Е!о. 7: Тташше ула4до\ оЕТауег гесиггеп{ то4е! 


«$ Меига! МебмогК Тгантито (ппбгайтоо!) — С х 


Меига! МеблогК 





Ащог тя 

ата Оллятог: Вапдот (Филд егапо) 

гало: [емепЬего-Магоиага* (гаи) 
РеПогтапсе: Меап 5дуагед ЁЕггог (пп5е) 
Сасшанопе: МЕХ 





Ргодгез5 

ро о [_— жене | пох 
ое 00| 
Ретоптпапсе: 6.79 [ПООиоиази |] 0.00 
бгадтепе 1.07е-- 03 [ ОООООООООООООООСООСООООИИИИ — —_| 1.00е-07 
Ми: 0.00100 1.00е+10 


Увйданоп Спескя м] 10 


Ес. 9: Тташае ула4до\ оЕЕПпап БасКргорагайоп 


Тре ууе1е 5 апа Мазез оГа рагатаег$ ууеге то 1Йе4 © Десгеазе Фе еггог Бебхееп {агоей уа[ез ап пебмогК 
ошриЕ ШгоиеВой Ше таше рвазе. ЕасВ пеигоп \’е1е( 15 аара%е Базе оп {$ парасё оп Фе пебмогК гези. ТБезе 
\е1215 апа Мазез \еге еуаае4 дигие атшие аз Фе \е12 5 оЁ Фе при тю Ше шаЧдеп 1Таует. 

Е1сиге$ (10-13) Шазгае Фе теап зацаге4 еггог уашез уз Фе питфег оЁ цеганоп$ Гог Ше Наштае ргосез$ оЁ е 
сотраге4 по4е[5, гезреспуу. У\УПеп изше Фе Еее4-Рюг\аг4 Баск-ргорагайоп то4е]|, фе Бе уайЧаНоп регюогтапсе \аз 
16.68 а еросб 1000. ОФегулзе, фе Ъезё уайданоп регогтапсе Юг Фе Тауег гесиггепь сазсаде апа Ейпап Баск- 
ргоразайоп по4е|5 \меге (11.52, 119,888 ап4 28,067) аё еросВ тего, гезресиуе[у. ш аП оЁ Фе сигуез$, ме обзегуе ве 
сопуегоепсе о? е {е5{ сигуез УИ Ше уа|4айоп сигуез, \шсь геуеа[е4 Ба{ Ше {е54 ап4 уаПдайоп сигуе$ аге уегу зипПаг. 
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2 Везх Уа1дайоп Ре {огтапсе 15 16.6883 ай еросВ 1000 2 Ве УаНааноп Регоппапсе 15 11.5246 а еросВ 0 

н Тгай н 

ы Е 

= = | 

о о 

Е Е 

5- 5- 

2 2 

& 10° Я 

Ф Ф 

> 0 200 400 600 800 1000 > | —ы 

0 200 400 600 800 1000 
1000 ЕросВ$ 1000 Еросвз 
Е!с. 10: Реюгтапсе оЁТее4д-РЮюг\уаг4 то4е1 Е!о. 11: Реюгтапсе о [ауег геситгеп по] 

5 ВезЕ Уа|Чайоп Регогтапсе 15 119.8888 аё еросй 0 Ве Уа|Чайоп Рейогтапсе 15 28.0672 аё еросв 0 

= Тгат $ 2 
— \УаЙда оп 1 0 Тгат 

о Тез! ви МаЙдаНо 
Е Вез1 Е Тез 
М>) 2] Вез! 

о М> 

Е 10°] 

3100 

> 5. 10 

= 2 

о = 
2 о 

0 200 400 600 800 1000 > 0 200 400 600 800 1000 
1000 ЕросВ$ 1000 ЕросВ$ 
Е1с. 12: Ре!оптапсе о{ сазсаде то4е] ре{югтапсе Е1е. 13: Ре!огтапсе о ЕПпап тоде! 


Ассогаштс 1ю Е1оигез (14 ап4 15), агайсла] пеига| пеб\могК$ аге еуаааже4 изтс а сот@айоп Рапсйоп (В) Бебмееп 
пебуогК геза уа[аез аза115 {агое{ даа ш ее 4-РЮюг\аг4 ап4 сазса4 те Баск-ргорасайоп то4е|5. ГооК1ие ай фе К-уаме оЁ 
фе по4е[5 дите Фе таште, уа|Чайоп ап4 {езНпо рВазез, ме пойсе а соо4 аэтеетепе Бебмееп Ше фагое{ уа[иез ап е 
пебуотК тези 5. У Пеге \ме Нп4 Ша Фе $оре оГ Ше Ппе Фа+{ гергезеп{$ Фе геаанопз$ р Бебмееп ааа ап гези$ 1$ 
арргохипау еама!| о |1 ш а %асез оЁ гашитте, ап 1$ шд1саез Ше гоби$тез$ о? Ше то4е[5 апа {Ве роз ИЩУ оЁР из 
Фет Фигте Фе ргеФсйоп ргосез$. 








Тгашше: В=0.99999 \Уа|дайоп: В=0.99999 Тгашшо: В=0.99997 \Уа|даноп: В=0.99996 
м 4000 — = < 4000 
_ _ Е + 
5 3000 5 ых % 3000 о 
ог г з з 5 
з 5 Е- Е 
== 2000 Е. 3 > 2000 я 
С — Г ] 
| 1000 р - Е 1000 Е 
Е Е Е 2 = 
Е = Е - © 
5 5 5 5 З 
1000 2000 3000 4000 1000 2000 3000 4000 2000 — 4000 о ы 
Тагое Тагое Тагое Тагое _ 
Тез: В=0.99993 АП: В=0.99998 Тез: В=0.99999 Тез: В=0.59997 = 
5 Я 
= 4000 $ 4000 2 4000 5 = 
+ - 5 3000 Е с 
В 3000 200 5 Я а 
: ы Е" 2000 - 5 
= 2000 == 2000 =. — т 
] — Г = 
| ] а 
5 1000 о ее” Е о 
а = 8 = з 
Е В 5 = => 
©) Е о © = 
1000 2000 3000 4000 1000 2000 3000 4000 2000 — 4000 . ыы о 
Татое Тагое! Тагое агое Е 
м 
и. 
[2 
Е!с. 14: Ееед-ГРюг\ууаг4 ое] геотезз1оп сва[$ Е!о. 15: Сазсаде БасКргорагайоп геотез1оп спаг$ а. 
к 
АЕ Фе Беошиис ог пебмогК [еагита®, аЙ ]еагилио $асез оЁ Наше, уа|даноп, апа {езйпе, гай ш ап ореп 1юор. : 
Опсе Фе ташитие рВазе 15 оуег, фе пебхогК регРогтапсе фагп$ 1 а с10зе4 |оор дигие Фе ши-5%асез оЁ ргед1сйоп. [9 
о 
-Е- 
и 
ых 
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Тре плоде] $ ба соггесЙу гергезет{ фе ауа|Пае даа аге з@есеа дигио фе уасез оЁР едисайоп, уеиНсайоп, ап 

{е5йпе, Базе оп Ше зай$Иса| рагатщетз, аз шепнопе ргемоиз[у. Ц \аз геаишед ю даегтше Ш Фе то4е]5 соч19 

Гогеса${ бе е]азис тои[а$ юг 25 зеотеп Ша Ва по{ Бееп чй|те4 дигти> то4е фгалите ап сотраге ет {о 1агое 

уа[ез, аюпо УИ Кпоммие Во\ шисВ Фе пебхогк буре аНеме из регЮюгтапсе ш {1епиз оЁ @азис тшодиа$ 

аззеззтепе. ТаШе 3 зВо\из Ше сотрат15оп ог Фе е|азис тшоди[а$ фагое{ уа[аез юг аП Тауегз (азрБай сопстще, Базе, апа 

зибога4е) оЁР 25 зесноп$ ап е Югесазипс гези {$ уааез оЁ Ше юг АММ (бурез (Ее4-ГЮюг\ага, 1ауег гесиггепф, сазсаде, 
апа Ейтап) Баск-ргоразайоп то4е[5. 

ТаЫе 3 

Тре сотрагзоп оЁ Фе ргед4сНоп$ ое ргорозе4 то4е[$ 


Ргефлае4 уа[аез о? Геед- Ргефлае4 уа[аез о? ]ауег Ргефлае4 уа[чез$ оЁ Ргефлае4 уа[аез оЁ 


Асва| тоди[а$ оЁ 


е1азисиу, МРа Гогууаг4 тоде|, Мра гесиггепе по4е|, Мра сазсаде БасКргор ЕПлап БасКргор то4е1, 


(151) (151) тоде|1, Мра (15п) Мра (15п) 





ВиЬ- 
ога4де 


иБ- иБ- иб- иБ- 


Гы ов [пет зто [лев 189 485 [ева [1877 эта | 18833 199 92 


р а 


А Б128 дестее о{ сопуегоепсе сап Бе зво\уп Беёбмееп Фе Готг пто4е]$ уВеп Феу аге сотраге4 {10 еасВ о ег ап4 
{о аса| даа, аз зВо\’и ш Е1еиге 16 (а, Ъ, с). [ 1$ а[5о ипдегюоо4 Нот Фе геза 5 фай Фе по4е]5 сап гедисе Фе еггог 
дигио Ще ргосе$$ о ргед1сипе Фе е[азИс то4ди1аз, уупеге Ше 1о\уе5{ уаше оЁ Ще сое слет оЁ шаре даегилтаноп 
Уаз “В? = 0.95” юг Ше сазса4е ФасК-ргораганоп то4е!. Но\еуег, Ше Ъезё регогтапсе ш Ше рге@сНоп ргосезз, аз 
звоуи ш Та Ше 2, \уаз Юг Ше Ееед-Ююг\иаг тоде| ут а 14-15-3 зтгасеаге, ушсЬ Вад Фе 1ю\е$ уамез Гог Фе зайзиса| 


БИр://уезиК-доп$ва.га 


РИ 
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рагатлеегз МАЕ, ВМЗЕ, апа МАРЕ, ап4 е 1агоез уаае ог 1е соеЁйсепе оЁ пиаре деегитаноп (В^) юг Ше Шгее 
рауетеп 1ауегз (азрБа®, Базе, апа заб згае). 





3500 400 
3000 350 
2500 3 - 
= 250 
Я 2000 = 
— З 200 
5 1500 = 6 
5 
1000 
100 
ы , ИИ 
| | ИИ 
35 7 9 [1 13 15 17 19 21 23 25 
Бесноп № Бесноп № 
м Асша! поди оРе!азисиу м Еее4-Гогугага плоде] (151) м Асва| тофи[л$ оЁейазнсйу ш Еее4-Гог\угаг4 тоЧе! (15п) 
м [ауег гесиггепе то (151) м сазсаде то4е! (151) м [ауег гесиггепе то4е! (15п) ш Сазсаде под] (15п) 
м Еплап то4е] (151) и Е|тап тоде! (1510) 
а) Ь) 


Е(зиботаде), пра 


о ИИ 
» ПЕНИИ 
‚ АНЕНЕЕЕНИИИИИИИЮЮЮЮЮНИЕ 





Зесйоп № 


ш Асса! тода[а$ о @азисиу ш Кее4-Гог\уагА то4е!(15п) 
м |ауег гесиггеп тоде] (15п) м Сазсаде по! (151) 
ш Е]тпаи и104е1(15п) 

с) 


Ео. 16 (а, Б, с): Шазганой о Фе сотрагзоп Бебуееп Ше уа[аез ое 1агое( е]азисиу тоди[а$ ахаш$Е 1е ехреще4 уашез оЁ Ше 
еуе!оред АММ тоде[$ ап еасВ о ег Гог а рауетеп: [ауег$ 


015си$$10п апа Сопс$1оп. ТБе аиог$ сопзгасцеа з1жееп по4е15 ш Ше зшау. Еуегу Ююиг то4е15 КЮПо\и Ючг 
Чегет агиЙсла1 пеига! пебхогК {урез (Ееед-КЮюг\ата, 1ауег гесиггет, сазса4е, ап Ейтап БаскК-ргоразаноп) ю сасиае 
Бе е!азис шоди[а$ оЁ Фе рауетеп: 1ауег$ зиБ]ес{е4 ю дупапис 1юа4д. ТВе Е\УШ {е31 уаз зеце4 о гергезеп{ ®е упатс 
|оа4 сепегжеа Бу Ще гоа4 ‘гаЁЯс ап ю теазиге е геасйоп оГ Фе рауетеп{. МаЧаб зой\аге \уаз изе4 ю стеме АММ 
п104е15 чи71ие деЙесноп аа Нот Фе (М4) еБ\ау ш фе Киззлап Еедегайоп Коа4 пебмотгК. 

Тре а1егепсез Бебмееп Ше гези {$ о? Ве тг Бе5{ по4е]$ Гог Фе Юю\г фурез оЁ зе а1еогИл$ \еге уегу зтай, 
аз Шеу зВо\е4 с1озепез$ Бебмееп ет ап4 Ше {агоей уаез. ТВе Без гези{$ \еге юг Фе Гее4-Ююг\аг4 тоде у 15 
пеигоп$ ш фе \о-В144еп 1ауегз №0 юг Ше эгасйше оГ Фе тоде 14-15-3, \еВ ргодисе4 ®е Ъез{ уашез Юг Фе 
анзИса1 соеРАслепк. 

ТБеге 15$ по сот@айноп Беб\ееп ап 1псгеазе ог а дестеазе о{Р пеигоп$ ш фе 44еп [ауегз ап ап пиргоуетеп т 
п1о4е15’ регогтапсе. шз{еаа, Ее Деслз1оп 1$ 4ерепдеп оп ила ап еггог. 

Егопа Фе геза5, уе Ни Фе сарабиу ор АММ то4е|$ © ргед1с фе еазис тои ог НехЫе рауетеп 
сотгесЙу ап4 зе # ш тапасше Фе рауетеп д&епогаНоп$. 
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АскпоУ\едотеп5. КезеагсВег МоватеЯ Моза ЕП 5Вату 15 Рап4деа Бу а зсво|аг$Мр ип4ег Ше ехесийуе 
ргосгат Бебмееп Ше Ага Керчб Ис о? Езур! апа Фе Кизз1ап Еедегайоп. ТВеу Вауе ту ФапК$ ап4 арргеслайоп. 
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Ттттодисйоп. Тга5$ умтасвиге$ аге ул4езргеа@ ш сопзбгасйоп 4ие №0 а патБег оЁ адуащасез, засВ аз есопоту, 
уегзай у, ап зсайа6Ииу. Ассогаше]у, Фет тодеппе ап са1сШайоп аге игоеп{ фазК$ ш Фе дезлеп оЁ Ба|аше 
угисигез. Ащотайс зоайоп {0 Шезе рго епа$ салзез ап шсгеазе ш 4ез1еп еаепсу, сайсайоп ассигасу, ап 
1о\уег с05{5. Тре оБ]есйуе оЁ Ше заду 15$ ю ехатше фе РапсНопаШу ап4 орегайоп а1еогит ог Фе зоЙ\аге тоде 
деуеореа Бу Ше аи ог$ Шай сепегаез Фе осеотекгу оЁ \о-Читепз1опа| $5 згасфагез Гог заб зеацеп( шо4ейпо. 
Маепта апа Мейо4$. ЕоПо\лио Ше гезеагсВ оЁ Фе ул4езргеаЯ 5$ сопйеигайопз, фе с|азз1Исайоп оЁ сВог$ 
ауаЙае ш Фе зой\аге ипдег сопз14егайоп 15$ слуеп. Те тео4 оЁ рагате{ет1тлие а {ги$$ згасаге 15 ргоу14е4. ТЬ1$ 
шефо шса4ез Базе сеотейлс рагатеегз ог Ше збисвте зисб аз 4ипеп$101$, то4е! сопзиасйой гШез, апа 
а4Ч1опа!| Геаагез, аз \е!П аз а сотргецепзуе а|еогийт. ТБе зой\уаге 1$ 4еуе]оре4 ш Лауазсире. 

Кезий5. 'ТВе зоЙ\аге тоде Ваз Бееп ицеогайеа шю а \еБ аррПсайоп Юг сасиайпе \о-дитепз1опа1 го згасвагез. 
То Шазиае ше псйопаЩу ог Фе зоЙ\аге, Фе ехатр|е$ ог изег имегЁсе аге олуеп аз \\еЙ аз ап ехатр!е ргоет. 
ТБе ехатр!е шс4ез сопЯсигайоп ап4 са1сщайоп оГ ап шсПпе4 #ги5$ згасге. ТБе гези $, засй аз заррой геасНоп$ 





ап п\егпа| огсез у ах1а| Юогсе Члаотат, аге ргоу14еч. 

О5сизяоп апа Сопсияюоп5. Озте 1$ зоЙ\аге тодие ут Ше Нате\хуогК оЁ Фе 1001 Юг сасшайпе гоа 
угасфагез аПо\уз Юг Фе зпирИйе4 ргосез$ оР тодейпе ап сася]айпе сотрМех {$$ згасбигез, 4ез1еп ите, апа 
гезоигсе гедисйоп. ТБе зоЙ\аге тодЩе ргоу14ез 1ю015 Юг зрес1Руше уатоц$ фурез о? згасвге$, арр!уше 10а4$ ап 
а$512пте ргорегаез оГа го4 зует, \мсВ таКез Ц а изе | шзиитепЕ юг 4ез1оп епотшеегс. 
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Введение. Ферменные конструкции получили широкое распространение в строительстве благодаря 
экономичности, универсальности и масштабируемости. В связи с этим их моделирование и расчет являются 
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актуальными задачами при проектировании строительных конструкций. Автоматизированное решение этих 
проблем приводит к повышению эффективности проектирования, точности расчетов и снижению затрат. Целью 
данного исследования являлось изучение функциональных возможностей и алгоритма работы разработанного 
авторами программного модуля, формирующего геометрию двумерных ферменных конструкций для 
последующего моделирования. 

Материалы и методы. По результатам исследования распространенных конфигураций ферм дана их условная 
классификация, доступная в рассматриваемой программе. Приведена методика параметризации ферменной 
конструкции. Эта методика включает в себя базовые геометрические параметры конструкции, такие как 
размеры, правила построения модели и дополнительные функции, а также комплексный алгоритм. 
Программное обеспечение разработано на ЛауаЗси1рё. 

Результаты исследования. Программный модуль интегрирован в веб-приложение для расчета двумерных 
стержневых конструкций. Для иллюстрации функциональности программного обеспечения приведены 
примеры пользовательского интерфейса, а также пример задачи. Пример включает конфигурацию и расчет 
наклонной ферменной конструкции. Приведены результаты в виде опорных реакций и внутренних сил, 
построена диаграмма осевых сил. 

Обсуждение и заключения. Использование данного программного модуля в составе инструмента для расчета 
стержневых конструкций позволяет упростить процесс моделирования и расчета сложных ферменных 
конструкций, сократить время проектирования и сократить ресурсы. Программный модуль предоставляет 
инструменты для задания различных типов конструкций, приложения нагрузок и задания свойств стержневой 
системы, что делает его полезным инструментом для инженеров-проектировщиков. 
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Шегодисйоп. Тгаз$ збгисеаге$ аге изе4 ш уапоч$ Ба|аше @етет. Эгасвага! зо[аНоп$ оЁ $5ез аге 
5ииполизред Бу а огеай уапеу, ЧНегто ш Фе 5е{ оЁ рагатаег$ Гог фе ие ап орип1тайоп. ш фе уойА ргасйсе 
оРгизшз згасигез, Феге аге тоге ап 100 уаточ$ зузетз Ч4Иегте ритагПу ш Фе дез1еп о? Фе го4 ище!асе. ТБе 
Баз1с суртса[ зиасвага| еее 15 а Баг а Ваз сописае Таепшо ул офег е]етеп{$ аё Фе подез. ТВезе 4ез1еп$ 
гедлиге сецат арргоасВез фо Фе са1сШайоп ап4 аззеззтепЕ ог фе гейабИиу. ш \огК$ [1-6], арргоасВез © зоте 
сайси]айоп ап4 ехрегилета| тео4$ Гог аззеззше фе гейаб Шу аге шоб >Ще4. Зоте арргоасВез$ ап4 сайсяайоп 
зспетез Гог аззеззше Ше гепабИиу ог дееснуе змисагез аз \меЙ аз зо]уше Фе рго ет оЁ гезюгше Фе ргорегйез 
ое шагу1Аиа| е]етеп{$ аге соуеге ш топозгарй$ [7-8]. 

Ко4 згасфагез аге зубетл$ соп515Ние ог Ппеаг го4$ соппее4 аё подез. Ву Фе пабге оЁ Фе соппесйоп о? 
е]етлеп5 ш Фе по4дез, зу$фетз УИ 11214 по4дез, шисе4 101065 оГ а| @етеп, ап фе сотЫше опез аге 
Ч1зНиеи1$Вед. Егот фе эеотейлс ро оЁ мех, р|апаг ап4 зраНа| го зузета$ аге 41зНие1$ Вед. 

Тга$5 згасге$ аге а зибсаесогу оЁто4 зузептл$ {а геташ сеотеблсаПу ипсвапоеа 11 11214 ] 018 ш по4ез 
аге гер|асе4 ул Бпое4 опез. ТВе 1оа4 1$ аррПе4 №0 #ги5$ згисагез ошу аё Ше по4ез, ап ошу {1епзПе-сотргез1уе 
Гогсе$ аг1зе ш феш @етеп [9, 10]. 

Тгаз$5 згасбигез аге улдезргеа4 аз @етептё$ оР Би 11$$ ап4 згасгез. ТЬ1$ роршагиу 15 4ие © а питбег оЁ 
аЧуатасез [11]. 

Есопоту — о\уег шаепа| сопзитриоп сотрагеа {0 ойег Ба|Пате згасагез. 

Уегзай Шу — Фе зсоре о{ аррПсайоп о{ ‘3$5ез 15 уегу ул4е ап шс[а4ез Биоез, гооР $$ зузе1тв, 
сгапе$, ес. 

Бсааб у — из$ 1еп2 1 сап уагу Кош а е\ шщег$ ю зеуега1 КПотаетз. 

ие ю ШФезе сВагасег$Нс$ оГ гиз5$ез, Фет то4деПпе ап са1сШайоп аге игоеп{ {а5Кз ш Фе 4ез1еп оЁ 
Бато згасвагез. Амютайс зо[айоп © ШФезе ргоеп$ саизез ап шсгеазе ш 4ез1еп еЁЙслепсу ап сасч]айоп 
ассигасу ап |о\уег соз5 [12, 13]. 

Тре оБ]есйуе оЁ Фе заду 1$ ю ехашше Фе РпсНопаШу ап4 орегайоп а[=отИт оЁа зоЁ\аге тодще фа 
сепегаез Ше сеотету оЁ Гмо-4итеп$1опа| {5$ згасбагез Юг зибзедаеп тшо4деПпе. ТЬ1$ ргоогат тодще 1$ а 
сотропеп{ ог Фе ЗАРКВОЗ зегу1се. 

Мажена[5 апа Мево9$. Эшсе 1ги5$ збгасвагез аге рапе] згисагез, 1.е., Шеу соп$15 оЁ папу зесйопз$ (рапе[5) 
1епйса| ог зипПаг ш сеотецу, и 15 роз51Ые о рагатеблсаПу сепеге а то4е] изше Фе геайуе соог4шаез$ оЁ Ше 


Информатика, вычислительная техника и управление 


377 


А4уапсеа Епятеетиз Кезеагсй 2021. Т. 21, №4. С. 376-383. 155М№ 2687—1653 





подез уПеп зрес1Руте Ше зесйоп. ТВе рагатёег$ Юг засВ то4деПпе аге Фе {оа| 4итеп$101$ ог Фе епаге згасаге 
(Бес ап4 5рап оЁ Пе 155) ап4 пе пиибег оЁ зесйоп$ изеа4. 

Тре ГапсНопаШу ог Фе зоЙ\аге тодше ип4дег сопз14егайоп ргоу1Аез Гог Ше то4деПпо оЁ {ги5$ згасфагез оЁ 
6 Гурез оЁсвог4$ ап ел зи Мурез, Череп 1х оп Фе 1асе изеа (Е1олте 1). 


Ее. 1. Згасвага| Фастата$ изед ш фе тодще 


Саззсаноп о{ Ше бурез оР1{ги$5ез изе4 1$ Фе оПохлие [14, 15]: 
1. тесфапоаг \уий рагаПе] сВог4$; 

2. тарето14а[; 

3. длапечаг; 

4. шсПпе4 у рагаПе] свог45; 

5. шсПпеа Нота Бо 514ез (агго\м фуре); 

6. ро]угопа|. 

Сазз1саНоп оЁ Ше бурез оЕ1авсез изеа 1$ Фе ЮПо\улипо: 

А. блапочаг; 

В. папощаг ми ириеб; 

С. длапощаг ул 515реп$1015; 

О. иру’ага Фасопа1| Бгасе; 

Е. до\пугагА Фасопа1| Бгасе; 

Е. м зиб-Н$$е$ [16]. 

Тауаз сре \уаз зе аз а ргоогатилис |апоцасе юг Фе деуеюортепЕ ое шодще. ш 1$ гегага, бе тодие 
ргоу14ез атр!е орропитиез Юг Фе паретещайоп оЁ сои сотрийпе ог Ба|4ше згасеагез ш Фе \меБ-базед 
ииегасе [17]. 

Тре зой\аге тодще 1$ Базе оп фе обБ]ес{ то4е1, \шсв 1$ ап аггау о об]ес{5 4езст1ие еасВ оЁ Фе з1х фурез 
ог $гасбиге$. Гепо@ / ап Бе1>Е Й аге изе4 аз рагатеегз Гог а {урез оЁ згасагез. Рог #га$$ фурез 2, 4, 5, Фе Вес 
Ч1егепсе ай 1$ гедитеа. Рог еасВ {уре оЁ згасйте, ап аггау оЁ оБ]есё5 дезстае Фе эеотефу оЁ 15 за Турез 15$ о1уеп; 
( шса4ез Фе роз1оп оЁ по4дез апа Багз, пез юг сопзбгасиие зесйопз$ (рагбу, шитогие ог згасвага| Вауез, апа 
изасе оЁ зерагафогз). 

Тре аеогИбиа оЁ фе зоЙ\аге тои[е сепегаез, ассог4ше о Фе зресНе4 рагатщегз, ап аггау о? абзо[ие 
соог4шаез оЁ по4ез оЁР а| {$$ го45. Тб1$ аа 1$ заЁРйсепЕ ю 4езсибе Ше то4е!| зеотешу. Еог а №П-Недэеа 
дейшиоп оГа згасваге то4е| уфы Фе тодще, ше ЮПо\злис а4Ч1юопа1| Геаагез аге ргоу14еа: 


ЬИр://уез-допза.га 


— 5ирроп$5; 
— №109е5; 
378 — ое! оЁ Фе 1ей поде гё]айуе © Фе опош о? соотдтаез оп Фе зспете; 
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ЕЕ ЕЕ ЕЕ ЕЕК ЕЕ ЕЕ ЕЕ ЕЕК 

— [оа4$ оп Фе иррег гасага! по4ез; 

— @етепЕ огоирше Юг Фе сопуешепсе о? зибзеаиеп{ ргосез51то. 

Тье под ие орегайоп а|2огифит 1$ фе ЮПоулие: 

1. Сасшаноп оЁ пода| соог4таез$ юг Фе ощег уегаса! Багз оЁ бе зиасваге 1 Шеу ех15ё юг Фе з4ецеа 
угасиге Туре. 

2. СайсШаноп о? Фе эл оЁ ша1у14ща1 зеснопз ап е соог4тае$ оЁ Фе по4ез$ оЁ Ше сешга| заррог Баг оЁ 
ги5$ буре 5. 

3. зедчепна| сайсу]айоп оЁ пода| соог4тайез оЁ еасв зесноп апа Фе Гогтайоп оЁтго4 аггау. 

4. Сасайоп о{ Фе соогдтаез оЁ по4дез юг а441юопа| уегаса1 Багз Бебмееп зесНоп$, 11 песеззагу, Юг е 
з@еце4 зиБбуре ое згас@те. 

5. Ргейплтагу у1зца|тайноп оРе згасфаге Чаотат 1п1$14е бе тодще. 

6. Аз51ептеп{ оЁ фе а4Ч1опа]| ргорег@ез. 

7. аа тапзЁег Юг дез1еп ап4 са|сшаяйоп. 

АЁег Ше даа Баз Бееп гапзРегге4, Ше аггау Ваз фо Бе Кифег ргосеззе ю а44 сепеге4 еетеп 0 Ше 
ех1зНпе соПесйоп$. Переп4ше оп Ше зреслИс гедитетет оЁ Фе \таррег аррИсайноп, 1$ ргосез$ сап Ъе уа5Йу 
Ч1егеп:. ТБе Йо\усваг о? Пе а1еогит чзе4 ш ЗАРКО$ зегусе 15 зво\ми ш Етечте 2. 





















1= 1..М ув М 
Бешс соПесйоп 


зиррог$? 
А94 чап поде Вива зиррог$ ю з1ае 
АЗЧ епа поде 


ль 
— [39 
(9) 






Стеже 


— «= > 


Стеа{е по4е этопрз 


— = > 


уе$ 
Стеа{е го этопрз 







я Ведгам 
рр:у 
ргеззигез? > 


Арр/!у Югсез №0 подез 


<=> 


АД поде’$ 0$ о Фе по 
отопр 


АЗС год 1 


уез 


АЗС год’$ ГО ю Те этопр 


по 
[5 ат рошй 
(0,0)? 


А94 5ап рошЕ соогатае$ 


10 по4е$ 


Е1е. 2. Ро$(-ргосеззте Йо\усвай 
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Аз а геза[, 2 соПесйопз аге сгемеа: СоПеспопМ№оаеу ап4 СоПесйопКоау$. Те даа зюге ш Фезе соЙесйопз$ 
1$ заЦае юг Аиег таттршайоп ап4 са1сШайоп$. 

Вези!65. \Беп Фе тодше 15 ицеогайе4 шю а меб аррИсайоп Гог са1сшайп» 6\о-дитепзлопа| го4 згасагез, 
Фе сБотсе оЁ Ше буре апа за Буре оЁ Фе зсВете оссигз ш Фе отаршса[ изег ищегЁсе. То ещег $5 4итеп$101$ а5 
уе аз а4Ч1юопа| рагатает$, Не]4$ ап з\исВез аге изе4. Ш Ше туре оР \ги5$ ап рагатеегз аге свапоеа, Фе 
ргейплпагу 4ез1еп зсВете 15 геБий{ аютайсаПу. 

То Шазиже фе го]е оЁ Фе тоде эт ш Фе НатехогК оГа \еб аррйсаНоп, \уе м детопзгае Фе 4ез1еп 
апа са1сШайоп оГап шсПпеа згасеаге 20 т 1оп2, 10 т 128 у 3 т ой5е{ ап4 соп$15Япе оЁ5 зесйоп$ (рапе]$). ш Фе 
поез о? Фе иррег сВот4 ог {1$ згасвге, сопсепиже4 Гогсез оЁ 5 АМ аге аррПе4 уегасаПу, Гог Ще [еЁ апа г1>0{-10105 
по4ез, 15$ Югсе 15$ гедисе4 Бу 50 % ачющтайсаПу ут Фе РошЕ Гоаа Еасб®юг орНоп. Тве тодще пмегасе ап4 Фе 
гези эбгастаге ЧФлаотат аге зпо\ми ш Е1сиге 3. ТВе згас@ге 15 сася]=же4 уу Фе Вице е]етепе тефоа [18, 19]. ТВе 
гези т» Чласталт о ах1а1 Гогсез 15 зЗпо\уй ш Е1этге 4. 


Тпу55 Оезкупет 


} 


Бегкли, | (т): 


неюм, В (т) 


ОН5е4, Ч (т) 


4 
О х 
- } 


М№илтфег о! Рае 


Рояги: и0а45 оп Мофез 


Маше Ох (Км);  \аше Оу (КМ) 


54 Ве Рони 1020 24 Фе Епд$ О ТИи5$ 


Рояги (сад Расфог а{ Не Егкф5 оЁ {Не Тгиз$: 





15 5 5КМ 
м 
ы 5КМ т :. 
УК 
4 
ы Км т в 
5КМ 14 22 —— 
ы 2.5 км 321 то 
ь № 6.4 Г\ 
11 -—4 8 о | \ 
57... | 
2 6 —— | \ 
10 ее / \ 
{ 7 \ 
| 
о } \ ] \ 
| | 
8 / \ | 5 \24 2 
] \ 19 |131 р \ 
] \ м 14 30 18 \ 
в | й э [29 ] \ 
| \ д 2 \ 
5 3 [з ] \ 
| \ | \ 
. \ \ 
3 / \ | Е ПЕ 3 
\ [5 20 ‚$ 
2 —— 
о $ 
| \ в ‚> 
[ 5 ‚$ 





Аха! Рогсе С\здгат (№) 
° я И ПО О ПО ПОНИ 
+ х -23.608 -18,886 -14.165 -9.443 -4.722 0 3.672 7.345 11.017 14.69 18.362 


Е1с. 4. Ахла| Югсе Частат 
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Тре уа[аез ог зиррой геасйоп$ ап ищегпа| Гогсез оЁ Ще са|сшаеа згасвге аге ргеземе ш Та ез 1 апа 2, 
гезресиуе[у. ш ТаЫе 1, К 1$ Ще уегиса| сотропепЕ о Фе заррой геасйоп, Н 1$ Фе Вопхоша! сотропеп оЁ Ше 
5ирро! геасНоп, М 15 ве Беп4 ше тотепй. 

ТаЫе 1 


Бирро! геасНоп$ 
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015си$$10п апа Сопс$10п$. Озтшсо 11$ зоЙ\аге тодШе ут Фе Ёате\уогК оЁ Ве 1001 Юг са1сшШайипх год 
убгастагез аПо\уз Юг Фе зпирИНе4 ргосез$ оР тодейпе ап сайси]айпе сошрМех $$ зиасбиагез, 4ез1еп ите, апа 
гезоигсе тгедисйоп. ТБе зоЙ\аге тодще ргоу14ез 1001$ Гог зрес1Руше уаптоиз$ бурез оЁ згас@ге$, арр!уше 10а4$ ап4 
а551епте ргорегае$ оГа год зубет, \”мсВ таКез Ц а изе | шзйлитеп" Юг 4ез1оп епотшеегс. 
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